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VORTRAGSBERICHTE!) 


A. Mechanik 


_ Untersuchungen über die Stabilisierung eines Doppelpendels 
_ mit periodisch erschütterten Aufhängepunkt 


Pe 


Von Wolfhart Haacke in Braunschweig 


78 ® 7 u Wird ein physikalisches Pendel in seinem Aufhängepunkt vertikal periodisch erschüttert, 


Dr 
E 


3 so ergeben sich neben der natürlichen zusätzliche Gleichgewichtslagen, deren Zustandekommen 


Bevor Erdelyi.(Z. angew. Math. Mech. Bd. 14 [1934], S.235) und Klotter (Z. angew. Math. 
re Mech. Bd. 19 [1939], S.289 und Forschg. Ing.-Wes. Bd. 12 [1941], S.209) mathematisch auf- 
E= geklärt werden konnte. Die voneinander abweichenden 
3% Ergebnisse dieser beiden Autoren lassen sich aufGrund 
& der Diskussion der Bewegungsgleichung als Grenzfälle 
E erkennen, die ein ganzes Stabilitätsintervall von Gleich- 
gewichtslagen begrenzen. Welche Lage dann tatsäch- 
BE lich eingenommen wird, soll demnächst in dieser Zeit- 
schrift ausgeführt werden. Bei der Untersuchung der 


entsprechenden Erscheinungen beim Doppelpendel soll 

im folgenden der Ergde&lyische Grenzfall betrachtet 
werden. Die dabei verwendete Methode besteht in einer 
Ausdehnung der Erd&lyischen auf ein System von 
Differentialgleichungen unter Einschaltung der Stö- 
rungsrechnung?). 

In der Figur stelle der Streckenzug ABÜ ein 
Doppelpendel dar, dessen Aufhängepunkt 4A peri- 
 odisch in Richtung der n- Achse erschüttert wird 
— (NM=4sin2rnt). Wir fragen nach der Existenz 
kleiner periodischer Schwingungen um Mittellagen: 

=” x < t 


3 =, ‚seien die reduzierten Pendellängen, m, und m, die Massen und ©,, © die Trägheits- 
omente der Pendelarme um A und 5. Mit der Erdbeschleunigung g und der unabhängigen 
nderlichen u=2rnt ergeben sich nach den übligen Vernachlässigungen als Lagrange- 
. Gleichungen für kleine Schwingungen _ 


«608 («—ß)y +2 0942 0x, 608 («—B)y’ rs (h+ y’) cos a-pP— x, (h+y”)sin a] 


1) Ausführliche Berichte über die Vorträge von Prof. A. Walther, Darmstadt und Prof. D. R. Hartra, 
‘werden in.einem späteren Heft erschinen. ER SuN 

se Untersuchung hier gilt also mit den Einschränkungen, die wir in der weiteren Arbeit ausführen 

ische Frequenz und der gegenseitige Einfluß der Pendelarme wird davon nicht berührt, während 

ngigkeit der Mittellage (mindestens beim einfachen Pendel) sich nicht bestätigen läßt. 
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Dabei ist 
FR MaSal re Im sı + mein a Zeil Ba Tr EL 
I Op+ mal O4 Mal? O5 l Ann? 
und 20. der Reibungskoeffizient. N 
— —tga und y= —-tg ß ist Partikularlösung von (2). Wir brauchen daher nur noch 


das homogene System (2), zu betrachten. 
Um zunächst einen groben Einblick in die Verhältnisse zu bekommen, soll eine Störungs- 


rechnung mit dem Ansatz m = — durchgeführt werden: 
'2 
Pla) = golu) + m gi) H mega) + 
ya)=ylu) rm pa) Hm put ©. 


Für Q,, y, ergeben sich lineare inhomogene Differentialgleichungen mit periodischen Koeffi- 
zienten. Die homogenen Bestandteile sind für alle » gleich. Durch @, = e”** y (u), y= ee" o(u) 
ergibt sich : 


4" | +7 (h+ y")cosalx—=0, o' —|0? a (h + y’’) cos B\ o==08 
1 2 


Nach dem Floquetschen Theorem existieren Lösungen der Form x = e’% f, (uw), ® = e’“f,(u) 
Dabei sind o, r Konstante, f, und f, in 2x periodische Funktionen. Für die Existenz der Gleich- 
gewichtslagen ist notwendige Voraussetzung, daß Ro = NRr= oe ist. Dann ist aber v= =. 
l \ 2 
Daraus kann man schließen, daß cos = ST cos & ist. In erster Näherung hängt also der Winkel 
1 

zwischen beiden Pendelarmen nur vom Quotienten der reduzierten Pendellängen ab, nicht aber 
von der Erschütterungsfrequenz. Die Partikularlösungen der einzelnen inhomogenen Gleichungen 
lassen sich so bestimmen, daß (1) erfüllt ist. 

Jetzt soll (2), unter den Voraussetzungen m< 1, !,—=1,, also «= ß untersucht werden. 
Da die Erschütterung nur in der zweiten Ableitung auftritt, wird die Erschütterung n durch 
zwei Parabelbogen so approximiert, daß diese dieselben Nullstellen wie 7 haben und stetig mit 
stetiger Tangente ineinander übergehen. Die Höhe A der Parabel werde so gewählt, daß beide 


Schwingungsformen dieselbe kinetische Energie haben. Dann wird A=0,% A,. Ist T die 
Erschütterungsdauer, so wirdin O<t< n TR — Saar, undin 2 stsTt J= + . 


Setzen wir diesin (2), ein, so ergibt sich ein System mit stückweise konstanten Koeffizienten, das 
wir getrennt in beiden Intervallen lösen. Die acht Integrationskonstanten bestimmen wir so, daß 
die beiden Lösungsteile stetig mit stetiger Tangente ineinander übergehen und die Floquetsche 
Form haben. Die achtreihige Koeffizientendeterminante dieses linearen homogenen Gleichungs- 
systems zur Bestimmung der Integrationskonstanten muß notwendig verschwinden. Dies ergibt 
nach längerer Umformung die ‚„‚Eigenwertgleichung‘‘ zwischen & und n: 


118228 Ä 
ee 
of2re=cosv: cos w 3 ae „j sin v sin w (3) 
mit 
“ 8A 84 
v— —n? ann =) cos a, w = —n? rl + = cos & 
und 
@, ats timam, 
2 2(1— x,%;) F °. 


‚Die Existenz der kleinen Schwingungen beweist man mit Hilfe der Verallgemeinerung des 
Floquetschen Theorems auf Systeme und des Satzes, daß die charakteristische Gleichung bei 
geraden periodischen Koeffizienten reziprok ist. Letzteres kann man mit Hilfe einer unendlichen 
Determinante zeigen. Die Lösung der Eigenwertgleichung (3) wird wie bei Erdelyi bewältigt. 


PR 7 £ 3 € 
Es gilt lim «= > Da@7 #6, ist, erhalten wir aber zwei Lösungen, die ein Intervallabgrenzen, 
No e H 
in dem die Mittellage sich befindet. Zugleich ergibt sich so eine Möglichkeit zur Abschätzung, 
ob m klein genug ist um eine Stabilisierung zu ermöglichen, a N. 


Anmerkung: Das vorliegende Autoreferat ist zugleich Auszug aus meiner Diplomarbeit, die zur Entnahme 
von Einzelheiten bei der Technischen Hochschule Braunschweig eingesehen oder angefordert werden 
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Über eine neue Iterationsmethode zur Berechnung der Knicklast 
Von H. Jung in Stuttgart 


‚ Ein ausführlicher Bericht wird in einem späteren Heft dieser Zeitschrift erscheinen. 


Zur Berechnung schwach gewölbter kreissymmetrischer Platten 
Von H. Münz in Braunschweig 


An einem Breitenkreisschnitt einer Rotationsschale betrachten wir die Schnittkraftkom- 
ponenten P, parallel und P, senkrecht: zur Schalenachse, die entsprechenden Verschiebungen Q, 
und @,, das Biegemoment P, um die Breitenkreistangente und die Verdrehung Q, der Meridian- 
tangente. Das Variationsproblem vom Minimum der potentiellen Energie, geschrieben in den 

. Verschiebungen Q, und deren ersten Ableitungen nach dem Breitenkreisradius o, läßt sich durch 
' eine.Legendresche Transformation in das kanonische Variationsproblem 
& 3 
: || 22, n)ao=0 EEE EVEN SE en (1) 
o=1 
umformen. Entsprechend der analytischen Mechanik sind hier die Verschiebungen @, die ‚‚Koor- 
dinaten“ und die Schnittreaktionen P, die „Impulse“. Q, ist in diesem Variationsproblem zy- 
. . . . . . - * ’ 0 14 
- Klische Koordinate. Folglich liefern die sechs kanonischen Gleichungen Q; = ne und P, = — u. 
(s—=1,2,3) zwei Quadraturen für die Schnittkraft- und die Verschiebungskomponente parallel 
zur Schalenachse, P, und Q,, und nur noch ein System von vier linearen Differentialgleichungen 
für die restlichen vier Größen. Sehen wir von den von der Belastung herrührenden Gliedern ab, 
so hat dieses System die Gestalt 


d d/ı 1 s 

1 Par 5, 5 = gm Pı+tQN a VAR 
d ? d /1 | 

1 PP)= Plan Qs = P3}, 2G = Ds REN (3), 


worin die Koeffizienten ß, &,&, durch die Form und das Material der Platte gegebene mit o 
variabele Größen sind, und £(o) die Meridiankurve darstellt. Streichen wir zunächst die mit £’ 
behafteten Glieder, so zerfällt das Differentialgleichungssystem in zwei voneinander unabhängige 
Systeme von zwei Gleichungen, die je für sich kanonische Gestalt besitzen. Die Bedeutung dieser 
Systeme wird besonders deutlich, wenn die Platte so wenig gewölbt ist, daß die in den Faktoren 


&, & auftretende bezogene Bogenlänge Y1-+ £’? gleich eins gesetzt werden kann. Das verkürzte 
System (2) stellt dann die Scheibengleichungen dar — P, und Q, sind die in der Plattenebene 
radial wirkende Schnittkraft und Verschiebung — und das verkürzte System (3) die Gleichungen 
der Plattenbiegung mit dem radialen Biegemoment P,und der Verdrehung eines Plattenradius Q;. 

Die vorher vernachlässigten Glieder werden durch eine in der Himmelsmechanik gebräuch- 


- liche Störungsrechnung berücksichtigt. Von der Hamiltonschen Funktion HZ des kanonischen 
. Variationsproblems (1)1äßt sich nämlich die „‚Störungsfunktion“ @ = £’P,@, abspalten. Kennen 


wir nun ein kanonisch konjugiertes Hauptsystem von Grundlösungen des ‚‚ungestörten“ Problems, 
so erhalten die hier zweckmäßigen Lagrangeschen Störungsgleichungen kanonische Gestalt. 
Sorgt man dafür, daß auch für das vollständige Gleichungssystem (2), (3) ein kanonisch konju- 
giertes System von Lösungen bestimmt wird, so lassen sich die Partikularlösungen für die jeweilige 


Belastung leicht nach der Methode der Variation der Konstanten berechnen. 


Der eingespannte, achsial pulsierend belastete Stab als 
Stabilitätsproblem') 


Von F. Weidenhammer in Clausthal 


Es ist bekannt [1], daß die gestreckte Gleichgewichtslage achsial pulsierend belasteter Stäbe 
(Bild 1) für gewisse Längsanregungsfrequenzen und Amplituden instabil gegen eine Auslenkung in 
Querrichtung sein kann. Eine Anfangsstörung in Querrichtung kann dann zu unbegrenzt an- 
wachsenden Querbewegungen führen. Für den Stab mit frei aufgestützter Lagerung läßt sich 
diese Stabilitätsfrage aus den Eigenschaften der Mathieu-Gleichung heraus beantworten. Für 


1) Gekürzte Wiedergabe eines Teilproblems aus der Dissertation des Verfassers (Clausthal 1950). 
r Fr 7 P} £ ® 17* 
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den Stab mit eingespannten Enden muß man ein allgemeineres Ditferentialgleic 
entsprechend vermehrten Instabilitätsbereichen diskutieren. ar Bene Ye 
In der Differentialgleichung der Querschwingungen des gedrückten Er ve 


EI Wann — Erw tu Bw =0 ... ... Ey nö 
mit den Randbedingungen ® — w, — O und beliebigen Anfangsbedingungen, ist %, e: Lö- 


sung der Randwertaufgabe für die Längs-: “, 


En R+hcoswt schwingungen eine bekannte Funktion: 


Brr, .:B BE er ir En 
an. EHER: ee 
Bild 1 are Tore Yu 
Hierin braucht das Resonanzverhältnis 7 = ae 1 nicht verschwindend klein zu sein. In 0m 
0 Längs x 
ist ein Dämpfungsglied eingeführt, um wenigstens näherungsweise den Dämpfungseintluß zu er x 22 


B 


fassen. Nach der Transformation w —= exp © PER ): "Wa, {) kann man der Ditterentialgleichung 
für W ein Variationsproblem ZUORÜRER und dieses wie in [2] durch einen Reihenansatz mit zeit- 
abhängigen Koeffizienten W = > Xn(X): T „(t) genähert lösen. Mit vorgegebenen, die Rand- 


bedingungen erfüllenden Finiktionenike erhält man dann die Eulerschen Gleichungen zur 
Bestimmung der Zeitfunktionen: 


62T HT] — Ecoss- Sn, 0, a: BR Ne...) 
Ra worin die Frequenzverschiebung durch die Dämpfung vernachlässigt ist. In (2) sind die An- 
j ne emieng ö und die Biegeeigenfrequenzen 6, durch ® = ö er als dimensionslose Zah- It, 
f 5 len eingeführt. Ferner ist auch durch = RERe: die Anregungsamplitude und mit gs = ot die 
EULER : 


Zeitzählung dimensionsfrei. Die Kochen Fyn bilden eine Kopplungsmatrix. 
Da die Lösungen von (2) periodisch 2. Art sind, kann man sie in der Form 


Pr, 
r ap 
« 


E DE meisten. RE SE IE TE SU EN 
y=—on Seh Kan Bun r 
Kr ’ ansetzen. Zur Entscheidung der Stabilitätsfrage muß man die i. a. . komplexen ‚charak- 
; teristischen Exponenten o„ kennen. Da die Lösung im wesentlichen durch Se Faktor“ 


exp IBETT a Er: Si bestimmt ist, Er es ns sofern ee : En It ne \ 


Zahlen unterscheiden, so hat die Determinante A(e) nur einfache 


2rR„sin Pr TE 
Ale)=1 Ds 
Re ’ (e) = +3 cos Ang a RATE KR 
Br A(e) = 0 liefert dann N Wertepaare -+ On. Zur Preis der 


rechnung erforderlich: FR 


Me: € (=) 
20% RR: 
BR „Se 
; | RR 2, 2 qm la Be 
Auch Sonderfälle hinsichtlich der q, lassen sich auf ähnlie 
gehörige Werte von öund &, für welche die gestreckte . 
in einer ö-e-Ebene nach Art der Struttschen 


Nähe der we =—, »P Br 0 Ban 
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 Instabilitätsbereiche zweiter Art in der Nähe der Schwebungsfirequenzen öd, = Ir rd ‚ von 

Fe: welchen durch die Dämpfung jedoch sehr viel größere Teile ausgelöscht werden als von Bereichen 

erster Art. Insbesondere ist interessant, daß solche Bereiche zweiter Art für eine ‚„‚gerade und 

= ungerade“ Frequenz (d.h. für Frequenzen, welche ent- & 

E sprechenden Schwingungsformen zugeordnet sind), 'nur 

i dann entstehen, wenn in Erweiterung der bisherigen TE 

E Untersuchungen nicht mehr n=.0 verlangt wird. Dies 92 

ist damit gleichbedeutend, daß man sich nicht auf extrem 

‚schlanke Stäbe beschränkt. ja 
- Vergrößert man die Längsanregungsfrequenz bis 07 

- zur ersten Längseigenfrequenz, so tritt ein Zusammen- 

wachsen zweier Bereiche ein. In Bild 2 ist eine solche 

Erscheinung für einen Bereich erster Art angedeutet. Zur 4 

genauen zahlenmäßigen Festlegung der Stabilitätsgren- bo Ö, Res 

zen in diesem Falle braucht man eine asymptotische Bild 2 

Integration von (2) für e>. Da hierüber bisher keine 

Ergebnisse vorzuliegen scheinen, kann man diesen Fall zahlenmäßig noch nicht zufriedenstellend 


instabil 


B: verfolgen. 

AR Literatur 

er 1] E. Mettler: Mitt. Forsch.-Anst. Gutehoffnungshütte 8 (1940), S. 1—15, 

2] E.Mettler: Ing.-Arch. 17 (1949), 8. 418- 449, 

Bes 83] H. v. Koch: Acta math. 16 (1892/93), S. 217— 295. Vgl. auch J.O. Fleckenstein: Comment, math. 
e helv. 15 (1943), 8. 367— 376. 
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# Zur thermokinetischen Theorie der Elastomere 
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Von Walter Lode in Leverkusen 


Nach einer Aufzählung der typischen mechanischen, thermischen und dielektrischen Eigen- 
schaften der Elastomere (kautschuk-elastischen Stoffe), ihrer qualitativen Erklärung aus der 
Thermokinetik und der in der Literatur eingeschlagenen Wege zum Ausbau der Theorie wird über 
einen anderen, analytischen Weg zu demselben Ziele berichtet. Die physikalischen Voraussetzungen 
- sind mutatis mutandis identisch mit denen in einer inzwischen erschienenen Veröffentlichung des 
-  Vortragendent). Während darin aber nur zylindersymmetrische Formänderungen berücksichtigt 
NS ‚sind, werden jetzt beliebige Verformungen behandelt. 

Zu diesem Zwecke wird die Verformung beschrieben durch einen jedes infinitesimale 
Volumenelement des unverformten Körpers auf seine Gestalt nach der Verformung abbildenden 
Tensor?) A. Die Kräfte auf jedes Flächenelement (je Einheit seiner Fläche vor der Yezlormung) 
werden beschrieben durch einen Tensor K. 


Es ergibt sich, daß es außerdem vom relativen Volumen » = || (Determinante) ab- 
 hängigen inneren Druck p eine Fuktion w gibt, die nur von der Temperatur und von der 1. In- 


die Beziehung zwischen beiden Tensoren 
Er K=uw-4.—p:|4|-4 


autet., Unter "2 aus den. ee der Theorie begründeten‘ Annahme, daß die 


2 Wird ein Stab durch die Zugkraft Z je Elfcheneinheit. def eenen Querschnitts 
i een isotherm u so.ist. >; 


% 


| BE ha a. und ER, bei endlichen 
EMöth, Mech, 29 (1949), 8. 65— 75. 


alien. Flispier ea a 455. _ Flory and Rehner: J. chem, Physics 11 


rianten des Tensors A A, (dem ‚‚doppelten skalaren Produkt“ A..A,) abhängt und mit der 
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keit von der Länge /, dafür aber eine Abhängigkeit von der molekularen Struktur geliefert hat. 
Wird ein Rundstab mit dem Halbmesser r durch das Moment M um den Winkel ö je Längen- 
einheit isotherm verdrillt, so kann auch hieraus die Funktion 


0040) -(3 an, faar 


berechnet werden. 


Die Funktion wist also eine charakteristische Eigenschaft jedes Elastomers. Sie gestattet 
einerseits die experimentelle Ermittlung des Gültigkeitsbereiches der Theorie, z.B. durch Ver- 
gleich ihrer aus dem Zug- und aus dem Torsionsversuch berechneten Werte, andererseits inner- 
halb dieses Bereiches die Beherrschung des mechanischen Verhaltens unter beliebigen langsamen . 
Formänderungen, außerhalb des Bereiches die Erforschung der noch nicht von der Theorie 
berücksichtigten Vorgänge. 


V-Verzahnung mit sroßen Profilverschiebungen'!) 
Von M. Bergsträßer in Bingen/Rhein 


Die allgemeine V;-Verzahnung mit Werkzeugen nach DIN 780 z.B. mit dem Flanken- 
winkel &, — 20° führt bei großen Profilverschiebungen zu Zahnformen, die bezüglich der Biege- 
beanspruchung, der Flankenpressung und des Zahnflankenschlupfes günstiger liegen als die 

ohne Profilverschiebung hergestellten 
ü % Re Nullräder (Bild 1). 


ee & 
N N In den allgemeinen Paarungs- 


m 
DL A vergleichungen ?) 
SI a nr 
[ALTER an g ED \ - 


7 


»%=Ü(l+ B,);: 
dg 
CO (y = — : C08 Ay, 
Ay 


- 


— 10,025 (a, — 20°) 3) 


Ss 


treten nicht die Profilverschiebun- 

2,20, 2,=40, m=15 mm, &,= 20°, 54 = #/40. - gen und z, einzeln, sondern nur 
a) Paarung mit Qy=L#ka- 450 mm m(X+%,)=-472mm ihre Summe als Unbekannte auf. 
b) Paarung mit a,=470mm , X:Xy= 25:2, , Miyex)“ +242mm Diese Summe ist derart aufzuteilen, i 
daß die Profilverschiebung des Klein- 
rades x, möglichst groß bis zu ,—=1 


Bild 1, Zahnräder mit verschieden großen Profilverschiebungen 


“zu machen ist, so daß mit einem Entwurf entsprechend der Gleichung 2: & = 2: 2, Ban 


schon sehr günstige Zahnformen erhalten werden. Hierbei werden die Möglichkeiten außer- 
halb der V-X-Verzahnung ®), die mit konstantem Achsenabstand tınd verringerter Zähnezahl- 
summe des Nullgetriebes arbeitet, erschöpft, indem der Achsenabstand vergrößert und mit 
normalem Werkzeug (x, = 20°) Stumpfzähne hergestellt werden, die stärker belastbar sind als 


die nach den bisher üblichen Methoden berechneten Zahnräder mit kleinen Profilverschie- 
bungen. | | 


Die natürliche Grenze dieses Berechnungsverfahrens liegt darin, daßder Überdeckunggrade 
praktisch nicht kleiner als 1,1 sein darf und daß die Kopfdicke des Kleinrades nicht zu kein, 
etwa nicht kleiner als m/4 werden darf.. eg | 


4 


‘) Eine ausführliche Behandlung dieses Themas erscheint demnächst in der VDI-Zeitschrift. RESENE 
°) Bezeichnungen: a, = rg, + gg Achsenabstand = Summe der Teilkreishalbmesser, au + a, Achsen- 
abstand bei V-Verzahnung, m Modul, t Teilung, S+ Flankenspiel, B = [(tgaw — au) - (tg, -— alte 
Achsenverschiebungsfaktor, r, Grundkreishalbmesser. Br, N Pod 
*) Bessere und für die Berechnung bequemere Näherung als die in DIN 870 angegebene Näherung: 
nn 4 N t . a ee } nr 
B]»=Yyl+7B,=Yl+ 13B für &, = 20°. Heide ko; Be) Sr ii ; 3 
*) H.M. Hiersig: Wege zur Weiterentwicklung der Stirnradverzahnung. 
(1949), S. 559566, Be 
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‚Beitrag zur analytischen Behandlung des neh 


En 3 BE Von ©. S. Heck in Wiesbaden 


Es wird gezeigt, daß sich die durch Inversion der Kegelschnitte entstehenden Kurven als 
Koppelkurven des gegenläufigen Zwillingskurbelgetriebes darstellen lassen!). 
Gehen wir von einem Kegelschnitt aus mit der Gleichung 


\ RES. 
FE £ cos @ Yo 
oder, bezogen auf den Punkt K (Bild 1), 


r}(e® cost —1)— 2rı {[(1—e2)2* — ep] cos 91 + y* sin pi} + (p+ eat)? — ar? — yr? = 
so erhalten wir durch Inversion in bezug auf K, d.h. durch die Transformation 
Br. Er ST ER PR 2 -(k >0, sonst beliebig), 
1 
| die Kurve 


| il 
N el ee ee lee a IE ER rg BE ZU Br 
7 | . +ecos?y, —1=0. 


4 
| 
| 
| 


Bildi- Bi Bild 2 


sich aus den Beziehungen 

3% a@=AB+CH—2AB. OE-0529, .@= AB: + CB +2AR.CR, 

er Be 2.08 —a@=2AB, ‚CB: a Va OE: cos? 9, 

A Be RA EIERN 

AB: cos mm UnEE 23 00 — + lcos (di +6) 

. rule tan Whssnit time tarmin tan 
) a Br EN 


N] tme[n Er ) 


m)? — a? — 0 


(3): 


"Die Gleichung der Koppelkurve Din gegenläufigen Zwillingskurbelgetriebes (Bild 2) ergibt 
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Durch die Substitutionen 


1 ei 1 
n KIN EEE = — I 
1a (Pr Ee)} ar y*?] Z mit L (+ m)? — a2’ | 
B: a RO ae 
Tr. ale Mora Be == —/[L, I = gL, &? u Im For Pı 9458 e- E. 


geht Gl. e: in (2) über, womit die eingangs aufgestellte Behauptung bewiesen ist. 

Ist ein Kegelschnitt und ein Punkt K gegeben, so erhält man die Daten des Zwillingskurbel- 
getriebes, das die zu dem Kegelschnitt in bezug auf K inverse Kurve erzeugt, durch TE 
der Gleichungen (4) nach a, c, I, m.. Mit er der Eulsgroben 


1 pr | 
IR ER *\2 DE *_. Br un 
I-,le +y (p-+ear2,  f = le —1) @* + ep], TEL 7 
und Me, Re 3 
2 
M=P+g+7, neh 7; ler ap: HAfg 
wird dann a er. 
anti an 
Se x Z 
(tm: ll HR), t=m’=- zb). 
Ist umgekehrt ein Zwillingskurbelgetriebe gegeben, so ist die zur Koppelkurve von P in 
“ bezug auf K inverse Kurve ein Kegelschnitt mit . Daten * = 
N Br We k:x mit «—=aLyL ei RN = 


“= —(fl-en),  y*=kgl, Be 


wobei sich ZL,f,g aus den a (4) bzw. (3) ergeben. 

Der aufgezeigte Zusammenhang zwischen den Kegelschnitten und en Zwillingskurbel- 
getriebe kann für technische Anwendungen von Nutzen sein. Z.B. gestattet er, die Piercy- 
Profile!) sowie als Grenzfall von diesen die Joukowski-Profile mit Hilfe von Kurbelgetrieben u 
erzeugen. Ohne Schwierigkeit lassen sich die Hauptabmessungen eines Piercy-Profiles aus den 
Daten der zugehörigen Hyperbel berechnen sowie umgekehrt aus vier gegebenen Hauptabmes- x 
sungen des Profiles, z.B. der Tiefe, der maximalen Dicke und der Wölbung des Profiles sowie er 
i dem Winkel an der Hinterkante, die inverse Hyperbel bzw. das Sasende Kubi er: GE 
Be . mitteln. EN x RER u; N 


Et Bei einer technischen Aukaahe ergab sich ken Prabien: Ge, 

ke: belastung p(z) des Randes der Halbebene > 0; EeSUEneS wird di . Vers 
als Funktion von x. Bild la zeigt Belastung und Versch 
nicht der Verlauf für ot ©, sondern im Gebiet 


Die Beziehung zwischen ia und 'v(e) lautet: 
nalen Potentialfunktion D(x, y), so sind p(x) die Be 


va) = Er Duni 2 


2) N. A. V, Pierey, BR, Wr Pipor und 1, H. Dr 
Ser. 7, Bd. 24 (1937) Nr. 161, 8 a | 
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> Cegeben sind die Werte der ableitung einer Potentialfunktion am Rande der Halb- 
ebene; gesucht werden die Werte der Potentialfunktion für diesen Rand. Die Lösung kann natür- 
lich nur bis auf eine beliebige Zusatzfunktion a,-+ a,x gefunden werden, da die Ableitung nach y 
dieser Funktionen den Wert Null ergibt. 


Nach der Potentialtheorie lautet die Lösung: 


ae. 


= 4% 
=, rom 


Dieser Ausdruck kann für verschiedene Werte x nach dieser Formel gefunden werden. Der Wert «a 
"ist beliebig und entspricht im wesentlichen der Integrationskonstanten a,. 
ne wollen die Ermittlung für die Stelle &= x, vorführen. Hierzu zeichnen wir die Kurve 


pin ® = 


für Z=g logakfthimisch unendlich. Diese geringfügige Schwierigkeit beseitigen wir dadurch, 
daß wir das Integral z. B. wie ae in Teile zerlegen: 


über £ auf, siehe Bild 1b, und bestimmen die Fläche. Der Integrand wird zwar 


nn 


ee 2 RR let M PO — rain? 


dE 


F v(E) he 


er E 


- [ron —las4 In I 


Mit anderen Worten: Wir FE von %— 4 B ta dem Integranden die Funktion 


& 


— p(z,) In Sn hinzu, um das Unendlichwerden zu beseitigen. Die abgezogene Funktion 


integrieren wir dann analytisch, sie gibt das 
vierte Glied der angeschriebenen Summe; 
. dieses Bu mit, —e=n: 


a 2) In. as) L In "an 
; = = Sa Pe) | (Inn —Ina)dn 
= = a plain Inn —nlo— aln 1a 


ER = —2p(2) a. 
"Bild lc zeigt den Integranden für den 


- graphisch zu ermittelnden Anteil. Dieses zwar 
SE Ben Verfahren erwies sich noch zu, 


Pin] 
erot) In] 


nme wurde. Da das Ergebnis nur PS) 
is auf elanglose Konstante berechnet 
ur es ee nicht o,— 


| PEI-PKY-) | 
AH ST. 


h 


Die Reibung wird als so groß vorausgesetzt, daß nur die Potenzen der Reibungskonstanten ı mit 


‚ wurzel hat, erhält man mit den Abkürzungen. 
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werden. Dicht vor der Stelle &= x, hat der Integrand ein Vorzeichen, Eine dieser ‚Stelle 
das entgegengesetzte Vorzeichen. Es liegt nahe beim Integrieren (Summieren der Flächen- 


stücke p(&) ; ER 248) von beiden Seiten sich der Stelle &—= «, zu nähern. Dieses erreichen wir, 
la 


indem wir für die gewählte Stelle = x, die Dr ve) —p2, —H)miti<a bilden, 
also von den Werten p(£) links von x, die Werte p(£)rechts von x, abziehen. Dann brauchen wir 


7 Er 


nur noch das Integral a Pe) — pl —2)] a d& zu ermitteln. 


Bild 1d zeigt den Verlauf von p(E) — p (2 x, — E) und Bild le den Verlauf des Integranden. 
Dieses Verfahren hat sich praktisch sehr gut bewährt. Erst dachte ich, daß es ein Ei des 
Columbus sei, später sagte ich mir jedoch, daß Boitor alles dagewesen ist. 


Der Massenpunkt mit konstanter Horizontalgeschwindigkeit auf dem 
beiderseits unbegrenzten elastisch gebetteten Geleis mit Reibung 
proportinal der Vertikalgeschwindigkeit 
Von Konrad Ludwig in Hannover 

Für die Geleishebung Z gilt die Differentialgleichung 


art x _ Be 

— EI For! er =M om 
Die Horizontalgeschwindigkeit wird mit vo und der Abstand 2 — vi vom 1 Massenpunkte mit u 
bezeichnet. 

= de ah 2 BT. % 
— BJ 7 ae a a rk ö 
Der Exponentialansatz 
. De Ae"u f 
führt zu der charakteristischen Gleichung 
PR EL 
MC ITRESRIN 2 BI, N 


größten Exponenten berücksichtigt werden. $ E 


v. mv 


vaulihn se m 
Die 4. Wurzel der charakteristischen Gleichung ist 5; ee ; Rh Ban 
EL E x . - ER a j 2 x - ae “x - I? £ 


- 


Das Gewicht des Massenpunktes Se mit 20% bezeichnet. Die Geishebung hinter dem Massen- E 


punkte ist E 
ok ur 
= el BR (> 2 =) 
ER Be er ZE 


ee 


ats N Br 


und vor’dem ed 


ale jr“ 


Gi De 


jugiert Ko Wurzeln und ” kleiner a 


. 
Z. angew. Math. Mech. ’ € 
Bd. 30 Nr. 8/9 Aug./Sept. 1950 B. Elastomechanik 243 


die Geleishebung 
Yn „7 En 

R K Yy | 

a + u 


E I en 
; VD EJo mw | | 


IM 
2 


mv? 
EJ 


4 
u— 
{9} 


mau 
& ee ee © 
e mv A. i 
12c+ — 2 Yyı cos Ar u— ; reel- 
"+ 
Die Abhängigkeit zwischen den dimensionslosen Veränderlichen 


m?v* Ba re 2 > % 
IHJr und 3 wird durch die Kurve 3. Grades 


Ber m? v* er REGIE L mente 
2cm 9EJe mv !3 '9EJe 


dargestellt. Diese Kurve wird durch,ihre Asymptoten 


0 3 6 3 RR m2v* 
Be: Dem 7 Tee 


in drei Äste zerlegt. Der Ast mit rein imaginärer Reibungskon- 
stanten erreicht das Maximum k—=0, wenn die Horizontalgeschwin- 
digkeit ihren kritischen Wert annimmt. 


Bild 1 


Neuartise Behandlung der Poissonschen und der inhomogenen 
Bipotentialgleichung bei rechteckigen Bereichen mit Anwendung 
auf Probleme der Torsion und der Plattenbiegung 

Von Th. Schade in Aachen 
Die hier benutzte Methode bedient sich der von Grammel (Ing.-Arch. 1948, Heft 4) be- 
handelten erweiterten Kreisfunktionen, Funktionen der algebraischen Gleichung: 
2 +y"=1; z=cos(n)); y=sin(n)v. 


Mit der schiefwinkligen Transformation 2—= 0: cos (n)v, y= osin (n)v 
erhält man für die Potentialgleichung A®=0 mit dem. Ansatz 
®,—=e"F;(v) die allgemeine Lösung: A 


F;(e) = [sin? (n)v + cos? (n)v}’? Sol 


artgtg (n)v. 


E \ Hand frei aufliegend 

FR Bild 2 

Re % a N R A i ; 5 

Ber Im Grenzfall des Quadrats (n— @) ergeben sich die bekannten Potentialfunktionen ®— »* 
4 2@=&-riy), und die Bipotentialfunktionen (AAy = 0) y— 22? (A = 0). Zur weiteren Rech- 


7 4 nung wird das 0, v-Koordinatensystem beibehalten. Teilt man das Gebiet des Quadrates in vier 
Gebiete I, H, III, IV (Bild 1), so hat man bei den behandelten Torsions- und Plattenproblemen 


rer 


= 


= 


Se : 
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dafür zu sorgen, daß neben den Randbedingungen auch die Stetigkeitsbedingungen an den 
Diagonalen erfüllt sind. Unter Vernachlässigung der Stetigkeit gewisser höherer Ableitungen 
lassen sich dann durch Koeffizientenvergleich leicht entsprechende Kombinationen von Poly- 
nomen finden, die den Randbedingungen exakt genügen. Die Rechnung ist für die Torsion des 
quadratischen Querschnittes mit einem Fehler von 0,2% für tuax durchgeführt. FAR 


Für die gleichförmige Belastung der quadratischen Platte ergeben sich folgende Trajek- 
torienbilder. (Bild2 an allen Rändern frei aufliegende Platte, Bild 3 an allen Rändern ein- 
gespannte Platte.) Se ( ER: j 

1 


De 
2) gerechnet 
00 


.. 


Unter anderm ist auch das Spannungsbild für die Belastung p—= % 
worden, sowie unsymmetrische Randbedingungen. 


X) 


Bild 4 


Eine weitere Anwendung auf die unendlich ausgedehnte Scheibe mit am Innenrand sym- 
metrisch belasteten quadratischem Loch ergab das Trajektorienbild (Bild 4) (Dissertation von 
Montag, Techn. Hochschule Aachen). : ZEN HE AR” 


Die gleiche Methode ist auf die elastisch gelagerte Platte sowie auf Platten mit unvoll- 
kommener Randeinspannung und auf Rechteckplatten mit dem Seitenverhältnis bis 1:1,5 an-\ 
gewendet worden, um vor allem für den Bauingenieur mit verhältnismäßig wenig Rechnung 
brauchbare Werte zu erhalten. $ \ EEE UT Re RT 


u-+54 Be 


u; , “ r net 


Allgemeine Theorie elastischer Flächenträger ohne : £ Na 


Bernoullische Annahme 


Von W.Zerna in Hannover 


(Schalen) dagegen ergeben sich Zweifel!). Es soll 
werden, in der die Bernoullische Annahme durch 
wird. ER EE" 


Die Mittelfläche werde Aurch AsenSlansluse® 


Abstand eines Punktes von der Mittelfläche, so wire 


!) Vgl. z.B. E. Reißner, Amer, J. Math. Bd. 64, 

?) Im folgenden durchlaufen kleine griechische Buch 
Buchstaben die Zahlen 1, 2, 3. Partielle Differentiation ı 
einer Tensorgröße durch einen senkrechten Strich angezeij ar 
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I Der Maßtensor Yir des Ran; der Maßtensor aag und die Koeffizienten bap der zweiten 
2; Grundform der Mittelfläche sind ı nun alles A Größen. 


Es Der Verschiebungsvektor“ v* werde mittels v; = “(L- charakteristische Längenabmes- 
sung) ebenfalls dimensionslos gemacht. ; 


Mit «'® als kontravariantem Spannungstensor wird der Spannungsausdruck 
2 Fre a 


3 R ar 
BE 
Be: angeht Darin ist N das Kronecker Delta, ae 4: 
© . © Mit (1) ne Fe Igor ,.@= laas| werden folgende Schnittgräßen definiert 
re = +1/2 +1/2 
FE ex (n) . q : m. g r 
er: na — Nor do, en j „grT’ao, Ra Lt 0 (2) 
Ben ir ; T 172 \ . —1/2 
3 und ferner die Belastungsgrößen 
pe he —abghe, mon Lie, Beh? , eler (3). 
WORDS REN 
5 NE . +1/2 k 7 I3 +1/2 } 
- kt - /E \ NE Nas era nern (4), 
ei -, a —1/2 a —1/2 
Wird > EB 
er (0) PERS (0) (3) 
Se To SW ae neh, wm, a EEE 
Br - 
= RE gesetzt, so lassen sich die Gleichgewichtsbedingungen wie folgt darstellen. 
rn = DE SER ne, Aber —=0, Sn bt Ag +4p?—=0, 
=: mr2| — g® + | , S re], Sir = m’e|a Ze) DEE (5). 


a: . Zur Aufstellung des Elastizitätsgesetzes werden nun Näherungsannahmen eingeführt. Zu- 
Ba: nächst wird ee daß die in der Definition der Flächenträger enthaltene Voraussetzung 


BR en | 2% ee ea) 


4. 


| erfüllt ist, Weiterhin werden die Rn 


N ee) 
eingeführt nad die leicht einzusehenden Annahmen gemacht 
RE er REN TEE, jo) (n+2) - [in (n+2) "Bl | 
. e BE : Ze Zilk il» Vz. & < %o&| > 2 SLBR 


x ‚Schließlich wird für den Benalt der Schubspannungen r?° über die Dicke der Ausdruck 
5 rin gl 4 AO ERBEN? 15 2.2 (9) 


_ Randbedingungen an den Te Hin genr erfüllt. Der einfacheren Schreibweise 
h en ERDE daß die Se ee Akt, Ah‘ und deren Ablei- 


(2) ut (1) (1) 
x Ze Ay AR EN (10). 
ME. 


nd ©. ®) 0) ergibt sich I mit Nernendung von (10) das Elastizitätsgesetz 


Be f Gahm + tb Bam) Wan aan Ra ar bm Esbuv u! 
B N 
= = Gel Bw 12 - : 
TE AG, stm Mr — Un): | (11). 


\ 


ieh a E 
N r “er Y7 f « 
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Darin ist @ — Schubmodul, «= Querkontraktionsziffer, 


Befw — ar a + ara® + es ara, 


Geßw — 2 la bau 1 aonhPv 1 ap bw I a hPa De - ee beP DE 


—_ Baba, — Errw pf = nn { Be 


Das Elastizitätsgesetz (11) ergibt zusammen mit den GlkichgewichieBEBEI 5 6) das 
vollständige System von Differentialgleichungen. 


Über die elastisch-plastischen Formänderungen der Platte 
Von W. Swida in Karlsruhe 
Die bekannte Plattengleichung: 57 
_ P(&Y) . | 
FE Be a ee ee) 
gilt nur für den elastischen Zustand. Wird eine Platte elastisch-plastisch beansprucht, so 
' ist die Gl.(1) für den elastisch-plastischen 
Plattenbereich ® (Bild 1) durch eine neue 
Differentialgleichung zu ersetzen, die den Ein- 
fluß der plastischen Zone (gestrichelt) berück- 
sichtigt. 
Bei der Aufstellung dieser Gleichung 
wurden folgende Annahmen zugrunde gelegt: 
a) Die Platte besteht aus einem elastisch- 
ideal-plastischen Werkstoff. Dabei 
sind die Zug- und Druckstreckgrenzen > 
gleich. n 
b) Die Punkte einer zur Mittelebene senk- 
rechten Geraden liegen nach der 
Formänderung wieder auf einer Ge- 
raden und bilden eine Normale zur 
Biegefläche. 


Auf Grund dieser. Aa Ee unter 
Benützung der Plastizitätsbedingung von 
Mises-Hencky: 


ER DEE ER j 
erhält man die Grunddifferentialgleichung für den elastisch-plastischen Bereich in der Kor 


u ee 
& 53 3 A IR h = 
Eee 


Bild 1 


ae] ER le) en \ 
Er) FEN ER sr 


Außerdem besteht zwischen der Halbhöhe £ der elastischen Zone und der Durehbiegung. w 2 
folgender Zusammenhang: CE EZ 


0o,(1— u?) BEE. 
ef ee m Eee 


Dabei bedeutet o, die Fließgrenze und u die Poissonsche Konsta 


s=-fl—utı; Pen ie; RL 


> n Hr bis etwa 1944. 
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Der elastisch-plastische Zustand ist in irgend einem Plattenbereich dann vorhanden, wenn die 
Momente in diesem Bereich folgende Plastizitätsbedingung befriedigen: 


h? ee En De FE 2 
7 > JM: EM} MM, F3M,>o, FRE (EIG 
6 
Für die gleichmäßige Biegung mit Torsion ergibt sich: 
1 
a 2: E) [M, — aM) 22 + 2(1 +), -2y+ (My, —uM,)g2. 
N 


Die Halbhöhe £ der elastischen Zone ist in diesem Falle konstant. 


Eine näherungsweise Integration der Differentialgleichung y"= f(y) 
bei Druckbiegeproblemen 


Von Josef Weinhold in Holzgerlingen 


Die analytische Behandlung des auf Druck und Biegung beanspruchten ursprünglich 
geraden Stabes führt auf eine Differentialgleichung von der Art y'’—=f(y), wenn die Form- 
änderungen im elastisch-plastischen Bereich liegen. Unter den Voraussetzungen der technischen 
Biegungstheorie folgt für konstante Längsspannung P/F im Stab (s. Bild 1) in der Regel ein 
einsinnig gekrümmter Verlauf y'’—y, wie es Bild 2 schematisch zeigt. Jeder Lösung für die 
beispielsweise hier in Betracht gezogenen Grenzbedingungen 


ya=l)=y, Ylıazl)=09 .....20.2 2220. ll) 


entspricht eine mögliche Gleichgewichtsfigur der elastischen Linie „=y(z). Von den möglichen 
Gleichgewichtszuständen ist hier hauptsächlich der indifferente, ‚‚kritische‘“ Gleichgewichts- 
zustand von Interesse. Für diesen hat die zweckmäßigerweise als Unbekannte des Problems 
gewählte Stablänge ! einen Höchstwert, den es zu bestimmen gilt. 


Df= konst 


Hartmann 
Ein=Sehnen» ; 
Verf: 


Bild 1. Biegelinie schematisch Bild 2. Verlauf von y"'—y Bild 3. Biegelinie aus zwei Ästen 


Die exakte Lösung dieser Aufgabe ist — soweit bekannt — bisher nur für den ideal- 
plastischen Formänderungszustand und den Rechtecksquerschnitt, sowie für den Kreisringquer- 
schnitt durchgeführt worden!). In allen anderen Fällen ist man noch auf Näherungsverfahren 
angewiesen, deren an Hand der eben genannten exakten Lösungen zu ermittelnde Fehler erheblich 
sein können. Schärfere, zugleich mit einem tragbaren Aufwand zu bewältigende Näherungs- 


‘verfahren sind aus verschiedenen Gründen erwünscht. Es ist für die praktische Anwendung 


derzeit erforderlich, die exakten Werte mit einer Genauigkeit von Zehntelprozenten eingrenzen 
zu können. Man kann dann z. B. aus dem Vergleich zwischen Rechnung und Versuch schließen, 
ob und in welchem Ausmaß sich vereinfachende Annahmen der Theorie auswirken. 


Die bisher meist angewandten Näherungsverfahren von Ros-Brunner und Hartmann 
ersetzen ins Analytische übertragen den Verlauf y’’—y durch je ein Geradenstück, nach Bild 2. 
Eine erste Verbesserung liefert die Sehne 1 2. Die zugehörigen Näherungskurven für die elastische 
Linie sind sin-cos-Linien, mit dem Scheitel bei &—=1. Es liegt nahe, den Verlauf y'—y-zur 


> 1) Siehe L. v. Wüllerstorff, Beiträge zur Frage der Traglasten außermittig gedrückter Stäbe bei 


AR ideal-plastischer Arbeitslinie unter besonderer Berücksichtigung des Kreisringquerschnittes. Bericht des 
Institutes für Festigkeitslehre und Festigkeitsprüfung an der Deutschen Techn. Hochschule Brünn, Veröff. 


in der Reihe Deutsche Luftfahrtforschung, FB Nr. 1912. Dieser Bericht enthält auch die einschlägige Literatur, 
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Die naher Biegelinie besteht jetzt aus zwei sin-cos-Ästen, die bet Ka 3 Bild 3) es 
chend dem frei gewählten y, mit gleicher Ordinate und Tangente ineinander übergehen. Mit 
den schon erwähnten Grenzbedingungen (1) ergeben sich dann die nachstehenden Bestimmungs- 
gleichungen für die unbekannte Gleichgewichts-Stablänge 7: 


n 


Pf = konst 


1 v4 IE 
ya\? | a3 ee 
= —1+ — la oe )=0 (2), 
Y; ö W13 sın W13%g = 4 
= Wgg (L— %;) = arc cos, h REES (3) 
% 3 - 
n " n " - } 
mit, > Ya Hs und uw. = 24a (4). 
j Y3—Yı Ya Yare 
% Für gleichbleibende Werte %, und %y,, aber für verschiedene ER 
Bild 4. Einschränkung des Verlaufes Annahmen von %, sind diese Gleichungen an Hand des Zusam- 


BR menhanges y'’’—y durch Probieren jeweils nach ! aufzulösen 


und /nax als die ‚‚kritische‘‘ Länge aufzusuchen. Man erhält so — wie leicht einzusehen ist —einen 
zu kleinen Wert von /max: Wählt man andererseits statt der Sehnen in Bild 4 die Tangenten- 
stücke 13’ und 3’2 als Ersatz für den tatsächlichen Verlauf y'’’—y, dann folgt in gleicher Weise _ _ 
wie schon beschrieben ein zu großer Wert von Imax. Diese beiden Näherungswerte schließen 
den genauen Wert bei den durchgeführten Rechnungen mit einer Schranke von Zehntelprozenten 
ein. Zum Vergleich seien die Fehler der einzelnen Näherungsmethoden angegeben, wie sie bei der 
Berechnung der kritischen Länge eines Stabes von rechteckigem Querschnitt aus idealplastischem 
Werkstoff folgen, wenn die mittlere Druckspannung gleich der halben Fließspannung ist und die 
Endkräfte in der Entfernung von einem Zwölftel der Querschnittshöhe außermittig angreifen. 
Der genaue Wert der auf den Trägheitradius » bezogenen kritischen Stablänge ist dabei 


= m _ 1,0520 r YEjor. E ist der Elastizitätsmodul, or die Fließspannung.) 


Methode | Hartmann Eine Sehne Zwei Sehnen Exakt Zwei Tangt. Ro$-Brunner 
Fehler v.H.| -— 2,14 BB BETZ ER RR 3. +30 


. 


Variationsmethoden zur praktischen Lösung von ebenen und räumlichen a 
Spannungsproblemen'!) a Er re 
Von Udo Wegner in Heidelberg RR 2 


Bei vielen Differentialgleichungen der Physik gelingt es einfach ein vollständiges System 
von partikulären Lösungen anzugeben. Die Bestimmung der Konstanten bei linearer Super- 
position der Lösungen zwecks Erfüllung der vorgegebenen Randbedingungen macht in vielen 3 
Fällen unüberwindliche Schwierigkeiten. Die bisher verwendeten Methoden zur Anpassung an 
“die Randwerte beruhten darauf, entweder ein endliches, lineares Aggregat von 'Partikular- ? 
lösungen an einzelnen Punkten streng den Randbedingungen genügen zu lassen oder die Rand- 2 
bedingungen „im Mittel“ zu approximieren. Gewöhnlich verwandte man das quadratische WER: 
Mitte] dabei. Die erstgenannte Methode ist praktisch nur bei partiellen Differentialgleichungen 2 
2. Ordnung und bei Benutzung von sehr vielen Randpunkten brauchbar. Die zweite Methode 
entbehrt der physikalischen Rechtfertigung und kann sogar — ER bei unendlichen 
intervallen — zu physikalisch falschen Schlüssen führen 2). Wir werden im folgenden ze 
daß man den Randbedingungen bei den verschiedensten ebenen und räur i 
aufgaben völlig äquivalente Variationsprobleme zuordnen kann, wi 
sehr einfach zu erzielende Approximationen für die Randaufgaben 
sind nicht nur auf Elastizitätsprobleme beschränkt, sondern gelten a 
tische partielle Differentialgleichungen. Die Form der Gewinnung de 
äquivalenten Variationsproblems ist die folgende: Man ehe zur 


2 Eine ausführliche Darstellung erscheint im Tue 
°) Siehe E. Freundlich, E. Hopss U. Men 


a Siehe dazu U. We gner: Bemerkung zur Lösun 
Randverschiebungen bei Scheiben. (Z. angew. Math. Meoi3 
Arbeiten des Verfassers. 


Be re at ah re 2 
$ 2”, er 


Mn: Math. M fe 
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gleichung ein zugehöriges Va üblen auf (dies BR möglich sein !), transformiere es in 
‘ ein freies V.-P. durch Einführung Lagrangescher Parameter, ‚stelle für dieses Problem alle 
Eulerschen Gleichungen und natürlichen Randbedingungen auf, beschränke jetzt die Menge 
der Vergleichsfunktionen auf den Bereich der Funktionen, die allen Eulerschen Gleichungen 
und natürlichen Randbedingungen genügen bis auf die Relationen, die den vorgegebenen Rand- 
bedingungen gleichkommen. Die Lösung des Ausgangsproblems liegt in dieser Menge. Elimi- 
nation der Parameter liefert ein im allgemeinen den Randbedingungen äquivalentes V.-P., wobei 
-— bedingt durch eine Legendre-Transformation — gewöhnlich nur eine Komponente des Diffe- 
= .. rentialoperators der Ausgangsdifferentialgleichung als Operator für die Differentialgleichung, 
die als Nebenbedingung für die zugehörigen Vergleichsfunktionen auftritt, dient. 
Beispiele: 1, 8=ebener, einfach zusammenhängender, endlicher Bereich mit stück- 
weise analytischer Randkurve EC: x(s), y(s). F= Airysche Spannungsfunktion. Spannungs- 


Een: AAF=0, mit Fle= f(s), |. ER). 
Äquivalent damit: 1/2 J (AF)?xdy = Extr. mit den en REN Äqui- 
valentes V.-P. für die Randbedingungen: 


1 [away + [(r2?-0:9)0= Br. 
® 1% ‚ 


und der Nebenbedingung 
A9—0. Dabei ist: ee 


; Praktische Auswertung: Ansatz: = ot, > + 34 pr, mit 9% 2 ir —b (2-4 in ne —in®B 


regulär analytische Funktion. c,, c; und £ Merden farch I. bestimmt. Bei der Auswahl der 9; 
und 9, beachte man die Symmetrien des Bereiches und der Belastungen und vermeide im Ansatz 
br nicht vorhandene Symmetrien des Problems. Man wird im allgemeinen 9, + ip, = et — 
- (2 + iy)® bzw. rke'*® wählen. Zur Bestimmung der Spannungen und Verschiebungen wird außer 


0) oF R : > 
p nur— und Er benötigt. Falls 9 + 9, = 2* ist, sei 


u 0x 5 

Er gr: a RR) EP 0% - x r 
3 en (ten) m, 
= RE Ten 
E: gesetzt. Dann mache man den Ansatz: 

2 

4 oH ERS 

nt Zat Dh, 

: > > Be , } oH son en IN 2 A 

Be | = | a u ee k Pk 

4 BE ? - ze 24 DEREN EN k=l. 


ud bestimme die d und d aus dem V. -P.: 


= +ler-]3 OH Like) — Po) ds = Exır. 


Die völlige Äquivalenz des V.-P. I.) mit der Ausgangsaufgabe folgt aus der Identität 


: ale ninay= z || amreayt raa+ —g° e)ds. 


ER I.) eignet sich auch bei unendlichen Bereichen; insbesondere zur approxi- 
Mask Bestimmung der Kerbspannungen. Konvergenzfragen des über unendliche viele Terme 
ausgedehnten Potentialansatzes für 9=0„-+-0o, können sofort hier funktionentheoretisch be- 
=... antw ortet werden. Konveregenzgüte sowie Kerbspannungsfaktor werden allein durch die Koeffi- 
en ee des ‚Äußeren des, nnentskreises auf das Äußere der Kerbe 


te _, : ro 


a ma 
er. 2 RT. 18 


A 


e 


nd 
r 


. 
> 
“ 
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mit 
uslc = u,(s) tangentielle 
_ und 
__ 772 | normaleVerschie- ey ‚08 
ul = ls) bungskomponenten. ° ds "ds “ 
u 00 du 
= — — j 1 m — ti 
D e _- 3 Dilatation R= eo —j. w. Rotation. 


Äquivalent mit dem Problem ist das folgende aanıı piohiem: re a 
le R) de dy = Extr. eu 


Das äquivalente V.-P. für die © Rand ingungen lautet hier: 


ll 


v 


)as — Extremum, 


mit D+ = = HUc-+ ıy) = regule analytische F 'unktion in ®. — Ansatz: : Sc :J 
Dar Santa z 3 

\ k=1 ni k=1 v> v+ x 

1—v 5 ES nn | >, ; ee 

3 R=o+ & CE Pk Sr . \ 4 

Bestimmung der c und © aus II. Bestimmung von wund vaus Dund R mit Hilfe eines Variations- = 
Problems siehe ZAMM.-Arbeit (loc. cit. 3). — Die Äquivalenz ergibt sich aus der Identität: Be 
VE A A Zw — Rje|deay x er : ü SE F 


if | @ 2 + = u) 5 E 


le emien, I E aa) . > E 4 


3. Gemischtes Problem: “ EEG, Es seien in & (ai) 0) die Spannungskom- 
ee X(s), R(s) gegeben, deiais 


aF 


t i : Arm 


2 =i ee 2er 40) ne 0 ER Er 


2 ö 
= x) dl; 3; ie ae 
ER e er Ys(s)) seien Wa — wo); Ufo = = U .() damit auch 2.uß) ) gegeben. ae 
Es seien x BET: 
Ei : RN 
er ot Dan | ER N ee ee 
Wir machen den Ansatz: FEN EN EEE IT GE RZ rt 


n L N. - \ 
Pecntzants ah und zu0r 


oF dx aFdy 
Flo: (ee reger ana we 


+AK)+B: + Ja a 2 
aF z 


es (1—») Dj er % 5 W+R: A 


Aa 0 > Te SEE 
EM ERE EI a TE Er 
Ce) * x 3” 

Eee j 
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oF 

(S, ist die Bogenlänge von@,). A, B, T,H sind unbekannte Konstanten, so daß Fe» und 
n |E, 
Nr FG; : Re» ter ee ‚6,, A, B, F, Habhängen. Zur Bestimmung der ee 


stanten c, € dient das Variationsproblem III: 


al Ierseus+ [p-38-0:0Ja+ {nat n)ar= stenum. 


Die c und 6 ergeben sich aus einem linearen, inhomogenen Gleichungssystem mit von Null ver- 
schiedener Determinante als lineare Funktionen von A, B,T,H. Setzt man diese Werte in f* 
und g* ein, so erhält man 4 reguläre, lineare, inhomogene Gleichungen zur Bestimmung von 


A,B,fT,H aus: 
f0) = *($ + 5.) 
681) = F*(S}) 
g(0) = g*(Sı + ©) 
951) = 9*(S1) 


Da F und OF/on nunmehr am ganzen Rand bekannt sind, verfährt man weiter wie in 1. 
4. Räumliches ee 


0D 
IS 
Aut 5, 
1 oD ou 9% ow 
NSTSTETTERN u PLABFTTAEF FE 
1 oD 


E Au a > 


. Randbedingungen: u/,—f(®); v/, = 9(w); w/,— h(o). 
® — räumlicher, einfach-zusammenhängender, endlicher Bereich, F =stückweise analytische 
Oberfläche von 8. — Mit der Aufgabe äquivalent: 


Sf | | Ip: = en (d% — Uy)? + (U, — w,)’ + (w, — %)”| dV/—= Extremum, 
8 


mit den obigen Randbedingungen. 
Äquivalentes V.-P. für die Randbedingungen: 


IV: ern rer mar 


+ Jaolı Nolan Er nyH un 


1— 2» de 1—% dy a 
EEE! a ._ — —— — Ext 
+ ho) H In) G 77 DD BR 
mit: 
2(1—»)0D OK o©H 232 9)0D2 00 DE 
1— 2» 95 0oy 9’ 1-2, 24 02.02’ 
2(1—v)8D 02H :G &G. OH. oR, 
EEE ne a oe 


Die allgemeine Lösung dieses Systems kann durch zwei räumliche Potentialfunktionen p,(z, Y,2); 
 P2(@, y,2) dargestellt werden mit: 


2(1—») dp ir " .0p,. ...2(1—») 97, Pr 
G= 1—2v dy i ee = 1— 2» Er 


gr 16: - P, (cos 6) ar ER « 1") cos md +” sin md) - Pl” (sin 9) ; 


ET & P, (cos ©) +. 23 IC cos m d-+ dm sin md)» P() (cos 9), 
18* 
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wenn z=rcosd®sin GO, 0< 0 <n, — P,„= Kugelfunktion x!% Grades, EeE 2 Ben. 
y=rsndsinO, —nSdSse — Pl” — Kugelfunktion x'%® Grades, m*® Ordnung 
z=1rc0s 9, r>0, 22 3 
ist. = Fa ee: ka e 
Die Aquivalenz vom V.-P. IV mit dem Ausgangsproblem folgt wieder aus einer analogen 
Identität wie die in 2). Aus p, und 9», können u, v, w wiederum aus einem Extremumsproblem 
gefunden werden. Es sei z.B. u= S+ P, wobei P ein Partikulärintegral von 


0° P, 
an — 0002 


ist, das leicht aus einzelnen Termen, aus denen p, besteht, aufgebaut werden kann. $ Deine 
sich als Potentialfunktion (Ansatz wie bei p, mit Koeffizienten e und ®) aus dem V.-P. IV’: 


IL + I + hleeuart [feo3Sao- worum. 


mit 2[o] = ulo] — P[o]. | R = Fr r 


Über die Berechnung der Gitterenergie endlicher Kristallstücke Be: 2 
Von OÖ. Emersleben in Berlin SE 


Bei einem Stück Kristallgitter vom Steinsalztyp, das in Richtung der rechtwinkligen 

* Koordinatenachsen abwechselnd gleiche positive und negative Ladungen aufweist, nenne ich den 

kürzesten Abstand zweier Ladungen die Gitterkonstante. Für unsere Zwecke genügt es, ein 

„Physikalisches Einheitsgitter‘‘“ zu betrachten, bei dem die Beträge der Ladungen und die Gitter- Kr 

konstante je gleich 1 sind. Da das Newtonsche Potential proportional der reziproken Ent- ° 

5 fernung ist, ist in diesem Fall das Selbstpotential eines Gitterstücks oder die Gitterenergie des N 
>. Kristalls die endliche Summe der mit abwechselnden Vorzeichen zu nehmenden reziproken 

{ Abstände aller Gitterpunkte voneinander. Um bei der meist sehr großen Zahl von Gliedern zum 

Erfolg zu kommen, hat man nach dem Vorgang von Madelung und Born das neutrale 5 

Kristallgitter allseitig unendlich ausgedehnt gedacht und den Wert der Gitterenergie pro e 

Ionenpaar berechnet. Den Betrag dieser Größe nennt man heute de Madefung-Kon- u 

stante des betr. Gitters. Sie sei A. Sei N die Anzahl der Ionenpaare in einem bestimmten TR 34 

2 Er - E- 


Gitterstück und , „d el Gitterenergie; so ist: 


N 


2N j 
TEE RE any Ach 
2 en 
falls der Grenzwert vorhanden ist. Dies ist nur der Fall, wenn man sn gleichartige Gitter- et = 
stücke zur Konkurrenz zuläßt. Durch eine interessante Anwendung der Funktionentheorie, it 
Hilfe der analytischen Fortsetzung, habe ich gezeigt, daß sich die Madelung- -Konstante 2 
. durch Werte einer Epstein schen Zetafunktion darstellen Jäßt, im Fall des Le I 
VE 
A=—Z ı 3 ıtqa). 
2.2 a 


an anerkannt. Dies darf uns ec nicht Sbachre tree ge en, 
der Oberfläche und der Größe des endlichen Gitterstücks zu erfassen. In dieser I 


in den letzten Jahren ein gut Stück weitergekommen. rasen 2 
_ Daneben können wir jetzt die Gitteremer gie, d ho BR pol: 
berechnen. . IE 32 


- Vollständig klar sind die Verhältnisse bei der Gi tt erg 
lineare Anordnung von n ee erhält. man 


Won lg+gttt.. 


! 


1) 0. En mersleben, .Das Selbstpotential der 
Bichten 3 (1950). Im Druck. — Insbes. Formel war 
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“ 


N i Hier tritt die endliche har Reihe von n Gliedern auf, „Z, bei der das erste Glied fehlt, 
. oder, ‚nach fallenden Potenzen von n als Eman. Reihe 


ee. 1 1 “ 
) Bi NVd— —(1— —t—— tt... 
> ® nlgn ( On +z BUS LSTOTE (3a), 
‘wo Ü die Eulersche Konstante ist, also 
A . (3b). 
EEE Be einmebenen Gitter hängen die Konstanten außer von der Art des Punktgitters auch von 
der äußeren Gestalt ab. Für ein quadratisches Stück eines quadratischen Punktgitters, bestehend 
: ; aus N BE Sitfergeraden mit ie n Gitterpunkten, erhält man 
BE UNSERER RT ld) 
mit 
; = > (Ar Sin 1 N ABO re lien 2... (da), 


Hierzu gelangt man, indem man die endliche Summe durch das Integral ersetzt und den Fehler 
abschätzt. Dabei zeigt sich, daß das zweite Glied eine Größenordnung in n niedriger ist als. das 
\ Haupsgbed! Setzt man die Gesamtzahl der Gitterpunkte n?—= N, so erhält man: 


Er a er 
wobei genähert 
E | Va Er ELEND IEERE Stra 9 (Ba) 
Er ‚Für Stücke eines kubischen Gitters ergibt sich . 
a EB Ne ee SS TR AUESETB): 


- Gestalt ab, z.B.?) bei einem kugelförmigen Stück kubischen Gitters 
.6 8 = i 
= — — EEE A te er 6 E 
| U y 6 0,967 (6a). 


Eine besondere Lücke schließt die Berechnung der Gitterenergie Dendlicher 

Stücke neutraler Gitter, bei denen ich mich auf rein alternierende Gitter gleich 
großer Ladungen beschränken werde. 

% "DiealternierendeGittergerade vonn Tönenpkaren führt, wie ihr homöopolares 

Gegenstück, zu einer einfachen Se Darstellung ?): 


= nli- +34) RER 9, 


i 


x N 
\ wo, gr: AN Aterwetende Harnnnisehe Reihe mit FE Gliedern sei, oder 


> a au tue m 


Br Für das schachbrettartig ahercierendö Gitter kann man rein numerisch für ars Anzahlen 
n\ ‚der ee an Nuadıatkente die DEN berechnen und mit dem Wert derMade- 


NE Tee? R STTERBN g REEL 
RR NR 3 EAN) 1) 1,6108 | at nn... (8a) 
; a Du 3] : 
n erhält so, ah für größere n genähert ER 


Du " i ar 
; ee a 


ARE. Dar. +0,37... —0,08/n +. 2) 


el aus leichschweren Massenpunkten, die sich 
1, Beitr., ‚Geophysik 61, (1,950), S. 163—1,80. 
(Boa). - 


tterene 'gie eines neutralen ebenen, insbesondere alter- x 


—609, insbes. 5. 594. 


wenn N wiederum die Gesamtzahl der Gitterpunkte ist. Die Konstante U hängt von der äußeren 


, 
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oder mit 2n?=N 


P 
g 
4 


BB—=—AN+ BUN H Bet} u ee 


wo die B, weitere Konstanten sind, mit vorläufigen numerischen Werten 


BB, 0,31... 0,1 Bi, = 017 Eee: 
Im Dreidimensionalen erhält man für das Gitter vom Steinsalztyp = 
BO— AN BNEBL EN N. N or 


vermutlich nach drittel Potenzen von N abfallend. 

Die Madelung-Konstante des Steinsalzgitters nach (2) hat den früher von mir be- 
rechneten Wert A = 1,74756..., der heute, meist auf 4. Stellen abgekürzt, weitgehend auch in 
die einschlägigen Lehrbücher eingegangen ist. Da das zweite Glied mit N2!3 relativ nahe an der 
Größenordnung des ersten liegt, erkennt man, daß in vielen Fällen dieses bisher nicht berück- 
sichtigte Glied keineswegs zu vernachlässigen ist. 


Anwendung der zyklographischen Abbildung in der Bruchtheorie 
Von (. Torre in Wien 


Nachfolgend werden wir den Inhalt des Vortrages nur andeuten, da die ausführlichen Ver- 
öffentlichungen des Verfassers®””) schon vorliegen. 

Die zyklographische Abbildung ist ein Zweig der konstruierenden Geometrie, deren Trans- 
formationsvorschrift darin besteht, daß man einem Punkt im Raume einen orientierten Kreis m 
— der Zykel genannt wird — in der Bildebene x zuordnet (die zugehörigen Abbildungen s.a. a. O.5) 
S. 326). Der Mittelpunkt des Zykels fällt dabei mit dem Normalriß M’ von M in x zusammen. 
Der Halbmesser des Kreises ist gleich der Kote M’M = z des Punktes M über x. Wenn man 
eine Gerade t des Raumes zyklographisch abbildet, dann besteht ihr zyklographisches Bild aus 
zwei Hüllgeraden i,, i, einer linearen Zykelreihe. Hiermit liegt eine Strahltransformation t—t,t, 
vor, die zweizweideutig und quadratisch ist. Ähnlich wird eine Kurve Kim Raume zyklographisch 
in der Hüllkurve K,K,= K, einer allgemeinen Zykelreihe abgebildet #), 2). 

Die Anwendung der zyklographischen Abbildung auf die Probleme der Festigkeitslehre . 
ist naheliegend, weil man die Spannungs- und Verzerrungszustände durch die Mohrschen?) 
Kreise darstellen kann. Eine wichtige Vereinfachung besteht darin, daß die Kurve K in der 
Ebene liegt. Der Halbmesser des Zykels ist dem!Halbmesser des Spannungskreises über (0, —0;) 
gleich, der die Größe der Hauptschubspannung 


Kir nn ns Shin aha Ba ra Fe 1 Me ed a nd rn a a ne LP 


ee 


(1 — 05) 
T, s .d) 
hat. Die Abszisse des Kreismittelpunktes entspricht der mittleren Spannung 
o—+0 
a ne 


Die Koordinaten der Hüllkurve bedeuten (im Druckgebiet mit guter Näherung) die Normal- o 
und Schubspannung r in der Bruchfläche. Ist eine Grenzkurve im Koordinatensystem (Om, T;) 


To TO re Ba ee eh BG) 
gegeben, dann erhält man durch die Transformation, die sich aus der graphischen Darstellung 
des Spannungszustandes ergibt, die Mohrsche®) Hüllkurve 


TEI) ER TE re et: 
Die Pknlattinken Ausdrücke für diese Transformation ereeben sich aus einem Hilfssatz über 


x E. Müller: Lehrbuch der Darstellenden Geometrie. Bd.I und IT. 3. Aufl. Verlag Teubner 1920. 

?®) E. Müller: Vorlesungen über Darstellende Geometrie. Bd. II: Zyklographie. Herausg. von J. L. 
Krames. Verlag Deuticke 1929. 

®») O. Mohr: Z. VDI 44 (1900) S.1524; Abhandlung aus dem Gebiete der Techn. Mechanik, 3. Aufl., ei 
5.192. Verlag Ernst & Sohn 1928. 

*) C. Torre: Österr. Ing.-Arch. I (1946), S. 36. 

e) C. Torre: Österr. Ing.-Arch. 1 (1946), S. 316. { A X. 

°) C. Torre: Der Spannungszustand in einem schweren Erdkörpor. Be der Akad. der Wiss, ER 
in Wien. Mathem,-naturw. Klasse; Abt. IIa, 156. Bd., 9. u. 10. Heft (1947), 8.583. er) 
% e FEDER: Schweiz. Arch. angew. Wiss. u. Technik 15 (1949), 8.116 u. 128 tern. Ing. Bi N a 


- Form erfaßt werden: 


2 wer. k 9 Me. ” N JE; En Br 
x h ee - 
ne wi, I - ar \ 


Fr 
Pr 
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die Gleichheit von Subnormalen der Kurven Gl. (3) und (4), die aufeinander zyklographisch 
abgebildet sind: 


Fr tree mit , —& ER REO 
m 
Mit Hilfe der Transformationsgleichung (1) und (2) geht Gl. (3) in 
BO ee ee (6) 
über. Aus Gl. (1) ergibt sich noch die Hilfsgröße 
,. (4 —1) 2 ‚do; 
7; = ee 1)" wo = 2, SEE Re (7). 


Aus Gl. (5) und (7) erhalten wir die Ausdrücke für die Berechnung der Normal- o und Schub- 
spannung 7 in der Bruchfläche, sowie des Bruchwinkels a, wenn Gl. (6) entweder experimentell 
oder analytisch gegeben ist: 


+00 a 70. a 2 
ira Te + er Ye;  tg=+Yyo .... (a,b; 9). 
Aus Gl. (9) ergibt sich noch die Differentialgleichung des Gleitlinienfeldes, die wir hier 
ohne Ableitungen angeben wollen (siehe ®, ?) 2; 
nal. en 7. 
a Ren 
Dabei bezieht sich y/(x) auf die Hauptnormalspannungslinie o, und o/(2) muß als Funktion 
von x bekannt sein. 


C. Strömungslehre 


Modellversuche mit wasserdurchströmten Zyklonabscheidern zur 
Vorausbestimmung der Abscheideleistung') 


Von W. Barth in Nonnenhorn (Bodensee) und K. Trunz in Karlsruhe 


Entstaubungsvorgänge verlaufen ähnlich, wenn neben der Forderung der Modellähnlich- 
keit die Reynoldssche und die Froudesche Kennzahl, sowie der Faktor y;/y; für Modell und 
Großausführung gleich ist?). Auf diese Forderungen kann bezüglich des Staubkornes verzichtet 
werden, wenn nachstehende aus den Bewegungsgleichungen für Staubteilchen abgeleitete Ähn- 
lichkeitsbedingungen erfüllt sind: 

” 
—— SEE EEE RE RE | 
5 const (1) 
2-8 
— — ee, EP 
RR cons (2) 
Diese Ersatzbedingungen gelten unter folgenden Voraussetzungen: 
a) Der Widerstand W der Staubteilchen kann durch eine Beziehung in der nachstehenden 


ur 


a er fu / 
ne 7ER Re = (3)- 


b) Der Wert Zr ist annähernd gleich 1. 
k 
Bezeichnungen: 
.. Fallgeschwindigkeit des Staubteilchens in ruhendem Medium, 
eine charakteristische Geschwindigkeit des strömenden Mediums, 
eine charakteristische Länge, 
Erdbeschleunigung, 
Durchmesser der Staubteilchen, 
Relative Anströmgeschwindigkeit des Staubteilchens, 
3) Vgl. auch die Ergebnisse einer demnächst erscheinenden Arbeit von K. Trunz, Karlsruhe, 
2) W. Barth: Die Anwendung des Modellversuches: zur Lösung strömungstechnischer Aufgaben. 
Z. 'VDI Bd. 92 (1950), S. 105— 110. 


son nes 
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y;  Wichte des Staubteilchens, IE 
yı  Wichte des Trägermediums, AN 
v Kinematische Zähigkeit, R der 
&, K Konstante, wobei O< K<]1. 


Die Durchführung von Modellversuchen zur Untersuchung von Entstaubungsvorgängen g 
scheitert in der Regel an der Forderung, daß entsprechend den Modellgesetzen mit einem sehr 
feinen Staub gearbeitet werden muß. Da bei diesem Staub einwandfreie Korngrößenbestim- 
mungen kaum möglich sind, außerdem der Staub an den Wandungen haftet und zur Flocken- 
bildung neigt, sind genaue Messungen meist nicht durchführbar. Diese Nachteile lassen sich _ 

. vermeiden, wenn an Stelle von Luft Wasser verwendet wird. Infolge der größeren Zähigkeit 
des Wassers kann bei gleicher Fallgeschwindigkeit- mit sehr viel gröberen Körnern gearbeitet 
werden. Außerdem ergibt sich der Vorteil, daß im Gegensatz zu Luft bei gleicher Modellgröße 
die Reynoldssche und die Froudesche Ähnlichkeitsforderungen leichter erfüllt werden können. 
Allerdings wird die Auswertung dadurch erschwert, daß Voraussetzung b) nur annähernd erfüllt = 
ist und der Faktor % nicht gleich 1 angenommen werden kann. Es ]äßt sich jedoch zeigen, daß in 
der Mehrzahl der Fälle praktisch durch diese Abweichungen keine größeren Fehler entstehen. 

Der in Bild 1 dargestellte Modellzyklon wurde mit verschiedenen Wassermengen und ver- 
schiedenen Fraktionen von Karborund (SiC) beaufschlagt. Die Abscheidegrade wurden in Bild 2 

in Abhängigkeit von der Durchflußmenge !V aufgetragen. ‚Der Ab- 

scheidegrad beträgt bei Durchfluß null 100%, er fällt zunächst mit 
steigender Belastung ab, um unter dem Einfluß der Zentrifugalkräfte 
allmählich wieder anzusteigen. 


Material: SiC I 7 DE SeRS u 
 Hörnungs-Nr 180 150 a ET 
mIem/s] 086. ın m | 


Abscheidegrad 


OEEEIT FRITZ ER NE 
Wassermenge ın l/S WERTEHS a ar 


Bild ı ee IE RE 


Da die Verhaltnis auf der Zylinderfläche F, (Abb. 1) für die Ähschsideläistäne Sutschielz EN, 
dend.sind, wurden die wichtigsten Größen auf diese Fläche bezogen. Bezeichnet man mit „de ä 
"relative Anströmgeschwindigkeit eines im Gleichgewichtszustand befindlichen mit der Umfangs- 
geschwindigkeit u auf der Zylinderfläche F, kreisenden Staubteilchens von der ale Be 
digkeit w;, so kann diese Größe wie folgt errechnet werden: ne ER je a En 


1: ER sn PA} > er ER, RR 
BR Far Br 5 BE 

17 CE ri d, gi Be DZ = ER De ve 
Bildet man den Quotienten w,/vg, wobei unter vg die mittlere Durchflußgesch 


durch die Zylinderfläche verstanden wird, so kann diese Größe über das Schicksal des 
einen gewissen Aufschluß vermitteln. Bei a ee | 


Kr 2—k > en 
va d;' Ei 


2) W. Barth: ‚Zur Probier Er a 
Beziehung (11). Dr ra 
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wovon der eine von der Zentrifugalkraft, der andere von der Schwerkraft abhängen wird, so 
Br gelangt man zu folgendem allgemein gültigen Ansatz für den Abscheidegrad 2: 


j 2—k k ’ 
re 
2 - ag / % (6). 


2 ala) 


SEE Der Auswertung der Versuchsergebnisse wurde 
Beziehung (6) zugrunde gelegt. w,/v; und w;/vg wurden 
« errechnet, wobei %/v, — 1,96 und k = 0,8 angenommen 
wurde. Die Größe w; wurde mit Hilfe von Sinkgeschwin- 
digkeitsbestimmungen ermittelt. Die Ergebnisse sind 
in Bild 3 zusammengestellt. Es zeigt, daß die Ver- 
suchswerte für kleine wy,/v; in eine Kurve fallen, 
und daß sich die Ergebnisse durch den Ansatz (6) 
= darstellen lassen. Die Kurve Bild3 kann auch als 
2 Fraktionsentstaubungsgradkurve des betreffenden Ab- 
scheiders für V=const gedeutet werden. Mit Hilfe Be 
dieser Kurven lassen sich die voraussichtlichen Abscheideleistungen für Großausführungen 
errechnen, wobei bei Gasen im allgemeinen k= 1 gesetzt werden kann. 


Abscheidegrad— ——- 


4. Wr 
Verhaltnis TEE 


Se im Übergangsraum der schnellaufenden Wasserturbinen 
Von M. Strscheletzky in Kreßbronn 


Die Strömung im Übergangsraum zwischen dem Leitapparat und dem Laufrad rd unter 

Berücksichtigung der Form des Turbinenhohlraums und der Bauart bzw. des Öffnungswinkels 
der Leitschaufeln untersucht. Die exakte rechnerische Ermittlung der Leitschaufelwirkung 

mittels konformer Abbildung (Kreisgitter — gerades Gitter) ist nur für zweidimensionale Strö- 
_ mung möglich. Für den Fall einer dreidimensionalen Strömung im Turbinenhohlraum (bzw. im 
Leitschaufelraum) wird eine angenäherte Lösung betrachtet, die eine für die Praxis völlig aus- 
reichende Genauigkeit aufweist. Damit wird der Zusammenhang der Zirkulation /”(b) um die 
; Leitschaufel mit der Zirkulation J'(b) um die Schaufel .des entsprechenden geraden Gitters 


A ER) 


_ annähernd bestimmt, wobei k(b) einen Beiwert und b- die Leitschaufelbreite bedeuten. 


a Die vom Leitapparat gesteuerte Strömungist im allgemeinen Fall nicht wirbelfrei. Für eine 
E:: IE einige wirbelbehaftete Strömung wird die Differentialgleichung 


nn +) + at )=0 REARSIL IE. TO) 


er: angeführt, die längs jeder zu den Stromflächen senkrechten Linie gilt!). 


Br: Die Integration der Gleichung 2 längs der betreffenden Normallinie ergibt die meridionale 
Komponente der lokalen Geschwindigkeit 


Sc SE Val if er 
; Re ‚ a—2]1 flb)e ° SUB ende (3), 


00 dr 
ob +8. db 


er 10) = eu 


_ bedeutet. Für die Auswertung der Gleichung 3 wird ein numerisches Verfahren angeführt. 


& ER Die für verschiedene Leitapparatöffnungen eines Turbinenhohlraums berechneten Strö- 

> mungen sind durch Versuche nachgeprüft worden, wobei die errechneten und gemessenen Werte 

ER ür die drei Komponenten der lokalen Geschwindigkeiten sehr gut zusammenfallen. 
.— 

Fr y Cu und cm bedeuten 2 Umfangs- und Meridionalkomponente der lokalen Geschwindigkeit, r — ist 

& en Punktes von der Tuzkippnsekes; und R ist der lokale Krümmungsradius der 


e- 
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Lösung einer Aufgabe über Strömung in Leistüngssystemen 
Von W. Richter in Frankenthal 


Für stationäre, turbulente Strömungen in einer überall dichten, einfachen Leitung, in 
welcher dem strömenden Mittel keine Energie zugeführt wird, besteht zwischen dem Förder- 
strom Q und der Gesamtdruckdifferenz P;, zwischen dem Zuströmpunkt z und dem Abström- 
punkt % die Beziehung 

Be 
Die Konstante R soll R- Wert genannt werden. 

Ein Leitungssystem haben folgende Form: im Zuströmpunkt des Systems teilt sich die 
Leitung in zwei Teilleitungen, welche sich im Abströmpunkt des Systems wieder vereinigen. 
Von einem Punkt der einen dieser Teilleitungen führt eine Querverbindung zu einem Punkt 
der zweiten dieser Teilleitungen. Die R-Werte der fünf das so beschriebene System bildenden 
Leitungszweige und die Gesamtdruckdifferenz P des Systems seien gegeben. Die Aufgabe, den 
in jedem Leitungszweig fließenden Förderstrom zu bestimmen, führt auf die Gleichungen 

ıRı @ +1RQ9, = 2,RQ,9; + P 
R9;, +:R = — 2:89,09 + R 
Rd — R,® +BR9%=0, 
für deren graphische Lösung ein Fluchtlinienmonögramm entwiekelt wurde, welches aus zwei 
Hyperbeln und aus Geraden als Leiternträgern besteht. An einem Beispiel wurde die Auslegung 
eines für ein solches gegebenes Leitungssystems geeigneten Lüfters besprochen. 

Wenn in einem der Leitungszweige des Systems ein Lüfter arbeitet, so ist der betreffende 
R-Wert nicht mehr konstant. Er ist als Funktion des entsprechenden Förderstromes durch die 
Lüfterkennlinie festgelegt. Das Nomogramm enthält dann zwei Scharen von Hyperbeln als 
Leiternträger. Auch dafür wurde ein praktisches Beispiel durchgeführt. Schließlich wurden 


Sonderfälle besprochen, die sich daraus ergeben, daß einzelne der gegebenen R-Werte null oder 
unendlich groß sein können. 


Das Problem ist anwendbar auf Pumpenleitungen, Tunnelbelüftung, Grubenbewetterung, 


Saugzug- und Unterwindanlagen und andere industrielle Ventilationszwecke, wofür seine LöBung 
Hinweise auf die Wahl der geeigneten Pumpen bzw. Lüfter gibt. 


Der Wärmeaustausch zwischen einem geheizten Band 
und einer Konvektionsströmung 


Von M. Herbeck, 
Max-Planck-Institut für Strömungsforschung, Göttingen 


Das Problem des Wärmeaustausches zwischen einem geheizten Band bzw. einer geheiztän 
Halbebene und einer Konvektionsströmung parallel zur Wandoberfläche ist von meßtechnischer 
Bedeutung (vgl. Reichardt[1], Schuh[2], Ludwieg[3]). Das Band bzw. die Halbebene 
sind in der z-Richtung unendlich ausgedehnt angenommen; die Konvektionsströmung erfolgt 
in der «-Richtung. Weiter wurden für die Lösung dieses ebenen Problems bisher folgende ver- 
einfachende Voraussetzungen gemacht: 

a) Die Strömung ist inkompressibel und ihre Geschwindigkeit « wächst proportional dem 


Wandabstand y nach dem Gesetz u—= y- — ‚ wie es einer idealisierten Reibungsschicht mit 


der Wandschubspannung 7, der Dichte o und der kinematischen Zähigkeit » des strömenden 


Mediums entspricht. 


b) Die auftretenden Temperaturdifferenzen sind klein, so daß die Stoffwerte als konstant 
angesehen werden können. 


c) Der Wärmeaustausch durch Leitung in STÜNSUDESDERELDE wird gegenüber der ‚Wärme- 


konvektion vernachlässigt. 
Bei der vorliegenden Untersuchung wird die letzte Annahme fallen gelassen, da sie im 
Gebiet kleiner Strömungsgeschwindigkeiten und für kleine Bandbreiten sicher nicht zulässig ist. 


Aus der allgemeinen Differentialgleichung für das Temperaturfeld T beim Zusammenwirken von - 


Konvektion und Leitung: 
“6 (w grad T) — div (Agrad T)= 0 
e spez. Wärme BE > 
7 Wärmeleitfähigkeit 


Ex 
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wird dann nach Einführung der neuen Koordinaten . 


E—= pl? n=ßRy mit = i (1) 


DIS? 
SUke 


die Differentialgleichung 
Ter u 25 Er nTe ee en, NE RE (2a) 


mit dimensionslosen Größen erhalten, wozu für die geheizte Halbebene die folgenden Rand- 
bedingungen hinzukommen: 


De] E>0 „= 0 
T-0 E<0 ME RN (2b).. 
El) E< 0 N—> © 


Ist T(«, ı „ die Lösung für die geheizte Halbebene, so lautet für das Heizband der Breite 5 die 
Lösung T= T(a, y)— T(2—b, y). 


Gleichung (2a) ist eine elliptische Ditferentialgleichung mit Singularitäten bi&=n=0 
und im Unendlichen, die durch die Randbedingungen erzwungen sind. Die schwere Singularität 
im Unendlichen würde bei der Behandlung des endlich breiten Bandes nicht auftreten; 
diese mathematische Komplikation wurde jedoch in Kauf genommen, um mit einer einzigen 
Rechnung beliebige Bandbreiten erfassen zu können, wie oben gezeigt wurde. 
Zur Bestimmung des Temperaturfeldes in der Nachbarschaft von &= n= 0 werden Polar- 
koordinaten r = (&?+ n?)!? und 9—=arc tg n/£ eingeführt und man erhält 
1 | 
In+7 7,+ 22 = sin PECOSUIE LE ECT 5 Dar else (3) 
Das gesuchte Temperaturfeld 7 wird als ein durch Wärmekonvektion gestörtes Feld eines ur- 
sprünglich reinen Wärmeleitungsfeldes betrachtet. Mit dem Lösungsansatz 
© =} Di) 
T=3rf,d)+hnr 3 r(9) und a1 Bar 
v=0—0 u=4 
einer Lösung für reine Wärmeleitung, die die Randbedingungen für „=0 erfüllt, und 
(0) = q(9) = 0 für » #0 und 9—=0 und d9=n, wird 
0, 
rss —. ra, sind + 12asin 29 _. > — — sin >) 
j 8 87 
sind  5sin = 
3. 4 4.9 —__ x 
3 dag: sin30-+re[a, sin a0+% sin? d 187 + 576% 
5rtlnr | 
a er SURET 5 in59®--2 int ® le 
768% sin 494-7? |assin 5 9 + 2/cos 0 sin 0) (4) 
en ; ‚__cos?dsin’d dsin69 ä 3 43 ) 
+17 las sin 6.9 ein Bm Eu do sin? d 1608 x 4a 
a; er Ze 
+ sin® a _ “+ Ta | a —5sintdcos®d + .... 
Die hier noch unbestimmten Koeffizienten a, müssen zur Erfüllung der Randbedingungen bei 
rsnd = n—% dienen, was in der Lösungsform (4) für T praktisch nicht möglich ist. Es wird 
deshalb für die Differentialgleichung (2a) noch eine 2. Lösung nach der Methode der Separation 
der Variablen &, n aufgestellt, die nach Berücksichtigung der Sagen und Benutzung 
- des Fourier-Integrals lautet: 
er miete 
T=— + an SR: (DR 
| 2 .2n Jvolal (1) 
| | ; Fri Hı)3 =; | a 
> wo HA), die 1. Hankel-Funktion der Ordnung 1/3 ist. 


E Anschließend wird für die beiden Lösungen (4) und (5) der Temperaturgradient auf der 
ERRTE Wandoberfläche 


Be 1 j $ m 
Pr a ee) nenn (6) 


Y# 


m urn 
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bestimmt und zur Übereinstimmung gebracht. Auf diese Weise erhält man die Koeftizienten Er . 
4 0,3813 Son, — | 0,0318 ee 
4; 2 —0,0098 Au 0,0014 a ER a Re - 
a,= 0,0005 : de 9,0001 USWFRT > RN: 


Damit ist die ‚innere Lösung“ (4) bestimmt. 3 Pr ae = 
Aus (5) gewinnt man ferner die asymptotische Lösung für große positive & und n = 
T =1—n (0,538 E12 22.0,132£°1°) s 
un T. (0,0398E-43 SE 0,098£-1113) 


. —n 5(0, 0155873 +0, 0848179). en = 
Zum Schluß sei die Nusseltsche Zahl Nu für den Konvektionsanteil des Wärmeübergangs ger © 
geben, d.h. 


Nu Hi ae 


ist die in der Zeiteinheit bei der Temperaturdifferenz 1 von der ‚Längeneinheit eines Heizbandes 
mit der Breite b abgegebene Wärmemenge vermindert um den reinen Wärmeleitungsanteil Ben 
dividiert durch die De 


Für kleines b V=- — b ist En / e = SE ne 3 N 3 


ur 
gs 


We 0625 32 4 0,0066 db — 0,0041 b* In d+ 0,0001 Dr 
Für 5>1 ergibt sich die Näherungsformel : FE - Sets 
Nu = — 0,807 b2B + 0,198 522 EB eh b -+ 0,510 + 0,043 BE  @b)- 


Man erhält also für sehr große b asymptotisch die 'bekannte Näherungslösung [1], [a]: Nuz x Be ER 
—0,807 b2B3, d.h. Nw ist proportional der 3. Wurze] aus der N ER Dagegen. Se 


ist für kleine b-Werte Nu proportional r. ae 
Die bisher von mir gemessenen Werte von Nu sind etwas BETT, x ans herein Werte. SET 
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Über ein Minimumproblem aus der Tragflü ügeltheorie 
Von K. Nickel in Tübingen JE 


- I 


s- Es wurden Ergebnisse aus einer Arbeit angegeben, die demnächst unter den Titel „Lösung. 
Eu eines Minimumproblems der Tragflügeltheorie‘ in det r- Zeitschrift erscheinen. ‚soll. 
S 


Der Einfluß von Staubgröße und Staubdichte Auf die Druc 
bei beschleunigten Strömungen 


Von Hans Calame in Karlsruhe 5 


Re. Beider Beschleunigung eines staubführenden Luftstromesi in der Ve gun 
A Stromröhre stellt sich ein höherer Druckabfall Ay, ein als bei eg ft; 
: digkeit (Apı). Eu. ; 

Die mitgeführten Staubteilchen werden von der Trägerluft.: 

infolge ihrer großen Massenträgheit (spez. Staubgewicht yx : 
geschwindigkeit angeströmt und durch die Reaktion ihrer 
mungsweges s (s— Länge der Verengung) auf eine eenias 

schleunigt. 2 
Aus dem RA SEN (geschwindigkeitsabhängigen 


ee 5 4. ys Pr 
EISEN ar} 2 
Au “ ; : ERR: e 
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A wobei der Staubgehalt RE (8 est) 9 
Ys 7 7m» Luft /sec A 


kann die Austrittsgeschwindigkeit ce, errechnet werden; die zur Beschleunigung der Staubkörner 
aufgewandte Druckenergie Ap, = Ap, — Apı ergibt sich aus der Messung. 

Handelt es sich um einen Staub kleiner Korngröße !), so zeigen Rechnung und Versuch, 
daß die Staubteilchen auf eine sehr hohe Geschwindigkeit, gröbere Teilchen dagegen bei gleicher 
Düsenlänge und gleicher Luftgeschwindigkeit nur auf eine mäßige Austrittsgeschwindigkeit 
beschleunigt werden.. 

Für eine Düse bestimmter geometrischer Abmessungen läßt sich der Geschwindigkeits- 
verlauf von Staubteilchen verschiedener Korngröße a, in Abhängigkeit von der Düsenlänge s 
und der Luftgeschwindigkeit v, bestimmen und als Funktion 


a AK 9 a De En:) 


mit dem Parameter «a, graphisch darstellen. h 

' Der aus (1) für s= s,errechenbaren Staubgeschwindigkeit O, ist somit nach (3) eindeutig 
die Staubgröße « zugeordnet. Damit ist ein Zusammenhang zwischen der meßbaren Größe 
Ay. = Ap,; — Apı und der Korngröße a gefunden, der exakte Rückschlüsse auf die Größe des 
mitgeführten Staubgutes, bzw. — bei bekannter Korngröße — auf den Staubgehalt des Luft- 
stromes gestattet. 

Umgekehrt kann auch der durch Staubbeimengung zu erwartende höhere Druckabfall 
eines Luftstromes in einer Verengung (Düse) bei bekannten Kenngrößen des Staubes (y,, yx, «) 
berechnet werden. 


Theoretische Grundlagen: Vernachlässigt man die Schubspannungen in der 
Strömung und ordnet jedem Düsenquerschnitt eine gleichmäßige mittlere Luftgeschwindigkeit 
nach der Kontinuitätsbedingung zu, — das spez. Staubvolumen ist verschwindend klein und 
kann vernachlässigt werden — so ergeben sich bei der zylindrischen Meßdüse und der 
Düse mit linearem Luftgeschwindigkeitsanstieg in Strömungs- 
richtung (s) verhältnismäßig einfache Beziehungen: 


Zylindrische Meßdüse: Die bei der zylindrischen Meßdüse auftretenden Ver- 
Fr hältnisse wurden bereits früher von W. Barth untersucht und im Ing.-Archiv 1949, Heft 3/4, 
S. 147/150 mitgeteilt. 


Düse mit linearem Luftgeschwindigkeitsanstieg längs s: Be- 
schränkt man sich bei dieser Düse auf die Untersuchung feinster Staube, d.h. auf Anström- 
Be, verhältnisse der Staubteilchen innerhalb des Stokesschen Bereiches Re <1, so zeigt sich, 
En daß auch die Geschwindigkeit der einzelnen Staubkörner linear mit s ansteigt. Der meßbare 
zusätzliche Druckabfall Ap,, läßt sich in dimensionsloser Schreibweise geschlossen ausdrücken 
durch: z 5 

EN LETTER ER RE (+) SEHR 3 
2 %q Vz So g 


Er = / Ur "Va 
ag Pe 0,5+] 0 , 


Versuchsergebnisse: Im Laboratorium für Strömungslehre der T. H. Karlsruhe 
- wurden 1949 im Auftrag und mit Mitteln des VDI Versuche durchgeführt, um die abgeleiteten 
Beziehungen zu überprüfen. 

Aus dem zusätzlichen Druckabfall Ap,, wird die Korngröße a rechnerisch bestimmt und 
3 mit der durch Mikroskop ermittelten Staubfeinheit verglichen. Die Ergebnisse sind unten 
zusammengefaßt. 

2 - — Für die Durchführung der Versuche wurde eine verhältnismäßig einfache Versuchsapparatur 
benötigt: Ein Ventilator saugt in einem Rohr Luft aus der Atmosphäre durch die Staubmeßdüse 
an, welche bei mitgeführtem Staub den Druckausschlag Ap,, bei staubfreier Luft den Druck- 

ausschlag 4»; liefert. 
Der Staub wird in einem Staubfilter hinter der Meßdüse vom Luftstrom getrennt. Das 

 sekundlich durch die Düse geförderte Staubgewicht @ kann durch Auswägung des Filters und 
Berücksichtigung der Versuchsdauer ermittelt werden. Die Luftmenge V (m?/sec) wird mit einer 
hinter dem Filter angeordneten genormten Meßblende gemessen. > 
Die experimentellen Untersuchungen erstreckten sich auf die oben beschriebenen Düsen- 
- formen; bei der Düse mit linearem Geschwindigkeitsgesetz wurde die Länge s, bei sonst gleichen 
2 Abmessungen in drei verschiedenen Ausführungen variiert. 


S 


4) Unter der Korngröße ist bei allen folgenden Betrachtungen der Durchmesser der als Kugeln auf- 
ten Staubkörner zu verstehen, die am Ende der Verengung die gleiche Geschwindigkeit ca wie die unregel- 
'iger geformten originalen Staubteilchen erreichen, diesen also widerstandsmäßig gleichwertig sind. 


» 


r > 


L 


I 


» aBSPE u 2) 
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Die Versuchswerte sind in untenstehender Zahlentafel zusammengefaßt. 


Auswertung der Versuchsergebnisse: Der zusätzliche Druckabfall Apyı 
wächst nahezu linear mit dem Staubgehalt y, an. 

Ein Vergleich der bei variierten Düsenlängen ermittelten Korngrößen bestätigt die für 
feine Staube gültige und in (4) enthaltene ale uepeöingung (w; * va)/(s  g) = konstant. 


— Die Korngröße feiner Staube von 
3 bis 5 u wird von der Meßdüse mit 
linearem Geschwindigkeitsgesetz noch 
verhältnismäßig gut erfaßt. Eine Un- 
terscheidung dieser Staube ist noch 
möglich. 

Die zylindrische Meßdüse mißt 
die Korngröße feiner- Staubteilchen 
nicht exakt; die Ursachen dieser Ab- 
weichungen können zum großen Teil 
erklärt werden. 

Die Zusammenhänge zwischen 
dem zusätzlichen Druckabfall Ap,, und. 
der Staubkorngröße ermöglichen die 
Anwendung eines neuartigen Verfah- 
rens zur Messung kleiner Korngrößen. 

Neben dem Vorzug eines verhält- 


. | Staubfeinheit, durch Versuche ermittelt 
z indrische Düse mit linearer 

3. Meßdüse | Geschwindigk- Charakterist. 
Sr EIEIEUE EEE 
en > 
lea» | 
Mai. 
Mais gr 
ee er 

nismäßig kleinen Apparateaufwandes 


Zusammenstellung der Meßergebnisse. 
Korngrössenermittlung. Karlsruhe 5 : ER: \ : ; 
et eröffnet es die Möglichkeit, Feinheits- 


Bild i bestimmungen von Stauben in einem 

bisher schwierigen Korngrößenbereich 

durchzuführen. Darüber hinaus kann die gleiche Meßanordnung für die Mengenmessung pneu- 

Bene geförderter Staube durch Bestimmung des Staubgehaltes y, der Strömung Verwendun 8 

inden 

In neuerer Zeit wurden gleichlaufende amerikanische Arbeiten vonCarlson,Frezier 

und Bengdahl (Transactions ASME 1948) bekannt, die sich der gleichen Meßmethode, 
jedoch zweier hintereinandergeschalteter Staubmeßdüsen bedienen. 


Ein einfaches Singularitätenverfahren zur Erzeugung von 
Schaufelgittern 


Von N. Scholz in Braunschweig 


Achsial durchströmte Leit- und Laufräder von Turbomaschinen können in bekannter 
Weise durch ebene unendlich lange Schaufelgitter idealisiert werden. Für die Beurteilung der 
Güte eines Schaufelgitters ist die Kenntnis der potentialtheoretischen Geschwindigkeitsverteilung 
auf der Schaufelkontur wegen ihres starken Einflusses auf das Verhalten der Grenzschicht von 


„ ausschlaggebender Bedeutung. Die bisherigen Verfahren zur Ermittlung der Druckverteilung 


am Gitterprofil beruhen zum größten Teil auf konformer Abbildung und sind in derrechnerischen 
Durchführung entweder recht langwierig oder aber zu ungenau. Wir haben uns daher die Auf- 
gabe gestellt, ein einfaches und bequemes Verfahren hierfür zu entwickeln, womit größere Reihen 
von Schaufelgittern potential- und grenzschichttheoretisch durchgerechnet werden können!). 

Hierbei schließen wir uns an die Vorstellung der Birnbaumschen Theorie der tragenden 
Fläche an, die bereits von E. Pistolesi?) auf den Fall der unendlichen Flügelreihe übertragen 
wurde. Wir ersetzen also das Gitter in der komplexen Ebene z=x-+-iy durch eine kontinuier- 
liche Wirbelbelegung y (x) auf der x-Achse sowie auf den Graden y=»vt—=0; +1; #25...) 
(t= Schaufelteilung) von = — b[/2 bis + 5/2 (b=Gitterbreite). Diese liefert eine komplexe 
induzierte Geschwindigkeit w(z)=w, —iw, von 


+ 
w(2) = 5. | y («’) ctgh E @—a)lar. 
«= —b/2 ; 
Im wesentlichen interessiert nun die Zusatzgeschwindigkeit am Profil selbst, also für x. 


!) Vgl. hierzu H. Schlichting und N. Scholz, „Über die theoretische Bergabnına der Strömung 
verluste eines ebenen Schaufelgitters“, (Erscheint demnächst im Ing.-Archiv 


v.) 
?) E. Pistolesi, „Sul calcado di schiere onfinite di ali sottili“. L’aerotecnica, Vol. XVII, fasc. 6, 1937. 


x 
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Diese besitzt die Komponenten +5/2 
x) 1 
wi(2; 0): en; ul )=— 5, | y (a’)etgh F@-2) da’. 


"= —b/2 


Für die Ermittlung der Schaufelkontur hat man die kinematische Strömungsbedingung 


(2) 
4%) Kontur 


wobei sich im allgemeinen Fall w,(x) aus der induzierten =, 
Abwärtsgeschwindigkeit w,; und einer überlagerten 
Translationskomponente wy,» zusammensetzt. %,& + 0 
bedeutet eine Staffelung des Gitters um den Winkel 
cetg Pk, = Ww, [Wg» (Ps Winkel zwischen Translationsge- 
schwindigkeit w» und x-Achse). Die Staffelung wird also 


durch Überlagerung einer 


parallel zur Gitterfront erzeugt. Dies stellt eine Näherung 
für nicht zu große Staffelungswinkel ß, dar, die in einer 
ausführlicheren Arbeit näher begründet werden soll. 

j Bei der Berechnung des Abwindintegrals ist die 
Singularität des Integranden bei @—=x’ noch etwas stö- 
rend. Die N ralnne der singulären Stelle ergibt 

2 


+b/ 


Wyi(z; =, [ KAUM, 


«= —b/2 
+b/2 


1 ; 
FRE / y(&) | ctgh 
«= —b/2 


Das erste Glied ist das bekannte Abwindintegral des 
Einzelflügels, während das zweite den Anteil aller übrigen 
Schaufeln darstellt. In dieser Form kann für vorgegebenes _ 7%, 
y(x) die Abwindverteilung und die Schaufelkontur leicht 


berechnet werden. . 
Als Beispiel haben wir 


y (z) in 
Erz ’ 
Wo HE) 08 


Translationsgeschwindigkeit 


! IR N 


aa) 


Bee]- 


ein Gitter mitelliptischer Zir- 


kulationsverteilung gewählt, (das bei verschiedenen Tei- mar. Schaufelgitter: verschiedener Teilung t/b 
lungsverhältnissen t/bgleiche Ablenkung bewirken soll(Ge- Mit gleicher Strömungsumlenkung bei elliptischer 


Zirkulationsverteilung 


samtzirkulation I’ der einzelnen Schaufel proportional der 4: zirkulationsverteilung. — b: Abwindanteil der 
Teilung t). Bild1 zeigt die Ergebnisse dieser Rechnung. Der zinzeischaufelw, = 9 w BE 

. . . . . . m y I Wyo + 6: Abwind- 
Gittereinfluß erreicht bei 2/b = 0,5 etwa die gleiche Größen- anteilaller übrigen Schauteln. _.d: Schaufelkon- 
ordnung wie der Einfluß der Einzelschaufel. Diesich daraus turfür Gleichdruckgitter bei 75° Umlenkung 


ergebenden Geschwindig- 


. keitsverteilungen auf der, *”r 


We 


Schaufelkontur, die sich mit 
Hilfe eines Hodographen be- 
quem ermitteln lassen, zeigt 
Bild 2. Bei einem gestaffel- 
ten Gitter wird, wie dort er- 
sichtlich, die Schaufelkontur 
nicht nur gedreht, sondern 
gleichzeitig längsverzerrt, so 


daß eine andere Profilform 


_ entsteht. Die numerische 
Durchführung solcher Rech- 
nungen gestaltet sich außer- 

-ordentlich einfach und 


. nimmt nicht mehr als zwei . 
5° Stunden für ein Beispiel in 
” Anspruch. 

2 a Abschließend sei er-. 
Br. wähnt, daß mit der gleichen 


Grundvorstellungdurch eine 
Quell-Senkenbelegung auch 
eine Dickenverteilung er- 
zeugt werden kann. 


4 


Da 
nu 
Kl 
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Bild 2; Geschwindigkeitsverteilung auf der Schaufelkontur für Schaufelgitter nach Bild 1 

a: Gleichdruckgitter w, = ws; ß, =102,5°; ß, = 78,5°. — b! Beschleunigungsgitter w,/w, = 0,66; 
ßı =60°; B2=35°. (Umkehrung der Strömungsrichtung liefert ein entsprechendes Verzö- 
gerungsgitter); (t/b»=0 ergibt die Stromiadentheorie mit ‚‚schaufelkongruenter‘‘ Strömung) 
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Die mathematischen Grundlagen zu einer Meßmethode 
des Schubmoduls zäher Flüssigkeiten 


Von E. Pier in Kassel | f Te 


1941/42 sind im Festigkeitslaboratorium der Technischen Hochschule Danzig von Herrn 
Prof. Dr.-Ing. H. Fromm Messungen zur Bestimmung des Schubmodüls zäher Flüssigkeiten 
veranlaßt' worden, die von Herrn Dr. Hartun 8 ausgeführt und ausgewertet wurden. Das 
Danziger Gerät bestand im wesentlichen aus zwei konzentrischen Zylindern, die durch einen 
schmalen Ringspalt von der Breite ö, der zur Aufnahme der Versuchsflüssigkeit diente, geuenit 
waren. Der äußere antreibende Zylinder wurde in Schwingungen nach 


yv=—a-.cosatl...... ; EEE 


versetzt, die durch die Zähigkeit und Elastizität der Versuchsflüssigkeit auf Er inneren frei- 
schwingenden übertragen wurden. Ausgewertet wurden die Versuche nach dem Maxwellschen 


Ansatz 


E25 


=. 1+g ee 


(7 Schiebungsgeschwindigkeit, n Viskosität, 2 Schubmodul und 7 Schubspannung). Weiter wurde 3 
angenommen, daß der Geschwindigkeitsabfall im Ringspalt konstant ist. Damit ergab sich ein 


einfaches Auswertverfahren. 
Nun hat später Prof. Fro mm!) den en Ansatz S durch die, ‚Tensorbeziehung 


x = Bot gie ee) 


(3, Deviator des Verzerrungsgeschwindigkeitstensor, B, Deyialar des 
Spannungstensors) wesentlich verfeinert und die sich aus diesem Gesetz 
ergebenden Formeln entwickelt. Da bei den Danziger Auswertungen Se 
Unstimmigkeiten auftraten, die damals nicht erklärt werden konnten, 9 
regte Herr Prof. Fro mm an, eine genaue mathematische Untersuchung rg 
der Danziger Versuche zu machen, die ich ausgeführt habe und deren: 
Ergebnis hier mitgeteilt wird. 
Bild ı Da der Ringspalt ö und die Amplitude a sehr klein gegenüber 
dem Zylinderdurchmesser bzw. Umfang waren, kann mit karthesischen 
Koordinaten gerechnet werden. Legt man das Koordinatensystem wie in dem Bilde, SO ef» 
hält man für den ‚Beharrungszustand das Randwertproblem: { 
oc do... | DIS, 


Stoffgl. EAAERE GE 
offgl. ERS 0 +0/9), N RR re e 


Y 


- 


ya rt 0 . 
(7-2: 0=0, = 0,5 ee Bu0z ke 2 : ® BR ER ® 


Pr Erg sdichte) BE RS 


Randbedingungen: (0, Vu a. d (6, 9 assinat a de BEE 


HE Mk 


(u reduzierte Massenflächenbelegung). 


Dieses Randwertproblem ist nicht linear. Daher wurde eine schrittweise Näherung 
Bis zu einer gewissen Größe von & — diese wurde bei den Danziger Versuchen nicht übers 
— kann o/@ neben 1 vernachlässigt werden. Dann braucht man die ‚zweite a. von 9 Br 
und die Stoffgleichungen gehen in den vereinfachten Ansatz (2) über. EREN 


Das Gleichungssystem für v und 7 ist dadurch linear. Ersetzt 
mittels (5) durch r, so erhält man ein lineares partielles Randwe: 
sich für den Beharrungszustand mittels des Ansatzes 


r=f(e) sin at + (2) cos RER 
lösen, Setzt man (Mi in das Randwertproblem ein, so erhält ana 
gleichung für , zwei gewöhnliche simultane lineare Differen 
konstanten Koeffizienten, die sich leicht lösen lassen. Es 
bis A, auf, die aus den Randbedingungen bestimmt: 
Gleichungen für die Konstanten. 
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‚Eine besondere Ausführung des Danziger Geräts zeichnete durch Überlagerung der Schwin- 
gungen des angetriebenen und freien Zylinders mittels eines gespiegelten Lichtstrahls auf einen 
P Schirm Lissajous-Figuren. Diese arteten in eine Gerade aus, wenn die Phasendifferenz der 
Schwingungen der beiden Zylinder 90° betrug. Die dazugehörige Frequenz soll kritische Fre- 
- - quenz genannt werden. Diese wurde in Danzig gemessen. 

Aus der kritischen Frequenz %,, kann man bei kleinen Frequenzen, die in Danzig nicht 
überschritten wurden, den Schubmodul @ der Versuchsflüssigkeit nach 


E%, N fr 1 RENNON: 5 
Fi 5 k ud een Vogt = . . . . . . . . . (8) 
i berechnen. Die Hartungsche Lösung erhält man,' wenn man in (8) o gleich Null setzt, denn die 
- Annahme eines konstanten Geschwindigkeitsgefälles ist damit identisch. Hartung wertete 

die Versuche also nach 


Ahr 
ne) 
aus. i 
(9) ergibt für Newtonsche Flüssigkeiten (#— ©) als kritische Frequenz ay—>o. Die 
Beobachtung einer kritischen Frequenz wurde also in Danzig als Beweis einer Maxwellschen 
| Flüssigkeit angesehen. Bei Beachtung der Flüssigkeitsmasse (0 #0) ergibt sich dagegen auch 
® für Newtonsche Flüssigkeiten eine beobachtbare kritische Frequenz und die Abweichung zur 
Maxwellschen Flüssigkeit ergibt nur eine Verschiebung dieser Frequenz. Das Kriterium der 
Maxwellschen Flüssigkeit hängt also nicht an dem Auftreten einer kritischen Frequenz, sondern 
in der genauen Erfassung eines gewöhnlich kleinen Frequenzunterschiedes. 
 Da’man in Danzig nach (9) auswertete, achtete man nicht auf eine peinlich genaue Fehler- 
diskussion. Deshalb ist es sehr schwer, diese jetzt nachträglich zu machen. Es ergibt sich aber, 
daß bei Newtonschen Flüssigkeiten V, statt eins auch zwischen den Grenzen 0,75 < Y, s 1,50 
liegen kann. Dieser Wert ist auf der sicheren Seite für unsere Betrachtung, d.h. die Schranken 
sind eher enger zu setzen. 

Untersucht wurden in Danzig Phtalate, Glyzerine, Öle und Teer bei verschiedenen Tem- 
peraturen (20 bis 60°C). Bei den meisten Versuchen liegt V ‚innerhalb der angegebenen Schran- 
ken. Hier läßt sich nicht angeben, ob die Flüssigkeit Newtonschen, Maxwellschen oder gar eines 
weiteren Typus ist. Fünf Werte zeigen ein kleineres V, als 0,75. Hier wird nach (8) @ positiv. 
Es kann also ein Maxwellscher Typus vorliegen. Drei Werte liegen deutlich über 1,50. Hier 
würde @ negativ werden. Da dies nicht möglich ist, muß ein weiterer unbekannter Flüssigkeits- 
typ hier bestimmt vorliegen. Die Untersuchungen zeigen also, daß es Flüssigkeiten gibt, bei 
denen das Newtonsche Gesetz 'nur eine Annäherung ist. 


; Über die Lösungen nichtlinearer partieller Differentialgleichungen 
-. - vom Reibungsschichttypus 


Von H. Görtler in Freiburg i. Br. 


00. Für die Berechnung laminarer Grenzschichten bei allgemeinem äußeren Druckverlauf 
ERX- besitzt man zwar eine große Zahl von analytischen und numerischen Verfahren — Reihenent- 
e' wicklungen, Differenzen- und Iterationsveriahren —, für keines dieser Verfahren jedoch eine 
N strenge mathematische Legitimation oder auch nur allgemeine strenge Abschätzungen für die 
Be Güte der jeweils erzielten Approximation. Völlig in der Luft hängen in dieser Hinsicht die 
2 schneller arbeitenden ein- und mehrparametrigen Näherungsverfahren nach dem Muster des 
- Pohlhausenverfahrens, die auch gelegentlich gänzlich versagen. Man ist daher auf die Erprobung 
= der Verfahren durch Vergleich mit einigen in Sonderfällen bekannten exakten Lösungen, ferner 
Sur auf die gegenseitige Kontrolle bei Anwendung verschiedener Verfahren 5 
und schließlich auf Rechenproben, wie die Erfüllung der Grenzschicht- 
bindungen oder gewisser Integralbeziehungen angewiesen. PIE 
Es liegt bei diesen Strömungsproblemen eine Rand- und An- | & 2 16, 
l 


gswertaufgabe vor, bei der — vgl. Bild, in der die positive 
Richtung die ausgezeichnete Sn erlun8 der ie 
i darstellt — in einem fangsquerschnitt z=x, das 

NV eitsprofil [7 Elemente a Geschwindigkeit in lleenden: Geinam 
De Yurgegeben ist, während auf zwei von dessen End- rg Nr 

ehenden Linien stromabwärts die Randverteilung von u vorgeschrieben ist. Der 
ierenden Strömungsgebietes B, auf dem damit insgesamt u gegeben ist, 
. (Es ist Berner De een ® der Geschwindigkeit auf einem Teilstück 
BA RE De 19 


- 
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von (, in dem Bild etwa auf dem unteren waagrechten Randstück, insoweit vorgegeben, daß sie 
aus u vermöge der Kontinuitätsgleichung %, + ®v,=0 durch Quadratur eindeutig ermittelt 
werden kann.) Einzelne Randteile von Ü können ins Unendliche hinausgerückt sein; so ist im 
Falle wandnaher Grenzschichten-der untere horizontale Rand als (abgewickelte) Körperkontur 
aufzufassen, während der obere Randteil nach y—= zu rücken ist, wo u asymptotisch in die 
gegebene äußere Geschwindigkeit U(x) übergehen muß; bei Nachlauf, Strahlausbreitung und 
Strahlvermischung sind die horizontalen Teile von Ü nach y=-- und — oo zu verlegen. 
Ein Schlußquerschnitt (©, (in Figur gestrichelt) in «= x, ohne Randvorgaben für die gesuchte 
Lösung schließt mit Ü das zu betrachtende Strömungsgebiet B ab; er kann unter Umständen 
stromabwärts nach = & verlegt sein. Das im folgenden betrachtete Gebiet B mit Rand C+C, 
möge in diesem dehnbaren Sinne verstanden werden. Die darzulegenden Zusammenhänge 
würden es erlauben, einen noch ganz erheblich allgemeineren Randverlauf zuzulassen, doch hat 
dies hier kein Interesse. Die Situation bei den wandnahen turbulenten Reibungsschichten und 
den Problemen der freien Turbulenz ist die gleiche, 

Im Bestreben, die Bewegungsgleichungen der laminaren und — nach bisher gemachten 
Ansätzen — der turbulenten Reibungsschichtströmungen gleichzeitig zu erfassen, legen wir die 


al 2a at mr fr a TE ee 


umfassendere, im folgenden als Differentialgleichung ‚‚vom Reibungsschichttypus‘‘ bezeichnete E 
Gleichung 2 

WU HU E NEU U U I ie Fan 5 
nebst Kontinuitätsgleichung u, + t, =0 zugıunde 1). Es wird dabei vorausgesetzt: 1. u(x, y) 1 


existierein B+ Ü + C, und sei dort stetig; 2. u(&, y)seiin B + C, zweimal stetig differenzierbar; 
3. Es sei u>0in B-+-C, (Ausschluß von Rückströmungsgebieten, für die im allgemeinen Ein- 
deutigkeit der Lösung schon aus physikalischen Gründen nicht erwartet werden kann); 4. die < 
Funktion f sei für alle vorkommenden Werte ihrer Variablen stetig und ferner bezüglich ihrer 
vierten Variablen monoton wachsend: f(&, % %, X4) — (1, 4 %, u) 0 falls X, — u, >0; 
anderen Bedingungen wird f nicht unterworfen. 

Der praktischen Berechnung von Näherungslösungen würde hier, wie auch sonst, am besten. 
gedient durch eine Anweisung der Art: ‚Wie auch die Näherung u, (x, y) des Geschwindigkeits- 
profils u(x,y) gewonnen sein mag, man prüfe sie einfach durch Einsetzen in die D.-Gl. (1) und 
schließe aus der Abweichung von der exakten Erfüllung nach einer gegebenen Vorschrift auf ö 
eine Abschätzung |u— w,|<'M des Fehlers der Näherungslösung selbst‘. Von praktischem 
Interesse wären verwandte Aussagen der Art: ,,Wenn man die Randvorgaben auf C, den äußeren 
Druckgradienten oder vorkommende Stoffwerte variiert, ändert sich die Lösung «(z, %), je nach 1 
jenen Änderungen, maximal um Soundsoviel in B+C+0,“. Dies und anderes leistet die : 
allgemeine Aussage, die ich im folgenden kurz der Struktur nach skizziere. Ihr Beweis kann, 
sobald man die geeignete Ausgangsungleichung hergestellt hat, in enger Anlehnung an den 
kürzlich von H. Westphal?) mitgeteilten überaus einfachen Beweis eines ein wenig spe- 
zielleren Satzes erbracht werden. ü 

Man betrachte unter den oben formulierten Voraussetzungen je eine Lösung der beiden 
Differential-Ungleichungen ae . 


ou; 9%; e On a 
MT day U a? Dye S& (= 1,2) 
u; ov; 


nebst Kontinuitätsgleichung ——=0 in B+(°). Nach den Randvorgaben sei 


02..:0% | F 
[4 — u] < A auf ©. Transformiert man auf die neuen unabhängigen Veränderlichen &= x und 
0 0 oe: | 
N= %g, Wo y,(x, y) die vermöge Se =, — Sn — v, erklärte Stromfunktion ist, so folgert 


man leicht, daß z(&, n) = u — u, einer Differential-Ungleichung der Gestalt 


ö 2? 
FF (me + a mie 2, An» 2m) 


in dem Bildgebiet B-+C, von B+C, in der &, n-Ebene genügt, wo F bezüglich seiner vier En. 
Variablen stetig ist und mit 2, monoton wächst. In F gehen ferner die Parameter e, und , ein. RL 


') Die Erweiterung der Betrachtung auf kompressible Medien wird, soweit ich sehe, keine wesentichen 
neuen Schwierigkeiten bringen. Die Übertragung auf die verwandten Probleme der Wärmeleitung und Diffu- 
sion ish fast unmittelbar möglich. RN SR RR a En ra Dean 
®) Math. Zschr. 51 (1949), 8. 690, SE RE RD 
. ._) Daß die Kontinuitätsgleichung nicht als genähert sondern als streng erfü 
keine Einschränkung. Man hat sich nur jeweils zu dem Geschwindigkeitsprofil u; (z, 
durch Quadratur aus der Kontinuitätsgleichung unter Beachtung der zur eindeutig 
Randvorgabe für vermittelt und dann diesesvi(z, y) in die Diff.-Gl. (1) zur Erm 
& = 0 eingesetzt zu denken. ER RE 
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Es gilt dann derfolgende Satz: Läßt sich einein B70710 stetIg®,siH BERG, 
zweimalstetig differenzierbare Funktion 2(&,n) angeben, für die 


a) lm +? Dare %, 29» mn) in B+ C, 
und \ ee 
b) 22, au 6 
sesorist-äuch A 
z2<z überallin B-+C. 
Bei im Unendlichen verlaufenden Teilen von Ü muß die Voraussetzung b) eine etwas modifizierte 
aber praktisch nicht engere Fassung erhalten. 

Es lassen sich ganz allgemein Majorantenfunktionen zZ angeben, jedoch wird man in jedem 
Einzelfall hoffen dürfen, bessere spezielle Majoranten zu finden, um u, — u, und analog u, — u, 
zu majorisieren und somit |u, — u,| in B ni C, abzuschätzen. 

Der Satz gestattet sofort den Einzigkeitsbe weis der exakten Lösung aller mit (1) 
erfaßten Strömungsprobleme bei den üblichen Randvorgaben auf Ct). Bei Annahme der Existenz 
zweier exakter Lösungen u, und % (4 =&,=0 oben) mit gleichen Randwerten auf C, also 
z= U — U —=0 auf (, ergibt sich 2 = @exp (—b£) als eine Majorante für jedes «> 0 bei geeignet 
wählbaremb. Indem man so z und —z majorisiert, folgert man |a|<a exp(-bE) in B+(,,und 
da b auch für a&—0 existiert, ergibt sich „— w,=0in B+Ü-+G(.. 

Mit diesem Einzigkeitssatz erhalten die oben formulierten Vergleichs- und Abschätzungs- 
sätze erst ihren Sinn, etwa: Vergleich einer Näherungslösung von (1) mit der exakten Lösung 
bei den gleichen Randbedingungen, usw. 

Eine ausführliche und präzisere Darstellung wird in Kürze veröffentlicht werden. Ihr soll 
folgen die Wiedergabe der inzwischen in Angriff genommenen Nutzanwendungen des obigen 
Satzes. 


Über ein gemischtes räumliches Randwertproblem der Potentialtheorie?) 
Von Hans Schubert in Rostock 


In einem Windkanal von kreisförmigem Querschnitt, dessen Achse wir zur &-Achse eines 
Systems von Zylinderkoordinaten &, r, # wählen, möge eine reibungslose inkompressible Flüssig- 
keit ein kleines Hindernis umströmen, das im Ursprung.angebracht sei; weit vor dem Hindernis 
habe die Flüssigkeit die konstante Anströmgeschwindigkeit v, = — V. Ferner soll der Kanal 
eine offene Meßstrecke aufweisen, die von den auf der Zylinderachse senkrechten Ebenen 2 = —a 
und & = abegrenzt ist; an der freien Strahlgrenze herrscht dann der konstante Druck des Außen- 
raumes. Zur weiteren Vereinfachung des Problems nehmen wir an, daß die vom Hindernis und 
den Kanalgrenzen hervorgerufene Störungsströmung ‚‚klein‘‘ gegenüber der Ausströmgeschwin- 
digkeit V ist und daher die Verformung der freien Strahlgrenze vernachlässigt werden darf. 

Als einfachste Form eines Störkörpers dient ein Druckdipol mit dem Beschleunigungs- 


potential 
- __Mcosd r 
y(%r,d) = SEE Sr, (158); 
’ dessen Normalableitung auf dem Kanalrand (r = 1) die Form E \ 
Be; ZR = Mm 'ı | 
ER es t ‚6 1,2 
B: [> = 6G,(8) cos 9 mi o(2) = ee) (1,2) 


annimmt. Die gesamte Störungsströmung denken wir uns nunmehr in die bekannte Strömung 
des Druckdipols und in eine noch unbekannte Zusatzströmung, mit dem Beschleunigungspoten- 
Be: tialx(% r, ®) zerlegt, ‘das so bestimmt werden muß, daß es allen Rand- und Grenzbedingun- 
Ic ‚gen genügt und im Innern des Zylinders eine reguläre Potentialfunktion darstellt. Für die 
„Fösung dieses gemischten Randwertproblems der Potentialtheorie machen wir den Ansatz 


% y x(2, 1,0) = L(x, r) cos 9 Bere eier... (21) 
4 & und schreiben die Randwerte der Normalableitung in der Form 
a gi ee a Fl A 99 
VER ERE ve = (4a) — — 6,(2)) cos d ae. (2,2). 


var Teer 
EN el) Meines Wieden ist für den ie, Problemkreis ein Einzigkeitssatz bisher nur für laminare 
andgrenzschichten. ausgesprochen worden (N. Piskounov, ©.R. Acad. Sci. URSS.27 (1940), S. 104 mit 
1 7# pen und mir nicht verständlichen Beweisskizze). 

B.ne Ausführliche; Arbeit des Verfassers über dieses Thema wird in einem späteren Heft der Zeit- 


Re 19* 
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Ts 


wobei die Funktion g(x) für |@| > a verschwinden muß. Zur Bestimmung von g(z) im Intervall 
— a<xsa dient die Integralgleichung erster Art ER AR : 


- | sol — 412-114 1,9} dı=9@) es, 


[77 v 
in der L(z,t) einen Zusatzkern mit stetigen partiellen Ableitungen erster Ordnung darstellt!) 
Hat man aus Gleichung (3) durch Übergang zu einem linearen Gleichungssystem mit überwiegen- _ 
der Hauptdiagonale die Funktion g(x) ermittelt, so erhält man das Beschleunigungspotential 
x(&, r, ®) der Zusatzströmung als Lösung eines zweiten Randwertproblems der Potentjialtheorie, . 
dessen Randwerte auf dem Zylindermantel durch (2,2) bestimmt sind. = Er 

Für den Abwind der Zusatzströmung in der Kanalachse findet man den Ausdruck 


sin q& 


wobei w(% (x) den Abwind in der Achse des durchweg geschlossenen Kanals bedeutet. In den 
Sonderfällen des Freistrahls (a = 00) und des durchweg geschlossenen Kanals (a= 0) ergeben 
sich für x («, r, 9) und w,(x) geschlossene Integralausdrücke, die mit den bekannten Ergebnissen 
von I. Flügge-Lotz im Einklang stehen. : = 


T 


Zum Schalldurchgang in zwei- und dreidimensionalen 
Düsenströmungen 


Von Hermann Behrbohm in Waldkirch/Brsg. 


Der erste Teil des Vortrags gab eine Skizze der Näherungsmethode zur Behandlung ebener 
isentroper Düsenströmungen in Schallnähe, wie sie vom Verfasser kürzlich entwickelt wurde?). 
„ Der zweite Teil des Vortrages befaßte sich mit einem Näherungsbeispiel zur isentropen drei- 
dimensionalen Düsenströmung in Schallnähe, unter der Annahme der Drehungsfreiheit. u, v, w 


Azk+17,; 4:-023; Yy»1 j a WER 

seien die x-, y-, z-Komponenten des Geschwindigkeitsvektors w,odi 
RU © w ; ERSTER N 

den Annahmen an 1: = <1, — <1 (c*= kritische Schal 
nuitätsgleichung approximativ SEELE 

ou 09V dw 
a0 ey 02 x 
1) Vgl. hierzu M. Hasse: Über eine singuläre Integral lei 

keit. Erscheint demnächst in dieser Zeitschrift, 1; 
w ®)H. Behrbohm, Näherungstheorie des unsymme 
2. angew. Math. Mech. Bd. 30 (1950), 8. 101— 112. 2 
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. Ersetzt man in der bekannten Näherungsgleichung für die Stromdichte 


0q ”»+1/q 2 ——— e 
ar 1=VaFRre. #2) 


- eq durch ou, so folgt aus (1) und aus den Gleichungen der Drehungsfreiheit für die dimensions- 


losen Störgeschwindigkeiten U —= = —1,V/= = Pe = das partielle nichtlineare Differen- 
tialgleichungssystem 1. Ordnung y ‘ 


oV oW- ou oV 2U oW >U 
Er Ger FRE, ER 3 (2) 


zur approxjmativen Beschreibung einer Strömung der betrachteten Art. 

Dieses System besitzt eine 14-parametrige Polynomiallösung. Beschränkt man sich aber 
auf Symmetrie sowohlzur Ebene y=0 als auch zur Ebene z=0 und verlangt man U=-V—W=0 
im Ursprung, so reduziert sich die Parameterzahl auf 4. Diese Lösung hat die Eigenschaft, daß 
im ganzen Feld U, = const = « gilt, und man kann noch auf «=1 normieren. Das damit dann 
noch wesentlich 3-parametrige Lösungssystem von (2) lautet!) 

„+ 


U-a+r+(%t —) u Veryt 2a + lin) y®-+3vya2, 


2 (3)- 
W=—pet @H1—2N get 374 rl _,)# 
Spezialisiert man hierin die Parameter A, u,» auf A — — ‚u=v=](, so entsteht der Fall 


einer schallnahen Strömung durch eine ebene Lavaldüse. Die kritische Isotache (Schallkurve) 

ist durch U=0, d.h. durch die Parabel + ı 
“+1 

6 4 


2 
Berg UNE vr ‚ so entsteht die schallnahe Strömung durch eine drehsymmetrische 
Lavaldüse. Die kritische Isotache ist das aus © + —— y?—=0 entstehende Rotationsparaboloid. 


y== 0 gegeben. Spezialisiert man in (3) auf 


48 et 
4 

Im allgemeinen Fall (3) ist die kritische Isotache durch die Fläche 2. Grad& 
zs+iy + = =) 220 gegeben. Für den Fall 0<A< as 
Halbraum x <0 verlaufendes elliptisches Paraboloid. Diesist eine in den prinzipiellen Zügen nichts 
Neues liefernde dreidimensionale Verallgemeinerung des zweidimensionalen ebenen bzw. dreh- 
| 1 
symmetrischen Problems. — Für den FallA<0 oder A> = dagegen ist die kritische Isotache 
ein hyperbolisches Paraboloid, eine Sattelfläche, und dies ist damit ein Fall, zu dem es kein 
zweidimensionales Änalogon gibt. Die Stromlinien sind im allgemeinen keine ebenen Kurven. 
Man kann sie noch auf mannigfache Weisen zu Stromflächen, d.h. Düsenkonturen zusammen- 
fassen, indem man beispielsweise einen Schnitt 2 —=const der Kontur beliebig (geschlossen) 
vorschreibt. In dem Bild wurde der Düsenschnitt x&—=0 als Ellipse vom Achsenverhältnis 2:3 
gewählt. 


ist sie ein ganz im negativen 


Die Stoßherzkurve bei allgemeiner Zustandsgleichung 
Von Hans Richter in Haltingen/Baden 


Bei der Behandlung von Stoßfrontreflexionen, -durchkreuzungen u.a. kommt es auf die 
Gleichheit von Drücken und Geschwindigkeitsrichtungen im Abströmgebiet an. Es empfiehlt 
sich daher die Verwendung der Herzkurve?) an Stelle der üblichen Stoßpolaren. 


Schreibt man die Hugoniot-Kurve in den dimensionslosen Variäblen &= (p— p})' = 


und B= 32 also u = u(£), so gilt für schiefe Stöße 


Y,' 
Ser : 1— 1 1 1 
BR EH E use ll) BREI Chan ta et. 0.5 0 40 (1) 
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= MWia%r=M—-E-l+u) .::2.2:.... Br See 


Sei Z. angew. EN oh? DER 
270 C. Strömungslehre 34.30 Nr. D Aue. Set. 1950 


und entsprechende Formeln für tg&, tg y usw. (1)ist die Gleichung der Herzkurve bei et 
Zustandsgleichung. Es soll untersucht werden, welche bei idealen Gasen bekannten Aussagen 
über die Gestalt der Herzkurve und das Verhalten der übrigen Stoßgrößen (%, y, ...) längs der 


Herzkurve bei Übergang zu allgemeinerer Zustandsgleichung erhalten bleiben, so daß Folge- 7 
rungen aus diesen Eigenschaften allgemeine qualitative Gültigkeit haben. - € 
1. Bei idealen Gasen hängt die Herzkurve nicht von p, und V,, sondern nur von M, ab. ’ 

Dies ist allgemeiner nur der Fall für Gase mit der Zustandsgleichung Re * 
U = const. VI*-f(s) . ER Ä 

U = innere Energie, f ()= — beliebige Funktion der ERS AR Si x 
2. Es wird nur vorausgesetzt, daß auf der Hugoniot-Kurve die En mit & wächst.  ” f 

Dies bedeutet, daß für senkrechte Stöße die Stoßstärke mit der Anströmgeschwindigkeit wächst, v 
was bei allen bekannten Zustandsgleichungen der Fall ist. e 
3. Die Herzkurve beginnt bil=0, 9—=0 mit ge — YM?—1/M! ‘und endet bei einem E 


en 


E = Cmax,; 9 — 0 mit senkrechter Tangente!). Es existiert mindestens ein Punkt (£„, 9m) maxi- 
maler Ablenkung. Für kleine M),— 1 gibt es stets genau ein 94; für beliebiges M} nur dann, 
wenn gilt: 


WRTRERT 2} 


d? Yu 1—u { | % 
TE Er SE. { 
(din? "ya | 9) 


also sicher bei Hugoniot-Kurven, die wie bei Abelscher Zustandsgleichung dauernd zur &-Achse 
konvex sind. 


4. Durchläuft man die Herzkurve in Richtung wachsender £, so wächst & monoton; % fällt 

zunächst und wächst dann; W, fällt monoton; M, beginnt mit M, >1 und endet mit einem 
Werte < 1. Es gibt also Stellen (£,, 9,) mit M,= =—=H, 

5. Für das anschauliche Verständnis der Labilität des starken Stoßes in einer konkaven. 
Ecke?) und für die Existenz der Crocco-Singularität®) ist entscheidend, daß gilt: <{„. Dies 
folgt aus der Voraussetzung, daß eine abgelöste Kopfwelle in dem Sinne stabil ist, daß infinitesi- 
Male Störungen nur eine infinitesimale Änderung der Konfiguration verursachen. Es zeigt sich 
nämlich, daß bei Anbringung einer infiritesimalen Störung im Bereiche M, >1 in kleiner, - 
jedoch endlicher Entfernung von der Grenze M,=1 ein endliches Stück der Kopfwelle 
sich aufsteilt und die Störquelle in den Bereich M, < 1 einbezogen wird, wenn £, >, wäre. 


6. Bei festen p, und V, hat die Gesamtheit aller Herzkurven eine eindeutige Einhüllende 


9-/i ET FR a ee er iD) 
mit zugehörigen Werten von M.: x Ar eo. = BER 
I+n | BR 
A, Z , £ 
hl en 


Die Herzkurven mit kleinen M,— 1 liegen zwiebelschalenförmig ineinander (bei Abelscher 


„Zustandsgleichung bis zu M) = /2); für M, > Y2treten stets Schnitte benachbarter Herzkurven }: 
‚auf. Die Elimination von £ aus (5) und (5a) liefert bei gegebenem halben Profilöffnungswinkel® | 
dasjenige M,, bei dem die anliegende Kopfwelle möglichst schwach ist; für kleine 0) tritt diesbi 
M,) »>2-+ const. d ein. EEE 
Expansion einer Gaskugel hohen Druckes 
Von Fr. Wecken in Haltingen 2 “ 

Zur Zeit =0 befinde sich eine Gaskugel vom Radius r, und unter hohem Imnendruck 
Pi >Ppa in ruhender Atmosphäre (u= 0, P = 9,); inr < r, seien bei t= 0 Anfangswerte p=p; ) 2. 
und w=u;(r) irgendwie vorgegeben (z.B. als Resultat einer chemischen Umsetzung); die. 7.4 
Entropie und chemische Zusammensetzung des Gases sei bei i=0 konstant. Dann ergibt sich für 
£>0 ein Ausgleichs-(Explosions-)Vorgang, der sich durch schrittweise Integration nach dem 
Charakteristikenverfahren®) unter Vernachlässigung von innerer Reibung, "Wärmeleitung Bu 3 a 
Difiusion berechnen läßt. Auf Grund einer für einen ‚bestimmaten Ba j hend Behagn 


1) Vgl. Bild2 a.a. 0, . rg. 
2) a. BR 

®) H. Richter: Ildistacco dell’onda anteriore sül profilo Denveesß: ae 
*) Kl. Oswatitsch: Z, angew. Math. Mech. 25/27 EN, 8. 5—20 
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numerischen Rechnung sowie theoretischer Untersuchungen lassen sich allgemeine qualitative 
Aussagen machen. 


1. Bei Punkt 1 (r=r,, t=0) beginnend erscheint der Bereich r>0, t>0 in 5 Gebiete 
geteilt: (A) = ungestörte Luft, (B) = komprimierte Luft, (C)= Gas unter Einfluß der Luft, 
(D) = expandierendes Gas, durch Außenluft nicht beeinflußt, (Z) = von der Expansion noch 
nicht betroffenes Gas. (D) umfaßt bald räumlich den größten Teil des Gases und erreicht einen 
Durchmesser d > r,. (D) läßt sich — ausgenommen die Abgrenzung gegen (©) — ohne Kenntnis 
der Zustandsgleichung der Luft isentropisch berechnen; die Expansion des Gases verläuft hier 
ebenso, wie es im Vakuum der Fall sein würde. Der Druck nimmt rasch ab; bei t— oo gilt für 


die Geschwindigkeit: u = =; für die Charakteristiken: = =uta, n == 2 Kan 1) 


a 


(a = Schallgeschwindigkeit, y = Adiabatenexponent). Ein einfacher Modellvorgang, bei dem 
diese Näherungsbeziehungen exakt gelten und außerdem p = const - f73? ist, wird zum Vergleich 
herangezogen und ermöglicht die Extrapolation für {— 00; hierbei empfiehlt sich eine Darstellung 
in einem Diagramm mit den Koordinaten r/t und 1/t. 

2. Grenzlinie zwischen (A) und (B) ist die Luftstoßwelle W,, zwischen (B) und (C) die 
Mediengrenze M, zwischen (C) und (D) imFalle regulärer Anfangswerte zunächst eine Charakteri- 
stik O,. Eine spezielle, der numerischen Rechnung zugrundegelegte, bei r = r, singuläre Anfangs- 
wertverteilung!) ergab dagegen anstatt C, eine bei Punkt 1 beginnende Stoßwelle W, als Grenze. 
Wegen der bei Punkt 1 vorliegenden Unstetigkeit für p und « war eine besondere Untersuchung 
des Verhaltens von p und win der Umgebung notwendig; als zweckmäßig erwies sich eine Reihen- 
entwicklung nach Potenzen von t (bei regulären Anfangswerten) bzw. von Yt mit Koeffizienten, 
Bl 


die von&= 2 und eventuell (schwach) von t abhingen. Als nullte Näherung ergab sich 


hier das vom analogen linearen Problem her bekannte Verhalten?), das man sich neuerdings in 
der ‚‚Stoßwellenröhre‘‘ zunutze macht °). 


Die numerische Rechnung bestätigt das Vorhandensein einer Stoßwelle W,, die späterhin 
sehr intensiv wurde. — Aber auch bei regulären Anfangswerten ist die Grenze von (C) und (D) 
im weiteren Verlauf stets eine Stoßwelle W,, da im anderen Falle die Überlegungen unter 1. das 
absurde Ergebnis liefern würden, daß bei t— oo überall p — Ogilt und das Gas den ganzen Raum 
erfüllt. W , verläuft anfangs auswärts und dicht hinter M und W, (jedoch einwärts relativ zum 
bewegten Gas), wird später bei r = 0 reflektiert und läuft dann hinter W, her. 

. 3. Die Grenze zwischen (D) und (E) ist eine Charakteristik O,, die alsbald (bei = t,, 
„Punkt 2°) den Mittelpunkt r = O erreicht. Hier sind die Zustandsgrößen stetig, nicht aber ihre 
Ableitungen; die bei Punkt 2 vorliegende mathematische Singularität wurde bereits von W. 
Hantzsche und H. Wendt, K.H. Weise sowie C. Heinz untersucht. Für die Umgebung von 
Punkt 2 wurden hier Entwicklungen nach 1/Inr bzw. 1/ln(t — t,) mit von (t— t,)/r abhängigen 
Koeffizienten verwendet; es wurde eine dritte Stoßwelle W, neu gefunden, die bei Punkt 2 be- 
ginnt, jedoch so schwach bleibt, daß sie bei der Integration vernachlässigt werden kann. Die 
numerischen Rechnungen sind noch nicht abgeschlossen und werden fortgesetzt. 


Über die Berechnung von kleinen stationären Störungen 
. in einer Semischten Überschall-Unterschallströmung 
mittels einer einheitlichen Charakteristikenmethode 


Von Manfred Schäfer 
(Max Planck-Institut für Strömungsforschung, Göttingen) 


Die Frage, in welcher Weise Störungen, die einer ebenen kompressiblen Strömung etwa 
- durch gewisse Bedingungen an einer festen Wand aufgezwungen werden, bei stationären Ver- 
hältnissen sich ausbreiten und dieses Strömungsfeld beeinflussen, ist rechnerisch bisher nur ver- 
einzelt (s. nächste Seite) t), 2), ®), *),°) angegriffen worden. Der Grund hierfür liegt zweifellos 
in der mathematischen Schwierigkeit, dieses Problem selbst unter den vereinfachenden Annahmen 
von Reibungs- und Wirbelfreiheit zu behandeln, da bekanntlich der Übergang von Unterschall- 
-u Überschallgeschwindigkeiten mit dem Wechsel der Strömungsdifferentialgleichung von ellip- 
tischem zu hyperbolischem Typus verbunden ist. 


3) G.Taylor: Proc. Roy. Soc. (A) 200 (1950), 8. 235— 247. 
- 2) H. Schardin: Physik. Z. 33 (1932), 8. 60-64. 
®) G. N, Patterson u.a,: Physic. Rev. 75 (1949), S. 1293, 
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Daß keine neuartigen Ergebnisse bei Beschränkung der Fragestellung auf reine Unterschall- 
strömungen zu erwarten sind, ist einleuchtend, da ja ihr Verhalten qualitativ dem der inkom- 
pressiblen Strömungen entspricht, und wird durch Vergleich mit einer einschlägigen Untersuchung. 
von A. Betz®) bestätigt. Andrerseits weiß man auch, daß sich in jeder Überschallströmung 
eine Störung längs einer sog. Charakteristik oder Machschen Welle fortpflanzt. Aber wie 2 
verhält sich nun die Störung, nachdem sie die Schallgrenze erreicht hat ? Die gegenseitige Beein- 
f lussung von Überschall- und Unterschallbereich erfordert, daß grundsätzlich von vornherein = 
das gesamte Strömungsfeld in die Betrachtung einbezogen werden muß und — von ganz 
speziellen Fällen ?) abgesehen — nicht etwa der Überschallbereich für sich allein studiert werden 
kann. 


Die Berechnung des Störverlaufes erfolgt im Unterschied zu einer früher» vom Verf. vor- - 
getragenen Methode ”) mittels direkter. Ansätze für die Störgeschwindigkeiten. Diese werden 3 
als Funktionen von „nebencharakteristischen‘“ Koordinaten betrachtet, die mit den geläufigen & 

„‚hauptcharakteristischen‘‘ Koordinaten des Überschallbereiches in einfachem Zusammenhang > 
3 
& 


stehen, ihnen gegenüber jedoch den wesentlichen Vorteil bieten, daß sie im ganzen Feld reell 
bleiben. Hinsichtlich der mathematischen Durchführung @ieser Idee muß auf die ausführliche 
Veröffentlichung 8) verwiesen werden; besonders angemerkt sei hier nur, daß für die Störfunk- 
tionen, die sehr einfachen Differentialgleichungen genügen, neben den ‚„‚Übergangsbedingungen‘“ 
an der Schallgrenze die Bedingung des Abklingens im Unendlichen der Strömungsebene zu for- 
dern ist. 


Für kleine Störungen liegt nun die Näherungsannahme nahe, daß die Charakteristiken IE 
der ungestörten Strömung unverändert für die gestörte Strömung beibehalten werden. Auf 
diese Weise kann für das von Taylor erstmalig behandelte Beispiel der periodisch gestörten 
ebenen Zirkulationsströmung !) ohne weiteres ein einfacher geschlossener Lösungsausdruck an- 
gegeben werden. Die volle Tragweite der Methode kommt bei der Behandlung des viel allge- 
meineren Problems zur Geltung, die Störausbreitung zu untersuchen, die durch ein im Überschall- 
bereich gelegenes hügelförmiges Hindernis auf der (krummlinigen) Wandkontur einer einseitig 
begrenzten Strömung hervorgerufen wird. Die Diskussion der Lösung zeigt, daß die Störung Bar 
sich längs der stromabwärts verlaufenden Machschen Welle unverändert fortpflanzt, ndn - 
Unterschallbereich sich dagegen diffus ausbreitet (damit hängt dort eine Wellung der Strom-_ 
linien zusammen) und rasch abklingt. Besonderes Interesse bietet der Fall, daß die Störung 
durch eine Krümmungssingularität der Wandkontur verursacht ist. Dann ergibt sich, daß die 
Geschwindigkeitsänderung längs der weiterhin glatt verlaufenden Berandung zwangsläufig ein 
singuläres Verhalten in gewissen Punkten stromabwärts aufweist, die sich an der Austrittsstelle 
des Überschallbereiches häufen. Die Lage dieser Punkte ist durch den geometrischen Verlauf 
der von der Störungsstelle ausgehenden und dann an Schallgrenze und Wandkontur wiederholt 
„reflektierten‘‘ Machschen Welle bestimmt. Entgegen demvonBusemann?),®)undGuder- 
ley?) gezogenen Schluß, daß die Häufung jener Singularitäten an der Austrittsecke im ll 
gemeinen ein unendliches Anwachsen der Störamplituden zur Folge haben müsse, wird gerade Ber 
die Stabilität der Störung streng ee | 2 a er 


& — PR 
L 


2 0. I. Taylor: Recent So on the flow of compressible tluids. I. Tondeit Mathem. Soc. Bd. 5 

- (1930), 8. 224—240. Deutsche Übersetzung in Z. angew. Math. Mech. Bd. 10 (1930), 8. 334—345. eg 

2) G.Guderle y: Die Ursache für das Auftreten von Verdichtungsstößen in gemischten Unterashall. EI 
UbeEPER EL ömungen (Monograph Transsonie Flow). Reports and Translations Nr. 110 British M. ©. B a FEN 
Olkenrode r r 
>) A. Busemannund G. Guderley: The Problem of Drag at High Subsonie Speeds., Reports, 

and Tr anslations Nr. 184. British M. O.S.(A) Völkenrode, März 1947. 

*) G. Guderley: Störungen in ebenen und rotationssymmetrischen Schall. und Überse 
strahlen. Z. angew. Math. Mech. Bd. 25/27 (1947), S. 190—195. 
EA rt Busemann: The Drag Problem at High Supersonic"Spoeds. 7: "Aoronanıt. 

6) A. Betz: Verlauf der eG a in der Nähe einer ' Wand bei 
der Krümmung. Luftfahrt-Forsch. Bd. 19 (1942), S. 129—131. > 
?)M. Schäfer: Über die Berechnung der ‚Ausbreitung von Säraiigent in eine 
nach der direkten Hodographenmethode. Vortragsauszug in Z. ber a Mir 
8) Erscheint demnächst in dieser Zeitschrift. 
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Überblick über Rechenhilfsmittel') 
- VonJ.C. P. Miller in London 


1. Einführung. Bedeutung, Notwendigkeit und Stellung der numerischen Rechnung in 
Mathematik und Technik. Auswirkungen der Entwicklung von Hochgeschwindigkeitsaufomaten. 
Rohe Einteilung der Rechenhilfsmittel. 

2. Erläuternde Beispiele. Herausstellung wesentlicher Gesichtspunkte bei der nume- 
rischen Rechnung an drei Beispielen: a) Auswertung eines Integrals, b) J. R. Aireys Verbesse- 
rungen der asymptotischen Entwicklungen bei ya uderiungeignen. c) Ausführung einer Division 
mit verschiedenen Hilfsmitteln. 

3. Mathematische Rechenhilfsmittel. Wichtigkeit Be neinatischer Methoden, Unter- 
scheidung zweier mathematischer Hilfsmittel: a) Spezielle Eigenschaften der betreffenden Funk- 
tion, b) Allgemeine Methoden wie Interpolation mit Differenzen usw. Sorgfältige Planung eines 
Problems. Anwendung mathematischen Scharfsinns. Beispiel: Berechnung einer Größe als 
Differenz zweier nahezu gleich großer Zahlen (Pendel). > 

4. Mathematische Tafeln. Großer Fortschritt in der Herstellung von Tafeln in unserem 
Jahrhundert. Zwei Richtungen: ein- oder zweibändige Sammelwerke und Erweiterung der 
tabellierten speziellen höheren Funktionen. 

a) Sammelwerke: «&) Elementare Funktionen. Siebenstellige Logarithmentafeln von 
Zahlen und trigonometrischen Funktionen, Barlows Tafeln, Chambers Tafeln, Fünfstellige Tafeln, 
Vierstellige Tafeln. 8) Höhere Funktionen. ‚Tafeln höherer Funktionen“ von Jahnke und 
Emde. Dreistellige Tafeln in une (Comrie-Miller). Tafeln von Hayashi, Anhang 

bei Peters und Stein. 

b) Tafeln spezieller Funktionen. Bedeutung und Leistung der BA. (British Asso- 
ciation, Mathematical Tables Committee 1871—1948, dann Royal Society Mathematical Tables 
Committee.) In Amerika NMTP (New York Mathematical Tables Project) und HCL (Harvard 
Computation Laboratories). 

Hervorhebung der Tafeln von Peters bei den speziellen elementaren Funktionen. Ent- 
wicklung der BA. Erwähnung älterer Arbeiten wie Bericht von Glaisher 1873, Faktor Tafeln 
von Glaisher und ‚‚Binary Canon‘ von Cunningham. Viele zahlentheoretische Werke. Gegen- 
wärtige Serie der BA.: Bd. Ibis IX mit zwei Teilbänden. 

Inhalt der Tafeln der NMTP: elementare Funktionen, Kreis-, Hyperbel-, Exponential- 
funktionen, Logarithmen, Reziprokwerte, Potenzen, Integralsinus, Integralkosinus, Exponential- 
integral, Wahrscheinlichkeitsintegrale, Besselfunktionen J,(z) und J; (z) für komplexes Argument, 
Kreis-Besselfunktionen, Funktionen mit gebrochenem Index +1/3, + 2/3, + 1/4, + 3/4 usw. 

HCL u.a.: J„(x) für n = 0(1)100, x = 0(0,001 oder 0,01)100 (13 Bände), Hankelfunktion 

der Ordnung 1/3 für komplexes Argument, verallgemeinerter Integralsinus und -Kosinus, inverse 
\  Hyperbelfunktion. 

‚  e) Reliefs und Schaubilder. Erleichterung für numerische Rechnung durch Reliefs, 
Nomogramme usw. Heranziehung bei Funktionen mit mehreren Parametern (Jahnke-Emde). 

5. Bemerkungen über die Ausführung mathematischer Tafeln. Hervorhebung 
gewisser drucktechnischer Gesichtspunkte, wie schlechte Lesbarkeit der Schreibmaschinen- 

typen usw. 

$ Tafelwerk nicht zu umfangreich, aber wenn möglich, lineare Interpolation. Kompromiß- 

lösung. Bedeutung und Empfehlung der Everettschen Formel mit modifizierten zweiten 
Differenzen. 

6. Rechenmaschinen und mathematische Maschinen. a) Instrumente. Gra- 

phische Darstellung, Rechenschieber, Planimeter, Integraph usw. TEE RESN Lösung 
'simultaner Gleichungen nach Mallock. Beschränkte Genauigkeit. 

‚ b) Ziffernmaschinen. Keine vollständige Aufzählung möglich. Gewöhnliche Rechen- 

maschinen: &%) Addiermaschinen, 8) Vier-Spezies-Maschinen. National-Buchungsmaschine. 
ET Hollerith- und andere Lochkartenmaschinen. Große Ziffernrechenmaschinen. 

BR 7. Allgemeine Bemerkungen. Keine Verdrängung der gegenwärtigen Methoden durch 
die großen Rechenautomaten. Im Gegenteil erhöhte Forderung nach ausgebildeten Rechen- 
Sn kräften, mathematischen Tafeln, gewöhnlichen Rechenmaschinen usw. Mathematische Tafeln 

” 2 als, „Gedächtnis“ für Maschine. Rechenautomaten müssen Lesen lernen. Vorläufige Verwendung 

= ‚von großen Rechenautomaten: Behandlung von Problemen, deren numerischer Aufwand nur mit 
sen großen Maschinen überhaupt bewältigt werden kann. 


1) Nach dem Manuskript von J. ©. P. Miller vorgetragen von H. Unger, Darmstadt. Es sei auf die 
\ che ‚Arbeit von Dr. J. ©. P. Miller verwiesen: General Survey in Modern Methods of Computation, 
h Kae Bi: The ‚Adrancement of Science 1948, Vol. V, Nr. 18, S. 79— 85. 
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Interpolation im Komplexen 
Von W. Schnell in Darmstadt 


Die Berechnung komplexer Funktionswerte, die ständig zunehmende Bedeutung erhält, 
ist meist recht schwierig, da nur wenige Tafeln von Funktionen komplexer Veränderlicher greifbar 
sind und diese sich wegen der doppelten Mannigfaltigkeit des komplexen Arguments mit großen 
Argumentschritten begnügen müssen. Ihre volle Ausnutzung erfordert eine Interpolation höherer 
Ordnung, für die sich drei Möglichkeiten ergeben. 

1. Methode: Durch Aufspalten der Funktion w = f(z) der komplexen Veränderlichen 
z—= x +iyin w=u-+- iv werden Realteil vu und Imaginärteil v reelle Funktionen von & und 
y, womit unmittelbar die reellen Interpolationsformeln für zwei Veränderliche übernommen 
werden können. Dann erfordert aber die Bestimmung eines Zwischenwertes bei Verwendung 
der Formel von Everett!) unter Mitnahme der zweiten Differenzen 44 Multiplikationen und 
30 Additionen. Nach Steffensen?)kann man diese Rechenoperationen zwar einigermaßen 
übersichtlich anordnen, doch bleibt dieser Weg für die komplexe Aufgabe ungeeignet. Durch 
das Aufspalten wird der wesentliche Charakter einer analytischen Funktion zerstört ?) (Cauchy- 
Riemannsche Differentialgleichungen, Laplace-Gleichung zeichnen Real- und Imaginärteil vor 
gewöhnlichen Funktionen von zwei Veränderlichen aus). 

2. Methode: Diese grundlegenden funktionentheorietschen Ergebnisse lassen sich aber 
in die reelle Aussage einarbeiten. Hierzu ersetzen wir in der zweidimensionalen reellen Newton- 
formel (Stützstellen auf einem Rechteckgitter mit Spanne h bzw. kin x bzw. %) 


m. 22 ; ; = R a 
F(p, g)= 35 > (?) (Ya A (0:0 mit p=°l ge ; Yon) 
ys=gıN 


u=0 


in der die zu interpolierende Funktion F nur die Zusammenfassung zweier reeller Funktionen 
U(p, q) und ©: V(p, q) bedeuten soll (?ist zunächst nur zusätzliche Kennzeichnung von V), alle 
Differenzen nach y und gemischte Differenzen nach x und y durch reine Differenzen nach x. 
Durch wiederholtes Differenzieren der reellen eindimensionalen Newtonformel kann man nämlich 
die Zusammenhänge zwischen Ableitungen und Differenzen bestimmen und damit die Differential- 
aussage vn Cauchy-Riemann und alle ihre Folgerungen in DIE angeben. 
Setzt man diese EB — die ersten lauten speziell für A= k 


4A,U=—A ‚/+5 l ay_ —4: U AV SAU „U = —- AU+HAU+LV 


A vn Sara. AsA, V=-A4U—2 RU AV.. 


usw. für höhere Differenzen — in die reelle Doppelsumme ein und betrachtet jetzt valsi auaemar 
Einheit, so entsteht eine DE komplexe Newtonformel | 


RO) D()ARO+ Run, 


=0 


die sich von gr reellen nur dadurch unterscheidet, daß jetzt ein komplexer Zuwachs des Argu- - 


a mE 
ments = —— auftritt. Wenn bei dieser Umformung die äquidistanten Stützstellen speziell 


auf einer Parallelen zur x-Achse liegen, so läßt sich aber dieselbe Überlegung auch für die y-Achse 

und darüber hinaus für jede beliebige andere Richtung durchführen, wobei dann die Spanne & 

selbst auch komplex wird und die Differenzen längs dieser Stützgeraden zu bilden sind. 
Dieselbe komplexe Interpolationsformel 


l 
Fo = I6)* F0)+ Ryy mit I= 


erhält man bekanntlich auch rein funktionentheoretisch über die Identität von Nico 1 e; doch 
ergibt sich bei dieser e leganteren Ableitung nicht ZWDR SSL eine Aussage über die Anordnung 
der Stützstellen, wie wir sie oben fanden. 
Besondere Bedeutung für den praktischen Rechner erhält diese Formel oder noch besser ER a 
— wegen der günstigeren Lage der Stützstellen — die entsprechende Erweiterung der Everett- 77 
') s. z.B. G. Schulz: Formelsammlung zur praktischen Mathematik. is Göschen us, 
Berlin 1945, 8.83. Re: 


’) J.F. Steffensen: Interpolation. London 1927, 178—223. AN, > Bu r 


°) A. Selow: Zur Berechnung und Anwendung der ynderinktin vom Index E ” ur kom- SR 3 
plexem Argument. Diss. Darmstadt 1949, 36—44, Ne 


u ch 
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formel dann, wenn man ihre Koeffizienten für eine Zehntel!) oder Hundertstelteilung im voraus 
vertafelt, da dann der gesuchte Zwischenwert mit 4 reellen und 4 komplexen, insgesamt gleich 
20 gewöhnlichen Multiplikationen und 6 Additionen zu finden ist, was eine Verminderung des 
Rechenaufwandes gegenüber der reellen Methode um über die Hälfte bedeutet. 

Auch läßt sich mit dieser Methode bei Aufstellung geeigneter Koeffiziententafeln die Inter- 
polation von Funktionen, die auf Strahlen vorgegeben sind (Argument 2—=r- e!P), verhältnis- 
mäßig leicht durchführen. 

3. Methode: Da allerdings bei Mitnahme höherer Differenzen nach der letzten Methode 
Funktionswerte in die Rechnung eingehen, die weit vom Interpolationspunkt entfernt liegen und 


- daher nicht dem eigentlichen Funktionsverlauf in der Nachbarschaft des interessierenden Punktes 


genügen, kann man eine weitere Vereinfachung erwarten, wenn man die Stützstellen statt auf 
eine Gerade in die Ecken eines Polygons legt. 

Koeffizienten der aus dieser Überlegung entstehenden Interpolationspolynome, die man 
auf den verschiedensten Wegen finden kann — z.B. komplexe Lagrangeformel; Ansatz mit 
unbestimmten Koeffizienten und Lösen des Gleichungssystems, das sich dann mit den bekannten 
Stützwerten ergibt usw. — und die sich je nach Wahl der Stützstellen unterscheiden, sind in den 
unten genannten Arbeiten z. T. vertafelt?)?), z. T. als Funktionen des komplexen Zuwachses 
angegeben?). 

Zusammenfassend sieht man, daß eine Interpolation im Komplexen — wie jede komplexe 
Rechnung — einen ziemlichen Aufwand erfordert. Trotzdem kann gesagt werden, daß sich nach 
den beiden letzten Methoden wesentliche Vereinfachungen erreichen lassen, besonders wenn man 
sich möglichst engmaschige Koeffiziententafeln im voraus zusammenstellt. Dabei gilt diese 
Bevorzugung der wirklich komplexen Rechnung nur zur Ermittlung einzelner, zerstreut liegender 
Zwischenwerte, während eine systematische Untertafelung auch im Komplexen durch Auf- 
spaltung in Real -und Imaginärteil mit den reellen Methoden, vor allem mit Endziffer- und 
Brückenverfahren nach Comrie, am schnellsten zum Ziel führt. 


Grundsätzliche Probleme der Abrundungsfehler 
Von A. van Wijngaarden in Amsterdam 


Untersucht wird die Frequenz des Auftretens von Fehlern vorgegebener Größe im Er- 

gebnis gewisser systematischen Rechnungen mit abgerundeten Zahlenfolgen. Durch passende 
Definition von Funktionen kann jede abgerundete Zahl als ganz betrachtet werden. Der ab- 
gerundete Wert Af und der Bruchteil Bf einer reellen Zahl f sind dann definiert durch 
Af+Bf=f; Af=0 (mod 1), |Bfl<1/2 wenn f=#1/2 (mod 1) und Af=0 (mod 2), |Bfl=1/2 
wenn f=1/2 (mod 1). Es liege eine große Zahlenfolge f;(k=1, 2,...) vor, und es sei 
f=F(f»..., fj+n-ı) zu bestimmen. 
Werden statt f; die Af; benützt und statt F, die etwa nur äbgerundet verwendbare Kon- 
stanten enthalten kann, eine ‚‚abgerundete‘“ Funktion AF, dann ist das Ergebnis der Rechnung 
g=AF(Af, ..., Afj+n-ı). Von der Diskrepanz y=f—g kann durch Wiederholung für ver- 
schiedene 5 die Verteilung v(y) definiert werden. Manchmal ist nur Af zu bestimmen. Die Rech- 
nung liefert Ag und die Frequenz der Diskrepanz P=Af—Ag sei 2(P). 

Wichtig für v(y) und Q(P) ist die. Verteilung und die Abhängigkeit der B/f,. Sind die 
f Werte für »—=kh einer Funktion f(x), dann tritt meistens praktisch Gleichverteilung der Bf, 
auf, was weiter auch vorausgesetzt wird. Ist f(x) n-fach stetig differenzierbar, dann ist bekannt- 
lich Ar f=Mf® (£), mit u <&< my. Istf"(x) beschränkt im betrachteten Bereich dann 
sind bei genügend kleiner h die 4” f vernachlässigbar. Allgemeiner sei die ‚kritische‘ Differenz 
Amf die Differenz niedrigster Ordnung deren Bruchteil praktisch konstant ist. Dann ist 


Bin = B(BAnf,— 3 (—1)mer1 (v)BR). also B/m+ı praktisch funktionsgemäß ab- 
k=1 


‚hängig von seinen m Vorgängern. Außerdem wird vorausgesetzt, daß weniger als m aufeinander- 


‚folgende Bf, unabhängig sind. 


< n 
Für den wichtigen Fall der homogenen Linearkombination f= > a; f; mit exakten a, und 
kl 


n<m wird v(y) mit Hilfe der Laplacetransformation bestimmt, während die Bestimmung von 
@(P) einfache zahlentheoretische Mittel erfordert. Für n >m werden v(y) und 2(P) durch mehr- 
dimensionale geometrische Betrachtungen gefunden. Im Ergebnis für 2(P) tritt für n<m 
2) W. Schnell: Interpolation und Untertafelung im Komplexen. Dipl.-Arb. Darmstadt 1950. 
2) A.N.Lowan und H.E. Salzer: Coefficients for interpolation within a square grid in the com- 


plex plane. J. Math. Phys. 23 (1944), 156—159. 
3) H.E.Salzer: Formulas for direct and inverse Interpolation of a complex function tabulated along 


_ quidistant Cireular Arcs. J. Math. Phys. 24 (1945), 141—142. 
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der kleinst mögliche gemeinschaftliche Nenner q ei a; auf. Bei kontinkerHrBeR Variation der 
a, verhält 4 Q(P) sich also ganz pathologisch. Be: 
Die Theorie ist auf die Bestimmung der Diskrepanz P=A4A"f— Ar Af andewande Eszeigt 
sich, daß außerhalb der Ordnung n auch der Wert n—m eine wichtige Rolle spielt. Ist n <m 
die Differenz also subkritisch, dann ist 2 (P) unabhängig von n und BAmf. Ist dagegen. nZm, 
die Differenz also kritisch oder superkritisch, so hängt 2(P) von n und BAmf ab. 
Eine zweite Anwendung liefern Interpolation, numerische Differentation oder Integra- 
tion in einem beschränkten Bereich. Sei z.B. J,(f4, -.., fj-rn-ı) der interpolierte' Wert. für 


M=EHh =] +»), dann hat die Frequenz Q(P) von P= Aly(fu ..., fi+n-ı) 


— Al, (Afp:- .,Afjtn-ı) eine Unstetigkeitsstelle für jeden rationalen Wert von p wenn n<m. 
So ist z.B. für p=0 offenbar 2 (0)—=1, während für p=0-e,2 (0)= 3/4 für Interpolation 
beliebiger Ordnung gilt. Für n>m liegt die Sache ganz anders. Speziell wenn n>m ist, 
z.B. bei Integration über einem langen Bereich treten ganz besondere Gesichtspunkte auf. 


Grundsätzliches zur numerischen Integration von Sewöhnlichen 
Differentialgleichungen 


Von W. Quade in Hannover 


sg. 


FIR WITTR BL ur 


Von den in den letzten Jahrzehnten entwickelten Verfahren zur angenäherten numerischen 
Integration von gewöhnlichen Differentialgleichungen sollen hier vornehmlich die der an- 
genäherten Quadratur mit nachfolgender Iteration betrachtet werden (Verfahren 
von Adams). Dabei möge zunächst nur die explizite Differentialgleichung erster Ordnung 
y =f(x, y) in Betracht gezogen werden. Ist y=u(x) ein durch den Punkt a, b gehendes 
Integral einer solchen Differentialgleichung, so besteht bei stetigem f(x, y) die Identität 


uatfnuwla. 


Den vorerwähnten Verfahren liegt der Gedanke zugrunde, das Integral rechter Hand durch an- 
genäherte Quadratur auszuwerten. Wird dabei der Inbegrand durch ein Kobrupe 2) an- 
genähert, so ist die Näherungslösung 


y=o(a) =b+ | Piyat 


MVRRIUTLTE, Br BE 


wieder ein Polynom. 3 
Zur Beurteilung der Annäherung von v(x) an u(x) dient die‘ folgende allge meine) Fehler- 

abschätzung!). In dem Gebietsstreifen & a 
G: asssa +0, — o<y<o (a>0) | a 

sei f(x, y) stetig und genüge dort einer Lipschitzbedingung mit der Konstanten L. Ferner be Sarg 
die durch a, b gehende N y=v(z) in<a,a-+ > ‚eine beschränkte, au Eee es 
bare Ableitung, und es werde TEN j Sa 
v’(2) — [a v(@)] = r(@) 


als der mit v(x) gebildete Rest der Differentialgleichung y’ = —=/®, 9 Pe Dann gi Ta 
< 4,0 -+- > die folgende allgemeine a N 


u) —v@)| < fer Ir@ldt. nn 


Bei den eingangs erwähnten Verfahren wird das Intervall <a, a+ er durch eine Zer- 
legung EEG ea 


ad zn Sera Br 


in n gleiche Teilintervalle von der Lange a S _ „ aufgeteilt; ferner wirdo(a Yı und f(% » Y%) 
= = gesetzt. EEE 


en der ersten Stufe, u Se 

Bei diesen wird das Polynom P(&) ausschließlich auf Grund 
stimmt, Hierher gehört u.a. das Verfahren von Euler- Cauchy, - 
Als Vertreter eines Interpolationsverfahrens möge hier nur der folgen 
1. P(x) ist ein Polynom vom zweiten Grade (Si mp: on) 

im Intervall < 2, %> durch 


P(&)=fo, - Pia) fi 
5) Vgl. hierzu auch W. Quade: Math. Z. 48 (1942), 
eine Fehlerabschätzung bei G. Hamel: Z, angew. Math, 
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festgelegt. Durch Quadratur ergibt sich 
AıT%= 75 Bft Se) > 


n-n= (ft +4ı+f2 


NY we zweite dieser Gleichungen ist mit der Simpsonregel identisch. Wird h so klein gewählt, daß 


(1). 


3 SL < list, dann können aus (1) die beiden Unbekannten y, und y, Baer Iteration bestimmt 


el. Der Rest r(x) ist im Intervall <x,, 2,> stetig und besitzt in den Punkten z,, &,, &, 
jeweils eine einfache Nullstelle. 

Der hierbei hervorkommende Wert von y, wird bei weitem ungenauer als der von %,. Diese 

.  Ungenauigkeit läßt sich vermindern, wenn im Intervall <&,, &> und nur in diesem die halbe 
Schrittweite verwendet wird (Anlaufrechnung), und im weiteren Verlauf der Rechnung nur noch 
die zweite der Gleichungen (1) benutzt wird. 

Interpolation mit Polynomen von höherem als dem zweiten Grade bringt im allgemeinen 
trotz vermehrten Aufwandes an Rechenarbeit keinen erheblichen Gewinn an Genauigkeit. Eine 
wirksame Steigerung der letzteren ist möglich 1. durch Verwendung einer kleineren Schrittweite, 
2. durch Benutzung der vollständigen Ableitung von f(x, y). Die zweite Möglichkeit führt 
auf die 

Verfahren der zweiten Stufe. 


Bei diesen wird vorausgesetzt, daß f(x, y) in @ stetige partielle Ableitungen besitzt, und 
daß a. vollständi ge ER 


2 fa y)=f (04) +f, (8 y) F®, y) a y) 


in @ einer ee mit der Konstanten L’ genügt. Bei. der Bestimmung des Poly- 
noms P(x) werden jetzt außer den Funktionswerten f; auch noch Werte g, herangezogen. Hier- 


her gehören die Verfahren von Duffing?) und Pflanz?). Bei weitem genauer als diese sind die _ 


% folgenden neuen Verfahren. 


- 2. P(x)ist ein Polynom vom dritten Grade (Euler-Maclaurin). 
P(x) wird im Intervall <x,, &,> durch 


Pia) = fo. P' (2) = 90, P(a)=fi P(a)=g9 

festgelegt. Hieraus folgt | 
| y-n=3(th-5 0-9). 

welches die ‚Euler- Maclaurinsche ee ist. Aus der letzten Gleichung kann y, 

: durch Iteration bestimmt werden, wenn Z 2 L Be 7 Lv < list. Der Rest r (x) besitzt in <z,, &,> 


_ eine stetige Ableitung und hat in x, > x, jeweils eine doppelte Nullstelle. 
FRE = 3. P(a) ist ein Polynom vom dritten Grade (Simpson verbessert). 
a Pia) wird im Intervall <E &,> durch 
Ben ea A) =f,, .. Balz: Po)=f 
> festgelegt, was auf; 
er 2% en DR rg UF ® 9% “ 6 0 + 20f, A: 


y— ren? 


i — 


pin Be): 1, 2) 


ge . Ing.-Wesen 224 (1920), 8. 2950. 
ih. Mech. 28 2. TS 107-172. | 


. % , are + 8% 
ie. 2 ie tla We Ky3: 
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bestimmt (Hermitesche Interpolation). Dies führt auf die Formeln 


h RAT: 
Yi7Z Ya MFH 128, +11%)+ 50 13 %— 409 — 39); 


h h? 
Ya 7 Io as YA+16f+7R) ER: (92 — 90) - 

A 1 
Hieraus können y, und ,y, durch Iteration bestimmt werden, wenn F hL-+ 3 h®L'<1 ist.: Der ? 
Rest r(x) besitzt im Intervall <&,, &,> eine stetige Ableitung und in jedem der Punkte x), &1, %, 
eine doppelte Nullstelle. Fa 

Eine weitere Erhöhung der Genauigkeit ist durch die Verfahren höherer Stufe möglich. 
Die hier für die explizite Differentialgleichung erster Ordnung angestellten Überlegungen lassen 
sich sowohl auf Systeme von Differentialgleichungen erster Ordnung als auch auf Differential- 
gleichungen höherer Ordnung ausdehnen. 


Bemerkung zur Theorie”der Formeln von Runge und Kutta 
Von E. Weinel in Jena 


Der Zugang zu einer allgemeinen Theorie der Formeln von Runge und Kutta ist bei der 
üblichen Art der Herleitung bekanntlich sehr erschwert durch die umfangreichen Rechnungen, 
die sich beim Angleich der Näherung an die Taylorentwicklung der genauen Lösung ergeben. 

In Vorlesungen beschränkt man sich daher gewöhnlich auf die Verifikation der besonders 
bequemen Formeln von Kutta!) und der entsprechenden Formelsätze von Nyström?) 
und Zurmühl?) für Differentialgleichungen höherer Ordnung. 

Indessen läßt sich — und das ist der Gegenstand dieser Mitteilung — die Theorie der 
gesamten Runge-Kuttaschen Formelmannigfaltigkeit unter einem einheitlichen Gesichtspunkt 
verhältnismäßig einfach entwickeln, wenn man statt der Taylorschen Formel einen Interpola- 
tionsausdruck zugrunde legt und die dem Runge-Kuttaschen Verfahren eigentümliche rekursive 
Rechenvorschrift (,‚Schrittrechnung‘‘) als eine geeignet korrigierte Euler-Cauchysche Näherung 
deutet. Aus diesem Prinzip ergibt sich zugleich eine Darstellung des Integrationsfehlers und der 
Fehlerfortpflanzung im Zuge der Rechnung. , ZEN, 


Zur Fehlerabschätzung bei der numerischen Integration (Differenzen- 
verfahren) gewöhnlicher Differentialgleichungen 

Von J. Weissinger in Hamburg | > 2 

Ein ausführlicher Bericht wird in einem späteren Heft dieser Zeitschrift erscheinen. u 


Zur Herleitung von Konvergenzkriterien für Iterationsverfahren 
bei linearen Gleichungssystemen | 


Von L. Collatz in Hannover 


In den letzten Jahrzehnten sind verschiedene Iterationsverfahren für lineare Gleichung- 
systeme aufgestellt worden, von denen die bekanntesten wohl die Verfahren in Gesamtschritten 
und in Einzelschritten sind®). Für die Verfahren wurden eine Reihe gewöhnlich hinreichender 
Kriterien aufgestellt und zwar wurde in der Regel der Beweis jedesmal eigens für das betreffende = © 
Kriterium neu erbracht. Man kann aber durch Zugrundelegung der genauen, d.h. der zugleich 
notwendigen und hinreichenden Kriterien und unter Verwendung der Sätze der Matrizenlehre 
die speziellen hinreichenden Kriterien auf einheitliche und kurze Weise (als Abschwächungen der _ 
genauen Kriterien) herleiten und dabei einige Kriterien zugleich etwas erweitern und auf eine 


für die Anwendungen bequeme Form bringen. SE Er NE 


')M.W. Kutta: Z. Math. u, Physik 46 (1901), 8: da. 2 
2) E.J. Nyström: Acta‘Soc.’ Sci. ;fenn..00.(1920), Zr = re 

®) R. Zurmühl: Z. angew. Math. Mech. 28 (1948), 8.13. ER EN, 
r 4) R. v. Mises und H. Pollaczek-Geiringer: Praktische Verfahre 
7. angew. Maöh. Mech. 9 (1929), 8. 58—77 und 152—164. — Weitere Literatur z . 
die Lösung von linearen Gleichungssystemen durch Iteration. Z. angew. Math. M di 
I. Cesari: Sulla risoluzioni dei sistemi di equazioni lineari per approssime ces; 
Accademia Nazionale dei Lincei, Classe Seienze fis, math. natur, vol. XXV, ser. 
E. Bodewig: Bericht über die verschiedenen Methoden zur Lösung eines Systems liı 
reellen Koeffizienten I-V. Koninklijke, nederlandsche Akademie appen 
1104—1116, 1285—1295, (1948), 8. 53—64, 211— 219. ER, TR 
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Vorgelegt sei das lineare Gleichungssystem 
Br ki=r, 

wobei A = (a;;) die gegebene (quadratische, n-reihige) Koeffizientenmatrix, x die Spaltenmatrix 
der Unbekannten &,,..., x, und r die gegebene Spaltenmatrix der rechten Seite.r,, u, Pa selk 
Es sei det Y #0. Nun wird Yin der Frm XA= 8-6 Beben und von einem willkürlichen 
Vektor x(%) ausgehend iterativ die Folge von Vektoren x, r'@),... nach 

Ber LEN —=r RU 
ermittelt, wobei det 8=+0 vorausgesetzt wird. Das Iterationsverfahren ist bei beliebigem Aus- 
nr x(°) dann und nur dann konvergent, wenn alle Wurzeln x der Gleichung 

det («B+&)= 0 

Beträge kleiner als 1 haben. 
Nimmt man die Hauptdiagonalelemente von X in die Matrix ® = (b;,), die anderen in G, 


also 
h; _ j@jr fürrerh "4 Oelür 27% 
Fell tor ek =, für j=#+% 
so hat man das Verfahren in ‚‚Gesamtschritten‘‘, für 
. RE Ak für N) DE | easlur,’ 7% 
Our <<. & \9 2, für.4<k 


erhält man das Verfahren in „Einzelschritten“; bei beiden Verfahren seien alle a,; #0. 

Dann ist also für die Konvergenz des Gesamtschrittverfahrens, bzw. des Einzelschritt- 
verfahrens, bei beliebigem Ausgangsvektor x(P) notwendig und hinreichend, daß alle Wurzeln x 
von 


%dıı dia at0.e din ä x 41 dis ... din 
dsı Hlsg dgn Do RAR % (sı # (digg nahe don Sa) 
RR ERTL HER ER ER 26 % Ayn 


Beträge kleiner als 1 haben. 

Nun kann man leicht die bekannten hinreichenden Kriterien verifizieren, z.B. sieht man 
durch eine einfache Umformung, daß bei einer reellen symmetrischen (allgemeiner : hermiteschen) 
positiv definiten Matrix X alle Wurzeln x der obigen rechtsstehenden Gleichung Beträge kleiner 
als 1 haben, daß dann also das Einzelschrittverfahren konvergiert; ist das ‚Zeilensummen- 
kriterium‘‘ 


n A e 
za < lan) G=1...,n) 
rd x 
oder das BOB NER SBEtUNG | 
Pal < |ar|' (Keil .sn) 
Find: - 


erfüllt, so sieht man ebenfalls leicht, daß dann alle Wurzeln x beider obigen Gleichungen Beträge 
kleiner als 1 haben, daß dann also Gesamtschritt- und Einzelschrittverfahren konvergieren. 
- Auch andere Kriterien, wie z.B. das E. Schmidtsche Kriterium 


Ay; |? 
| x S- = 1 
FR iz a Ei] 
lassen sich bequem erfassen. [ 
222 ‚Est bei einer Matrix A ? > 
: aeTm Er. ? ©; <lay| für delt, 
- 5 be : £ „für Butene ein j 


für a 
für mindestens ein k. 
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Diese abgeschwächten Kriterien reichen noch nicht für Konvergenz des Gesamtschritt- F% 


oder Einzelschrittverfahrens aus, wie .das Beispiel der Matrix ; = A 

Ze 2 > RE 

0 1 1 ee 

Be 1 a 

ee dazu braucht man noch, daß X nicht „‚zerfällt‘‘. Wenn man bei einer Matrix X die Indizes 
1,...,ns0 Mit 01... Op 01 +: On, (mit 1 Sr Sn— 1) numerieren kann, daß £ 


a, >0 istfü v=el..,r; vwl..,nr 


so heißt X „‚zerfallend“. Dann gilt: Ist die Matrix X nicht zerfallend und außerdem das schwache 
Zeilen- oder das schwache Spaltensummenkriterium erfüllt, so konvergieren Gesamtschritt-und 
Einzelschrittverfahren. Dieser Satz läßt sich bei Anwendungsbeispielen oft benutzen; z.B. bei Bi 
genäherter Lösung von linearen Randwertaufgaben bei gewöhnlichen und partiellen Differential- Fr 
gleichungen mit dem Differenzenverfahren kommt man häufig auf Gleichungssysteme, die sich 
dem eben genannten Typ unterordnen. Die Konvergenz des Liebmannschen Mittelungs- 
verfahrens bei der ersten Randwertaufgabe der Potentialtheorie ergibt sich so als Spezialfall 
unseres allgemeineren Satzes. A 

Eine ausführlichere Darstellung mit Beweisen soll in der Mathematischen Zeitschrift er- 
scheinen. 


Zur Iteration bei linearen Gleichungen 


Von H. Sassenfeld in u 
Für das Gleichungssystem 


mit der Iterationsvorschrift i Sr 5 a 
"a Yırd LEW a0 1 Re er EN Ze 2 
Br strebt die Folge x für jeden Anfangsvektor x gegen den Lösungsvektor r, wenn n die charakte- 
ristischen Zahlen A von = B-1@ im Einheitskreis liegen. Schreiben wir mit einer His 

matrix ®, deren charakteristischen Zahlen=+# —1 seien, die Identität Br u 
CHSBPLHFE-BDIET 60 7 a 
so erhalten wir entsprechend (1) und (2) die Iterationsvorschrift | EEE 
SBLOLFHE-BPMPLHR=0 2.2.2.2... 


‘Das neue Iterationsverfahren (4) konvergiert immer, wenn die charakteristischen Zahlen vn = 
3+r— (B+€E)1(8— $P) im Einheitskreis liegen. Durch passende Wahl der Hilfsmatrix 2 
läßt sich also der Wirkungsbereich eines Iterationsverfahrens von der Form (2) erweitern, wozu RR 
gegebenenfalls BP = ®®) als vom Iterationsschritt abhängig genommen werden kann. er a = Br 


Für ® wird man zweckmäßigerweise eine einfache Matrix, etwa eine ea 
wählen. So lautet für die Skalarmatrix B = € die verallgemeinerte Iterationsvorschrift 


Te (HB HC-aB)P+a=0 .... 


Vergleicht man für,z(0) = z(9 die Näherungsvektoren wu und x”, so ergibt sich | 


hang e 
n 2 z ENT T 2 
ar Cie Fu? - Ir 


d.h. die 2 sind die Eulersche Transformation der ww m 
index at). Aus I* — TE e 
al). Aus I*= (B+ E)1(8 dr a =. Igt u 


ristischen Zahlen x von 3* und A von 8 in 1 


I IR I | 
und es gilt: Zu jeder Iterationsmatrix 8, für ee. 
charakteristischen Zahlen A enthält, den Punkt A —= 
Werte a, für welche die charakteristischen Zahlen 

!) Vgl. K.Knopp, Theorie und Anwendung de 


berg (Springe -Verlag) 1947, 863. In der Reihenlehre ; 
2? —1l=x steht. 


Paz fr VER 
N fe 
Z. angew. Math. Mech, g we, 3 ; 
Bd.30 Nr.8/9 Aug./Sept. 1950 . D. Praktische Analysis 81 
FPF]F]7]€7€— € 1 — 


solche «& ist das Verfahren konvergent. Das ursprüngliche Konvergenzgebiet Al<1 wird dann 


"und nur dann in |<|<1 überführt, wenn « reell und nieht negativ ist. 


Man erhält aus (5) ein verallgemeinertes Gesamtschrittverfahren, 
wenn man dort ®'gleich der Diagonalmatrix D®= (a,;) setzt: 


en 1 . 
wi = ze DIEN HF DIA) +] a un 29 


Diese Gleichung zeigt deutlich den geringen zusätzlichen Rechenaufwand gegenüber dem 


‚gewöhnlichen Verfahren, das sich für « = Oergibt. Für (8)ist durch passende Wahl von «immer 


Konvergenz erzwingbar, wenn die charakteristischen Zahlen u von ®-1X in der rechten Halb- 
ebene liegen (insbesondere also, wenn D-1X- positiv definit ist). Es läßt sich das folgende ein- 


 fache hinreichende Konvergenzkriterium, das allerdings nur charakteristische Zahlen im Strei- 


fen 0<NReu<2 erfassen kann, aussprechen: Ist 
Sr Gr | pi 


SaN.<1 und Max 
ki Gi Mk a 


ki ki 


5 ME 3.220), 


= Max 


‘so konvergiert für M<1 das gewöhnliche Gesamtschrittverfahren, für M 21 das Verfahren 


M’—1 
2(1—N) 
Bei der durch (9) vermittelten Schätzung für halbe Abszisse und Betrag der charakteristischen 
Zahlen A ist am günstigsten MN 


Be A TR DIE 
Praktisch gibt man & einen runden Wert zwischen den Werten (10) und (11). 


(8) sicher für jedes reelle 


> 822 (10): 


Nichtlineare Behandlung von Eigenwertaufgaben 
Von A. Unger in Darmstadt 
1. Newtonsches Verfahren!) bei einfachen charakteristischen Zahlen 
Gegeben sei das lineare spezielle Eigenwertproblem 
ee A-FAHL-NAr-=0 mt NA)=mmy..s MM) 22er: (1). 
Gesucht: Die charakteristische Zahl A; mit dem zugehörigen Eigenvektor ı; = {2,, &,, . . ., Zn}: 
Setzt man 2, —1 (gegebenenfalls eine andere Komponente), dann sind die n Komponenten des 


Vektors D = {A;, &, %, ..., 2,} unbekannt. Man kann (1) als „‚nichtlineares Gleichungssystem“ 
in den Komponenten von d auffassen und damit eine Reihe von Näherungsverfahren zur Be- 


' stimmung der charakteristischen Zahl und des Eigenvektors heranziehen, von denen hier aus- 
- schließlich das Newtonsche Verfahren?) behandelt werden soll. Geht man von einem bereits 


bekannten Näherungsvektor v, aus, dann erhält man nach Newton die Verbesserung 
I, —{AN, dam, ..., da} für v=1,23,3,... durch 


| ög\ | og 3% 
ur), nt mt mm) 


ED 


Die Determinante der Funktionalmatrix ©, ist im Falle einer einfachen charakteristischen 
Zahl A; nebst dazugehörigem Eigenvektor x; von Null verschieden. 


Ist kr die Abweichung der /-ten Komponente in b, vom richtigen Wert, dann gilt: 
BRITZ ENDE mit  MIE<SEN für I=12%..un 2.0. 8). 


(0 eine abschätzbare Konstante; im Falle eines linearen Gleichungssystems ist PM —= 0.) Für 
OK <1 liegt Konvergenz, und zwar quadratische Konvergenz vor. Eine Fehlerabschätzung 


kann ebenfalls angegeben werden. | 2% 
Beim Vorgehen mit konstanter Funktionalmatrix in (2) geht der quadratische Charakter 


_ der Konvergenz verloren; dafür wird der Rechnungsgang besonders einfach. 


2. Weitere Anwendungsmöglichkeiten des Newtonschen Verfahrens 
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aber auch insbesondere die nichtlineare Eigenwertaufgabe, z.B. MA) = AA? + B/A-+E mit 
dem Newtonschen Verfahren in gleicher Weise behandelt werden. 


3. Kompiexe charakteristische Zahlen bei der speziellen Eigenwertaufgabe und 
reeller Matrix 


Mit A=e-+ io und r—=y-+ iz kanneine Umschreibung auf ein 2 n dimensionales reelles 
Problem vorgenommen werden. 


2 3=—--004 
ee) a, 


Die Elimination von Real- oder Imaginärteil führt auf ein und dasselbe nichtlineare Problem: 

[(A—oE’+0]9y=0 bzw [A—-eC%+0]3=0 ........0). 
Das nichtlineare System (5) wird mit dem Newtonschen Verfahren behandelt und liefert Y 
oder 3. Der fehlende Imaginärteil 3 oder Realteil y kann dann nach (4) ermittelt werden. Für 
o=+0 sind 9 und 3 linear unabhängig. . 


Zur Praxis der Eigenwerteingrenzung 
Von N. J. Lehmann in Dresden 


Nach neueren Arbeiten ist es möglich, numerische Randwertberechnungen in der Form 
exakter Eingrenzungen durchzuführen !). Einige Hinweise sollen die a der Verfahren 
erleichtern. 

Es wird eine lineare, selbstadjungierte Rand wertaufgabe 


L(y)=ASy, Ss>0 U,(y)=0 BRENNEN (1) 
betrachtet, für die eine quadratisch integrable Gre e.nsche Funktion existieren soll. 


Um den Anschluß an frühere Arbeiten ?) zu gewinnen, wurde im Vortrag zunächst gezeigt, 
wie für (1) der Templesche Quotient 


Ei Ze 2 
7(y, oe (y) ST da— ef Liy)yda®) 
[ L(iy)yde— t [Sy?dx 
und die zugehörigen Einschließungssätze mit einer einfachen Spektralverschiebung (entspr. 


G1.(51) der in !) genannten Mittlg. I d. Verf.) aus dem ots 
werden. IRA y)ydz 


gewonnen 


Die weiteren Überlegungen sind auf die möglichst einfache Bestimmung grober Eigenwert- 


schranken gerichtet. Solche werden nämlich in den Verfahren !) zur Festlegung von Parameter- 
werten Zund damit zur Anwendung der Methode benötigt. Dazu sei angenommen, daß Gl. (1) 
höchstens endlich viele negative Eigenwerte besitzt. Dann gilt nach Courant mit V als 
Klasse der Vergleichsfunktionen (die Bildung von L(y) = quadr. integr. erlauben und die Rand- 
bedingungen erfüllen) 


{ L(y)yds 

Max Min (m — SLyyda _ .... 

v Mu\ =) [Sy da h Ay: ) 

Am Beispiel (Schornsteinschwingung) : 5 
yY+2[(1—-a)y=ıy, y0)=y()=y"’(d)=y”"()=0 ..:.. 06) 


sei gezeigt, wie man bisher über Gl. (2) mit den Vergleichssätzen der Variationsrechnung ?) zu 
(groben) Eigenwerteingrenzungen gelangte, und welche Schwierigkeiten dabei auftreten können. 


Zuerst wird der Quotient »(y) für yEV durch partielle Ines eon umgeformt und dann 
durch Rn enne Ne enponls eingegrenzt: 


| [wraa Iwrae—2] (1 a)(y')® da ] Wraa— 2] (Yde 
mi (y) = — Dm(y=I Eu ea Zar 


0 


probleme I u. II (Diss. 1948). Z. angew. Math. Mech. Bd, 29, 8. 341—356 u. Bd..30, 8 ‚1 16. 


*) Integrationen sind überall über das Grundgebiet von (1) zu erstrecken. 


2. Aufl. (1931), 496 S., insbes. S. 345 u. f£. 


°)R. Courantu.D. Hilbert: Methoden der mathematischen Ehyeik I. Berlin: Julius Springen. = 2 


‘ . 
:) H. Wielandt: Ein Tinschließungsentz für charakteristische Wurzeln a, Matrizen. Arch. 
Math. 1 (1949) 8.348—352. — N.J. Lehmann, Beiträge zur Ama EZ linearer = Eurawn 
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In dieser Gestalt sind bei m(y) und damit auch für die anderen Quotienten zur Bestimmung der 
Extremwerte (Eigenwerte) nur die geometrischen Bedingungen y(0) = y’(0)—= 0 notwendig!). 
Die Nebenaufgaben werden über die zugehörigen Differentialgleichungen, die man durch Vari- 
ation bestimmt, gelöst und liefern die Eingrenzungen 4 >A, >A?. Die Schwierigkeiten. be- 
ginnen mit der Berechnung der A}, der Rechenaufwand wırd zu groß. 

Diese Mängel lassen sich vermeiden, wenn die Nebenprobleme so gewählt werden können, 
daß ihre Eigenwerte durch eine bequem lösbare Hilfsaufgabe bestimmt sind. Das gelingt aber 
über einen neuen, sehr einfachen Hilfssatz (formuliert in bezug auf Gl. (2)). 

Ist r(x) eine für & > A, stetige und monoton wachsende Funktion von x (also etwa r’(x)> 0 
für 2>A,), so wird 


v 


Max Min IR) ra) A sen (5) 
yer 


Die Begründung ergibt das monotone Wachstum von r (m) bei m >A,. Wegen m(y) > 4, folgt 
daraus Min r(m(y)) = r (Min m(y)), entsprechendes für Maxima und schließlich (5). 
Die Anwendung wird mit dem Beispiel (3) gezeigt. Als bequem lösbares Hilfsproblem 
bietet sich an: ; 
1 
[(w")da 
Max Min m (= — — 2(=13,%); 2(=485); 23807); +” . . (6) 


Dez yEeV- [y de 
0 


(V nach (3), Lösung über Diff.-Gl.) Damit soll m(y) aus (4) nach unten abgeschätzt werden; 
yEV und übliche Integralungleichungen ergeben 


1 1 1 af 1 \ 1/2 
f (1 — a)(y)%de= — fta —a)yYydas f y? ee] R {1 — 2)y”—y’}? ie] 


1 reffı “11 1 1/8 1 y2a(fı 1/72 1 y2) 
<|r was] "ra - war] "+ [roras]"}< [fra] rorarl + |Fwrael”), 


0 


1 rı Trac 1/2 1 1 1/2 
[was s [{ was) f waa] + lv’ A)y@)|, vl] = Y y’ del < fr wel > 


9:04 1/2 
< af 02% ie] 
0 


v(A)| = = vlrwo as) da 
0 10 


% 
Sy da 
° 


% x l2 
<|f fa fwrerae| da 
0 0 0) 


und man erhält 


my) <miy) <my)—2Ymylır YE + ——)= rm) - - - - 8. 
| 3 my) 


Die Eigenwerte der Nebenaufgaben gibt (6) und Satz (5), dessen Voraussetzungeh mit 
r’(2<>12)>0 erfüllt sind, sie grenzen 4; ein: 


3 1 PR>1 222415380724 2 3581 
>> Zayaelı V: Ar *eheml 22 
| ie =) 485 > 4, > 402 


Mit diesem Ergebnis überblickt man auch den Einfluß des zweiten Gliedes der Aufgabe (3), 
entsprechend der Ordnung kommt es mit YA? zur Wirkung. Hierin liegt ein besonderer Vorteil 
dieser Abschätzungen, die auch für die asymptotische Eigenwertberechnung mit Vorteil be- 
nutzt werden können. ‘ 


(8). 


Zum Graeffe-Verfahren und Horner-Schema bei komplexen Wurzeln 
Von Rudolf Zurmühlin Darmstadt | 
‚Das Graeffe-Verfahren zur Auflösung algebraischer Gleichungen liefert im Falle konjugiert 


komplexer Wurzeln unmittelbar nur den Betrag. Die Bestimmung des Realteiles ist auf rein 


reellem Wege möglich unter Benutzung des reellen Horner-Schemas nach Collatz?) für 


1) R.Courantu.D.Hilbert:a.a. O. insbes. 8.179. — Vgl.auch E. Kamke: Über die.definiten 
 selbstadjungierten Eigenwertaufgaben bei gewöhnlichen linearen Differentialgleichungen. IV. Math, Z. Bd. 48 
(1942/43), S. 67—100, hier werden die geometr. Randbedingungen ‚wesentliche‘ genannt. 

- 2%) Z. angew. Math. Mech. 20 (1940), S. 235— 236. - 
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den 


komplexes Argument 2, —=&-iß. Dabei Division des Polynoms fie) d Qu 


ausdruck g(x) = 2 — px — q mit p= 2a, g= —r? gemäß ; 


She +ArtB fa)=An+B. 


D 


urch 


y v 2 hr x 


Pr 


Ist ©, eine Wurzel, so verschwinden A und B gleichzeitig. Ist nun gnach Graeffe Arachene 5 
genau bekannt, so lassensich Aund B auffassen als Funktionen der einen reellen Veränderlichen 9. & = 
Durch überschlägiges Berechnen von A und Bfür eine Anzahl runder p-Werte aus dem Intervall a: 
— 2r<p <-+2rverschafft man sich einen Überblick über den Verlauf der beiden Funktionen 
A(p) und B(p) und gewinnt durch Aufzeichnen einen Näherungswert für ihre gemeinsame Null- 
stelle 9, = 20. > * ee ER Se 

Eine Verbesserung nach Newton erster oder höherer Annäherung ist dadurch leicht möglich, 
daß die Ableitungen der Funktionen A(p), B(p)in einfacher Weise aus dem fortgesetztenHorner- 
Schema ähnlich wie beim gewöhnlichen Schema zu entnehmen sind, s. Gl. (4) bis (6). Damit 


2 


erhält man zur Berechnung der Korrektur 6p zwei Taylor-Entwicklungen, etwa bis zu Glieden 1 
2. Ordnung: ER 2 Sy er RU fe > “3 
A+Aöp+Anlöp?—=0, B+Böp+Bulp?=0......l. 

Ist q genau bekannt, so liefern beide Gleichungen gleiches öp (Kontrolle für g!). — Ist hingegen 
auch q nur angenähert bekannt!), so sind A und B aufzufassen als Funktionen zweier Verände- 
licher p und g. Auch dafür sind die Ableitungen leicht aus dem fortgesetzten Horner-Schma 
abzulesen, s. Gl. (4) bis (6). Die beiden Korrekturgleichungen 2. Grades lauten dann: rg i 
Ark ör Härte] + Anm? HaAnspöchan =). 4 
B+|B,65p + B,ög| + Bu (6öp)?+2 Ba» dp ög+ Ba (ög’—0} = 

Auflösung zweckmäßig als lineares Gleichungssystem nach dem umrahmten linearen Tilunter 
iterativer Berücksichtigung der quadratischen Glieder. Koeffizientendeterminane D=AB— 


A,B, = B}-+p B,A,— gq4% = positiv definite quadratische Form in A,, B, mit Diskrimi- FR 
nante 82 >0. Sie verschwindet genau für f’(&,) = f(x) = 0, also bei Mehrfachwurzen. 


_ ae az 


Vollständiges Horner-Schema: 


Ay An-ı .Im-2 .°*° 17; GG: dA BE 


“ Rs ar er , ’ A: 7 BE Se eye Be FUSS ns 34 
q Ian 1% 790,4 da Tag 
' ’ = ı ‚m , e st re Sr } s 


P An 


Taylor-Koeffizienten (Indizes p und g bezeichnen partielle Abl un, 


. 


A 2m. as N en De 
4, =A =) BET 
A, Et As Er 

1 DZ 

31 Ap = Aı = = 

1 2 DH 77 2 

91 Apg —=4, = a, % 

1 1 ER - BR 

57 Ag = Ay = % Er IR 


e2 


1, 7,B, bei Auflösung. auf instrumentellem Wege; v 
g en‘. 
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= 1 1 
> at re ie 2% 31 Appp = A = ga;" 
{ 1 DI 1 1 Pe 
sts Binz Br 31 Bypa = Bi = ga + 3% 
1 Er 1 ERTIERE . (6) 
3! Apga == An= 2 3] Boga = Pas = ga! | 3 a, 
e e u 1 11 
5 Ey. Ana, 51 Ba=Ba=g9% + a, 


Be „ei 0... 8 


Methoden der Theorie der Integraltransformationen 
Von L. Schwarz in Oberwolfach 


— 


s Wir geben im folgenden, um Einsicht in die funktionentheoretische und algebraische 
Struktur von Funktionaltransformationen zu erhalten, zwei Methoden an, die geeignet sind, die 
Lösung singulärer Integralgleichungen (1. Art) zu erleichtern. 

# 


I. Die singuläre Integralgleichung .- 
b 


g(y)= sW x) f(a)de 


mit &y,)=a(y—)!+ln|y—a|- Er x)+ M(y,2); L, M ganz, sei vorgelegt. 
Betrachtet man die Umlaufrelationen der Integralgleichung beim Umlauf von y in der kom- 
plexen Zahlenebene um « und b, so wird Det auf folgenden Lösungsansatz geführt: 


N) = JR dal) da 
mit Ro) =&- (1-9 Wa) Win 4 2a bin DRAUE EL IOLIUEE 
L,M ganz, W(a)=Yb— 2) («— a), Way)=—aey+ — > (a+b) (+y)— ab, 
&a bi 0 Wlognat |. 6-0) (y— 2) [Wa y) + Wow 


Er: a L ist lösender Kern einer leicht angebbaren Volterraschen Integralgleichung mit dem Kern L, 

R die u. U. leichter zu lösen ist als die Ausgangsintegralgleichung. Die Bestimmung von M bleibt 

. offen. Ganz allgemein empfiehlt sich die Methode der Umlaufrelationen für verzweigte Kerne 
und liefert für den lösenden Kern Struktureigenschaften. 


II. Die zweite Methode besteht in folgendem: Es sei & der Modul der zerfallenden Funk- 


Er tionen > F(y) Gi (2), Dy der durch 2 =, und y erzeugte Körper der Differentialoperatoren 


>in Y: entsprechend ®,. Gilt dann = 
a y:K(y,«) = M.: K(y, «)mod 3, 
"52 RE D,:K(y,«) = N„:K(y, z)mod 8, 
Bi wobei MN ES niesende Differentialoperatoren aus ®, sind und MN „ihre Adjungierten, so 
- wird durch die Zuordnung „— M, D,— N, ein Automorphismus von ® definiert, der Auto- 
a  morphismus des Kerns K(y,&) oder "der zugehörigen Fredholmschen Integraltransformation 
heißt. Bekannt ist der Automorphismus der Laplacetransformation (d.h. des Kerns exp (— y%)): 
 y—D,„ Dy,——x. Der Hintereinanderschaltung zweier Transformationen entspricht das 
= Produkt der Automorphismen, der inversen Transformation der inverse Automorphismus. 
Gestattet also eine Fredholmsche Integralgleichung einen Automorphismus, so bestimme 
zu dem inversen. Automorphismus eine zugehörige Transformation, was sich in Beispielen 
schwierig herausgestellt hat, wende sie auf die vorgelegte Integralgleichung an 
‚eine Integralgleichung von identischem Automorphismus (y„— x, D,— D,), d.h. 
OR deren a genau u ist. Es gilt nämlich der Ss atz: K(y,®) 
da 


on . B Z. angew. Math. Mech. - 
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Dabei ist A eine quadratische Matrix rationaler Funktionen, p(y) eine Zeile, y(x) eine Spalte 
gleicher Länge und, wie man sieht, adjungiert zu o(x). 
Die Auflösung der Integralgleichungen vom identischen Automorphismus ist ein allgemein 
noch ungelöstes Problem. Mann kann indessen hier wertvolle Aufschlüsse mit Hilfe der ersten 
Methode gewinnen. Auch für numerische Lösungen empfiehlt sich oft die Reduktion auf eine 
Inte gralgleichung von identischem Automorphismus. Näheres siehe an anderer Stelle. 


E. Mathematik 


Schmiegsungsverfahren der konformen Abbildung 
Von J. Heinhold in München 5 


Unter Schmiegungsfunktionen 9,(2,) verstehen wir Funktionen, die 1. im Innern eines im 
Einheitskreis der z-Ebene gelegenen schlichten Bereiches $,, der den Nullpunkt als inneren Punkt 
enthält, eindeutige analytische Funktionen von 2, sind und %, durch 2,4, = 9,(2,) auf einen 
Bereich ®,+, derselben Art in der 2,+,-Ebene abbilden, 2. eine in |2,+,| < 1 eindeutige analytische 
Umkehrfunktion 2,—= @!(2,+,) besitzen, welche dort den Bedingungen des Schwaärzschen 
Lemmas genügt, |9(2,41)|< 1, 951(0)= 0, und 3. die Bedingung 9,(0) > 0 erfüllt. 

Aus diesen Bedingungen folgt, daß die Abbildungen 34, = 9%,(%), v=0,1,2,..., aus 
®, eine Folge von Bereichen ®, erzeugen, deren Ränder einen- Minimalbestand R, vom Null- 
punkt besitzen, der mit »— oo monoton gegen 1 wächst. Hieraus ergibt sich dann, daß 
lim 9 (Pn-ı (- - - (Po(20)) - » -)) = f(2) = w eine in ®, analytische Funktion von 2, ist, die das 
no 
Innere von ®, umkehrbar eindeutig und konform auf das Innere des Einheitskreises der w-Ebene 


abbildet, so daß das positive Richtungselement im Nullpunkt fest bleibt, woraus die Eindeutigkeit 
der Lösung dieser Abbildungsaufgabe folgt. 

Bei der praktischen Anwendung dieses Schmiegungsverfahrens handelt es sich darum 
Schmiegungsfunktionen anzugeben, die A. so einfach sind, daß sie eine bequeme graphische oder 
numerische Durchführung der Abbildungen gestatten und mehrmals hintereinander angewandt 
werden können, B. schon nach wenigen Schritten eine brauchbare Annäherung an den Ein- 
heitskreis liefern, sei es unabhängig von der Art der vorliegenden Bereiche oder durch eine 
besonders günstige Anpassung der g, an die auftretenden Bereiche 8,. 

Koebet) verwendet beim Beweis des Riemannschen Abbildungssatzes als Schmiegungs- 
funktionen o, die Abbildungen einer zweiblättrigen Einheitskreisscheibe mit einem nicht im In- 
nern von ®, gelegenen Verzweigungspunkt a, auf die einfache Einheitskreisscheibe. Diese Schmie- 
gungsfunktionen haben die Eigenschaft A, jedoch nicht die Eigenschaft B. Ringleb?) benutzt 
daher neben den Koebeschen Schmiegungsfunktionen noch die Abbildungen von Kreissicheln 
auf den Einheitskreis. Damit die Bedingungen für Schmiegungsfunktionen erfüllt sind, hat man 
hier die Kreissichel als Teil der Einheitskreisscheibe so zu wählen, daß %, im Innern der Kreis- 
sichel gelegen ist und nach Abbildung der Sichel auf die Einheitskreisscheibe letztere 

- noch so auf sich abzubilden, daß die Zuordnung der Nullpunkte mit den zugehörigen positiven 
Richtungselementen hergestellt wird. 

Durch die Koebeschen Schmiegungsfunktionen wird z.B. der von a, nach e'®e a, gerad- 
linig aufgeschlitzte Einheitskreis (X,) auf eine Kreissichel abgebildet. Man erhält eine bessere 
Annäherung des abzubildenden Bereiches an den Einheitskreis, wenn man (K,) sofort auf den gan- 


zen Einheitskreis abbildet®). Das kann am einfachsten mit Hilfe der Abbildung w -, (z + =) 


und ihrer inversen, sowie einer Ähnlichkeitstransformation geschehen und ist graphisch oder | 
numerisch (unter Benutzung einer Tafel für den Hyperbelcosinus komplexen Argumentes) bequem a 
durchzuführen. Man hat dabei, um die Bedingungen für Schmierungsfunktionen zu erfüllen, 
den zum Einheitskreis orthogonalen Stachel mit der Spitze a, so zu wählen, daß er keine inneren 
Punkte von ®, enthält, da sonst der Bereich ®,+, über den Einheitskreis hinaustritt#). Durch 


& ERS Koebe: Abhandlungen zur Theorie der konformen Abbildung I. J. reine angew. Math. 145 (1918), 
. N Ri BET Numerische und graphische Verfahren der konformen Abbildung: "Habilitationsschrift. 
Heide rg 
®) J. Heinhold: Ein Schmiegungsverfahren der konformen Atbildung. Se Bayer. Äiad.. PS 
Wiss, Math.-nat. Kl. 1948, 8. 203— 222. 
OBER *) Auch in einem solchen Falle liefert die Einziehung eines "außerhalb By gelegenen Geraden 
H EN eine Annäherung der sämtlichen Randpunkte des Bereiches ®, an den Einheitskreis, man Ben jedoch keine 
Ei im Äußeren des Einheitskreises gelegene Stacheln einziehen, : De 


Ra 
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- Zwischenschaltung einfacher Bilder, z. B. linearer, lassen sich auch krummlinige, z.B. Kreis- 
bogenstacheln ‚‚einziehen‘‘. Weiterhin kann man unter Verwendung des Bildes w — 2” und ihrer 
inversen auch mehrere regelmäßig am Einheitskreis angeordnete Geradenstacheln gleicher oder 
abwechselnd gleicher Länge gleichzeitig einziehen. Stachel- und Kreissichelbilder zusammen 
liefern ein Schmiegungsverfahren, das sich den abzubildenden Bereichen gut anpassen läßt und 
gerade für Bereiche, die sehr von der Kreisform abweichen, z. B. mehrfach verästelte, schon nach 
wenigen Schritten bei passender Wahl der Stacheln und Kreissicheln eine gute Annäherung an 
den Einheitskreis ergibt. Dabei hat man, nachdem die größten Abweichungen vom Einheitskreis 
beseitigt sind, bei Verwendung von Stachelabbildungen mit wachsender Zahl der Schritte 
auch die Anzahl der gleichzeitig einzuziehenden Stacheln zu vermehren. 

Ist der Rand R, des Bereiches ®, dem Einheitskreis schon sehr nahe gekommen, so wird 
die weitere Durchführung des Schmiegungsverfahrens unrentabel. Es empfiehlt sich dann mit 
einem speziell für solche Bereiche geeigneten Verfahren, z.B. dem Verfahren von Theodorsent), 
die Annäherung an den Einheitskreis zu verbessern. 


Eine reelle Deutung komplexer Vektoren 
Von Frank Löbellin München 


Die auf Felix Klein zurückgehende Behandlungsweise der Geraden des hyperbolischen 
Raumes mittels binärer quadratischer Formen kann zu einem Kalkül ausgebaut werden, der mit 
der euklidischen Vektorrechnung völlig übereinstimmt, sofern imaginäre Vektoren zugelassen 


werden. > 
Eine ausführliche Darstellung erscheint in Bd. 52 dr Mathematischen Zeitschrift (1950), 
5. 758.ff. 


Da hiernach die reellen Bewegungen des nichteuklidischen Raumes durch homogene ortho- 
gonale Substitutionen mit komplexen Koeffizienten dargestellt werden können, eröffnet sich eine 
neue Möglichkeit, die Zusammensetzung der Vierpole zu behandeln. 


Die Entwicklung willkürlicher Funktionen nach den Eigenfunktionen 
. des Drehparabols 


Von Herbert Buchholz, Darmstadt 
Führt man in die Wh.D. Gl., die bekanntlich die Normalform hat, die neue abhängige 


Veränderliche Y = y- z!!?ein und setzt hinterher z= —in mit n> 0 und den vorderen Para- 
meter «—=t-A mit A2 0, so entsteht die Eigenwert-D.Gl. 
d ( je (Ü u2/4 ") 
se 3m Y a We N ER RFER TEnE 1 
Bien Tre L (1) 


mit der für u > —1 im Nullpunkt endlich bleibenden oder verschwindenden Lösung 


Min.un (—in) _ im, Hinz ( BR S:} ) 

EEE FEN I TIEREN . ey | en EN 
a en (12) 

Diese Funktion ist zwar für reelle Werte von « und n nicht selbst reellwertig, sie wird es aber 

nach Multiplikation mit exp (+ T »-(1-+u)|. Das Verhalten von (la) für A— ersieht 


man auch bei komplexen Werten am einfachsten aus der für alle Werte von x gleichmäßig 
konvergenten Entwicklung _ 


M „(—i RR 
ey Di PP). Tote Tram POW-I uw. 2) 
(—in)utae = era Pem)=--—-! 

” Hierin sind die Pl gerade Polynome in n, die sämtliche geräde Potenzen von n mit dem kleinsten 

#E Exponenten 2 Ei und dem höchsten Exponenten 2» enthalten. Sie sind allgemein definiert 

F 2 - . h 

5 durch das Umlaufsintegral = - 

FR 04) 00° 

a 17: Wi, (ein) + inj2- [&tge— (Sint ER di 


EN Th. end raen: Theory of wing sections of arbitrary shape. N. A. C. A. Report Nr. 411 (1931), 
 Theodorsen and EA: General potential theory of arbitrary wing sections. N. A. Ü. A. Report Nr. 452. 
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Es wurde bereits in dem Bericht des Velaseers in der ZAMM, Berl 1 23 (1943), darauf = 
hingewiesen, daß die zu verschiedenen Werten von A gehörigen Eigenfunktionen (la) z.B. im 
Grundgebiet 0...no ein vollständiges und orthogonales Funktionensystem bilden, dessen Eigen- - 
werte Ay zugleich die einfachen Nullstellen der Funktion M, » (— in.) in bezug auf den vor- 


Ey 
deren Parameter # bilden. Nach dem Sturm-Lionvilleschen Satz folgt daraus ohne witers 0 
die Möglichkeit, jede willkürliche, integrierbare Funktion in dem Bereich 0 <as Ei s no auf er 
Fouriersche Weise nach Eigenfunktionen von (1) zu entwickeln. 3 = = e 
Eine besondere Untersuchung erfordern jedoch die Fälle, wo entweder die re Ge 
a dieses Bereichs gleich Null wird oder die obere n, ins Unendliche rückt oder aber heides 
zugleich eintritt. Um die Möglichkeit und die Art der Darstellung einer willkürlichen ERnEUGES 
in diesen Grenzfällen nachzuweisen und herauszufinden, hat der Verfasser den von TILEeN=AE 
marsh [1] eingeschlagenen Weg beschritten. Die Ergebnisse sind in Kürze die folgenden: _ 
1. Liegt als Entwicklungsbereich das Intervall 0OZn <n vor, so gilt für In) die Datz ae 
stellung er 
ker: wer No x 
Cim. FM 


1 > er. ; i2, Fe . fe : 
2, MALHBUR LE oM u ino) TE fray Fr 
A = Be h 0) = mil e N 4) & 2 g 


„>—1) | En 
falls die in (4) vorkommendex, Integrale integrierbar sind. Di Ay sind darin die unendlich B. 
einfachen und reellen Nullstellen A, <A, <A,... von N, (in) hinsichtlich A. Die Ent- 


wicklungskoeffizienten von (4) stimmen dann überein mit dene? die sich aus dem formalen 


Ansatz SR Sn B 
= Mi,, (— ı n he Le e re ni, 2 : 
NOE Aa  EEZ EEE 
Be Mr: En 
nach einer Behandlung auf die Fouriersche Weise ergeben. a er ee 
Als Beispiel einer solchen Entwicklung werde angeführt: GN Es: Bee: 1 
1+u x BEN - 5 7 a ü ER 
S re: = 
b 2 2 


ae (N w>— BEE 

7 = 1 Er &, TE - ; 
ee Se: 

v & ar) IE : 
2. Liegt der Bereich 0 n <x vor, so N wie in der Theorie der Z,ylinderfunktionen, 

an Stelle der Reihenentwicklung eine Dar lunz durch ein .Dopgehnteßzeh en der G Sr 


ce e 

ultin) M Rz m EN 
f nf rer er rer. = vr san) in 
| ws ws: 


„ Diese Darstellung ist zuerst von A: Erdel yi cn len und fi 
“ leitung in der zweiten Auflage der Formelsammlung von Magnus "und 
Die Originalarbeit ist jedoch anscheinend noch nicht erschienen. 
> 3. Eine noch wieder. andere Entwicklung gilt im; Intervall nz <o@. 
ürzung ER 

; a 

. 5 4 

= vBic —i9)= 


nr 


4 al. 


ee 
W;. a v0)" Zr nn) 2 


—o 3 x = 


\M. 
Auch i in n (6) ist wie in (4) die Auswahl der. 
getroffen Woran — für NZ == "10 Verschwinden 
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Auflösung von Polynomgleichungen auf elektrischem Wege 
Von Wilhelm Bader in Stuttgart 


Bei Schwingungs- und Regelungsaufgaben, in der Elektrotechnik und im Maschinenbau 
überdies bei der Ermittlung von Anordnungen vorgeschriebener Frequenzabhängigkeit und in 
„vielen anderen Fällen sind Polynomgleichungen höheren Grades zu lösen. Im folgenden wird 
über zwei Verfahren berichtet, nach welchen man die Nullstellen auf elektrischem Wege finden 
kann. In beiden Fällen wird aus rein technischen Gründen die gegebene Gleichung n-ten Grades 
mit z multipliziert; so entsteht 


> n+1 L 
v=i 


Bei der Entwicklung des ersten Verfahrens glaubte ich, daß es keine physikalische Erscheinung 
. gäbe, bei welcher die Nullstellen als unmittelbar meßbare Größen erscheinen. So blieb nichts 
übrig, als versuchsweise 2 zu verändern, die Summe w auf elektrischem Wege zu bilden und jene 
Werte z festzustellen, für die sie verschwindet. z besitzt als komplexe Zahl zwei Freiheitsgrade 
und in der Summe müssen gleichzeitig zwei Komponenten, etwa Real- und Imaginärteil, ver- 
schwinden. Die Zerlegung von z selbst in Real- und Imaginärteil führt nicht zu Ausdrücken, 
. die sich elektrisch leicht auswerten lassen. Geeignet ist nur die Polarform 2= re‘? und es ist 
daher nicht verwunderlich, daß alle Vorschläge, Gleichungen durch versuchsweise Einstellung 
von 2 zu lösen, in ihren Grundzügen einander ähnlich sind, während sie natürlich in der Wahl 
der einzelnen Bausteine und damit im Aufwand sowie in der Schnelligkeit und der Genauigkeit 
der Auswertung voneinander sich unterscheiden. Wir schreiben statt (1) 
n+1 F 
N een een (2) 


er v=1 
oder auch 


W— 3 a,r’ cosyp ti a,r’sinvp. \ (2b) 
nn nn - nn nn FR Ä £ 


= U +1 ® 


Es gibt nun drei Möglichkeiten mit unterschiedlicher Schnelligkeit der Auswertung. Man kann r. 
und plangsam von Hand verändern und unter viel Zeitaufwand jene Wertepaare r und suchen, ' 
bei denen die elektrisch gemessene Summe w verschwindet (vgl. [1...4]). Man kann weiterhin 
den Betrag r langsam von Hand und den Phasenwinkel 9 selbsttätig sehr rasch ändern, so daß 
. „man scheinbar nur eine Größe, nämlich 7, abzusuchen hat. Es ist also = ot, wobei t die Zeit 
und w eine genügend hohe Winkelgeschwindigkeit bedeutet. Dann entsteht nach (2a) bei fest- 
 gehaltenem r eine geschlossene zykloidenartige Kurve. Jedesmal, wenn bei allmählicher Steige- 
"rung von r ein oder mehrere Äste der Kurve w durch den Nullpunkt gehen, ist der Betrag einer 
oder mehrer Nullstellen gewonnen. Man muß noch auf irgendeinem Wege die zugehörige Phase o 
ermitteln. Schließlich kann man selbsttätig den Betrag r rasch und irgendwie periodisch und o 
‚noch um eine Größenordnung rascher ändern, so daß z eine eng gewundene Spirale in der Gauß- 
‘ebene durchläuft, die man durch den Strahl eines Elektronenstrahloszillographen schreiben läßt. 
‚Der Strahl soll aber nur dann aufgehellt werden, so oft im Lauf der spiraligen Abtastung das 
Polynom w zu Null geworden ist. Dieses Verfahren wurde in verschiedenen Spielarten von den 
Herren Tischler, Rasch und Glubrecht vorgeschlagen; die Veröffentlichung [8] läßt nicht 
erkennen, ob ein Gerät fertiggestellt und welche Genauigkeit erzielt wurde. 
0. Der Verfasser hat sich für die zweite Möglichkeit entschieden, weilsie offenbar der Forderung 
nach hinreichend rascher und genauer Gleichungsauflösung am besten Rechnung trägt (vgl. [5, 
 6,7]). Hierbei benötigt man n+1 phasenstarre, durch Siebe genau sinusförmig geformte Wechsel- 
_  spannungen der Kreisfrequenzen ®, 2, ..., (n+1)o. Die Amplitude der Wechselspannung 
Nr.» e=1,..., n—+-1) wird zunächst von Hand mit Drehknopf auf den Koeffizienten a, ein- 
gestellt. In n+1 mit einander gekuppelten und sehr genauen, d.h. mechanischen Potenzier- 
ben werden die einzelnen Amplituden weiterhin mit dem Faktor r” versehen. Die so ge- 
wonnenen Wechselspannungen ‘werden in je einem Netzwerk in zwei um z/2 gegeneinander 
versetzt Teilspannungen a,r” cos vot und ayr’sin ywt aufgespalten. Nun werden gemäß (2b) 
it g t) die beiden Summen u und » elektrisch gebildet, die nach genügender Verstärkung 
nander senkrechten Ablenksystemen einer Elektronenstrahlröhre zugeführt wer- 
Schirm zur w-Ebene, auf der die Kurve w(re‘?) für langsam von Hand ver- 
ER zunehmenden u p aufleuchtet. Wenn die Kurve durch 
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Bild 1. Bestimmung einer Nullstelle bei einer Gleichung 
8. Grades 
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Bild 2. w=2-+2° für |2|=1; Lichtkante noch abseits 


den Betrag einer oder mehrerer Nullstellen. Die noch unbekannte Phase bestimmt man mit 
einem auf beliebigen Winkel @ einstellbaren induktiven Abnehmer, der einen Impuls und damit 
eine plötzlich ansteigende und allmählich abklingende Aufhellung der w-Kurve bewirkt, so oft 


Techn. Hochschule Sturfgeer 
Prof. Bader 


Biid3. Fertiges Gerät 


Instiit für Theorie ver Eiektrofsehmw 


ein mit der Wechselstrom- 
maschine 1 umlaufendes 
Eisenstück an ihm vorbei- 
streicht. Indem man den Ab- 
nehmer so lange verstellt, 
bis diese eindeutige Licht- 
kante entlang der w-Kurve 
in den Nullpunkt gewandert 
ist, erkennt man. an einer 
Winkeleinteilung den oder 
die zu r gehörigen Phasen- 
winkel der gesuchten Null- 
stelle. 


Bild 1 zeigt die Auflö- 
sung einer Gleichung 8. Gra- 
des. r ist richtig eingestellt 
und liefert den Betrag zweier 
konjugiert, komplexer Wur- , 
zeln, da zwei Äste durch den 
Nullpunkt laufen. Auch die 
Lichtkante des Phasenab- 


nehmers liegt schonim Nullpunkt. Für die bisher erzielte Genauigkeit diene folgendes Bei- 


spiel. Gleichung: 


— 0,211 + 0,524 2 — 0,96 22 + 0,84 23 — 24 + 0,576 25 — 0,746 28 + 0,407 2? + 0,571 2° 0. 


A Nullstellen: 
Messung 

T P 
1 0,5° 
2,04 180° 
0,752 70% 
0,752 — 76,5° 
0,945 121,5° 
0,938 —121,5° 
0,606 43° 
0,606 — 43° 


Rechnung 

# p 

1 0 
2,00 180° 
0,75 75°, 
0,75 — 75° 
0,95 120° 
0,95% — 120° 
0,6 45° 
0,6 — 45° 


Anlage gemacht worden. 
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Bild 2 läßt die zur Justierung benützte Kurve w= 2-+ 2° für r—= |z]—=1 erkennen. Bild 3 
zeigt das fertige Gerät, das nach den Plänen des Verfassers von Herrn K. Franke [7] ausgeführt 


wurde; es wurde auf der Tagung vorgeführt. 

Das zweite vom Verfasser vorgeschlagene Verfahren, welches mit anderen Lösungen offenbar 
keine Verwandtschaft besitzt und sich noch in Entwicklung befindet, nutzt die Tatsache, daß 
in (2b) Au=0 und Av=0 ist. Im kreisförmigen elektrolytischen Trog werden auf dessen 
Rande r=1 Realteil v und Imaginärfaktor v des Randwertes von w= u--i» durch passend 
eingeschaltete Randpotentiale erzwungen, .die zur Unterscheidung voneinander mit zwei ver- 


‚schiedenen Frequenzen schwingen. Damit sind « und v auch im Innern des Troges erzeugt. 


Führt man einen Taststift im Troge herum, so findet man die Nullstellen von w dort, wo die 
beiden akustisch oder optisch erkennbaren und durch ihre Frequenz voneinander unterscheid- 
baren Spannungen, also Real- und Imaginärteil, gleichzeitig verschwinden. Der Bereich |z| > 1 
wird bei beiden Verfahren durch einen Transformationsschalter in den Arbeitsbereich des Gerätes 
eingestülpt. 
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Automatisches lichtelektrisches Kurvenabtasten bei Integrieranlagen 
Von H.-J. Dreyer in Darmstadt 


Eine Integrieranlage !) liefert die Lösungen gewöhnlicher Differentialgleichungen 
ym =fl(y",...,y‘, y, x) im allgemeinen als gezeichnete Kurven. Die Veränderlichen x, y, y’, ... 
werden in der Anlage als physikalische Größen dargestellt, z. B. als Längen, Wellenumdrehungen 
oder elektrische Stromstärken. Es liegt nahe, auch eingehende veränderliche Daten, etwa ver- 
änderliche Koeffizienten der Differentialgleichung, mit Hilfe gezeichnet vorliegender Kurven 
zu verarbeiten. Dazu wird gewöhnlich das Kurvenblatt auf einen ‚‚Funktionstisch‘‘ gespannt, 
welcher es entsprechend der unabhängigen Veränderlichen in Abszissenrichtung bewegt; durch 
Drehen an einer Handkurbel hat man einen ordinatenparallel verschieblichen Fahrstift laufend 
auf der Kuvre zu halten. Diese Tätigkeit ist aber*fauf Dauer sehr mühsam und auch ungenau. 


Hazen, Jaeger und Brown?) haben deshalb schon zu der Integrieranlage von 
V. Bush), dem ‚Differeftial Analyzer‘‘, eine Einrichtung für automatisches Kurvenabtasten 
angegeben. Merkwürdigerweise sind aber diese Einrichtung und auch spätere an anderen Stellen 
entwickelte wieder verlassen worden, wahrscheinlich wegen technischer Mängel, Unzuverlässigkeit 
oder unzureichender Genauigkeit. Vielmehr hat sich bei Integrieranlagen die Gepflogenheit 
herausgebildet, das Abtasten von Kurven nach Möglichkeit zu vermeiden. Der dazu meist ein- 
geschlagene Weg, das ‚Erzeugen‘ von Funktionen als Lösungen von Sonderdifferentialglei- 
chungen, erfordert oft viele zusätzliche Rechengeräte und ist einer der Gründe für den großen 
Umfang des Bush -Differential Analyzers und seiner Nachbauten. 

Die Integrieranlage IPM-O tt !)in Darmstadt greift wieder auf das unmittelbare Abtasten 
von Funktionskurven zurück, und zwar mit einer eigens entwickelten automatischen, licht- 
elektrisch arbeitenden Kurvenabtasteinriehtung. Eine Photozelle tastet die eine Kante des ent- 
sprechend gezeichneten Kurvenstrichs mit einer Genauigkeit von + 0,1mm ab. Damit ist eine 
bei Geräten mit Zeichnungsgrundlage natürliche Genauigkeitsschranke erreicht. Die Abtast- 
fehler liegen weit unter 0,1%, des verfügbaren Ordinatenbereichs und auch unter derjenigen der 
übrigen Rechengeräte. Die Funktionstische sind so zu vollwertigen Rechengeräten der Integrier- 


\ 


1) Ygl. z.B. A. Walther u. H.-J. Dreyer: Die Integrieranlage IPM-Ott für gewöhnliche Diffe- 


rentialgleichungen. Naturwiss. 36 (1949), S. 199—206. 


L.Hazen, J. J.Jaegerand G.8.Brown: An automatic curve follower. Rev.sci. Instrum. 


7 (1936), 8. 353—357. 


3») V. Bush: The differential analyzer. A new machine for solving differential equations. J. Franklin 


Er 5 Be Inst. 212 (1931), 8. 447—488. 


FASEU 


F Mathematische Maschinen und Insert 


Die Zweckmäßigkeit des Kurvenabtastens ist offenkundig bei Differentialgleichungen mit 
empirisch gewonnenen Zusammenhängen zwischen den Veränderlichen und ihren Ableitungen. 
Beispielsweise bereitet es keine Schwierigkeiten, wenn in der Differentialgleichung die als Hyste- 
reseschleife gegebene Beziehung zwischen Strom und Spannung an einer eisenhaltigen Drossel- 
spule auftritt. Noch häufiger kann man das Kurvenabtasten für solche formelmäßigen Aus- 
drücke verwenden, die in verwickelter Weise von nur einer einzigen Veränderlichen oder von 
einer durch Rechengeräte bildbaren Kombination mehrerer Veränderlicher abhängen. So wurde 
in Darmstadt die Differentialgleichung!) 
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gelöst. Sie beschreibt thermodynamische und chemische Vorgänge in der Flammenzone einer 
Gasflamme; a, b, c, d, psind Konstante. Nur 2 Integriergeräte, 1 Multipliziergerät und 2 Addier- 
. geräte waren erforderlich, ferner 2 Funktionstische mit Kurvenabtastung: der eine zum Eingeben 
der Funktion f(T) = T3!2 e-XIT, der andere zum Potenzieren des durch Rechengeräte aus 7’ | 
und 7 gebildeten Klammerausdrucks. Das Ersetzen dieser beiden Funktionstisch oe durch 
„erzeugende‘ Rechenschaltungen würde weit mehr zusätzliche Rechengeräte erfordern, als sonst £ 
zur Lösung der Differentialgleichung überhaupt erforderlich sind, so daß diese Gleichung nur 
mit einer sehr großen Anlage angreifbar gewesen wäre. — Diese beiden Beispiele zeigen, daß 
durch automatisch arbeitende Funktionstische der Anwendungsbereich einer Integrieranlage auf N 
viele wichtige Differentialgleichungen erweitert wird, die sonst gar nicht oder mit großem geräte- NE 
mäßigem Aufwand gelöst werden können. 


Programmgesteuerte Rechenmaschinen in Deutschland 
j Von K. Zuse in Neukirchen 


Der erste Versuch, eine Rechenmaschine zu bauen, die eine Folge von Rechenoperationen 
nach vorgegebenen Programm durchführt, wurde von dem Engländer Ch. Babbage in der ersten 
Hälfte des 19. Jahrhunderts unternommen. Seine Arbeiten führten damals mangels Eealgucie 
technischer Mittel zu keinem Erfolg. x 


In neuester Zeit wurden vornehmlich in den USA, in England und Deutschland Entwick- 
lungsarbeiten geführt, deren Ergebnis die ‚„‚programmgesteuerten‘‘ Rechenmaschinen sind. Eine 
solche Maschine führt die arithmetischen Operationen in einem ‚‚Rechenwerk‘“ durch und kann 

: die für die Rechnung benötigten Angaben in einem ‚‚Speicherwerk“‘ speichern. Die Reihenfolge 
\ der entsprechend der jeweiligen Aufgabe vorzunehmenden Rechenoperationen, sowie der Trans- 

_ port von Angaben vom Speicherwerk ins Rechenwerk und umgekehrt wa durch ein BES Er 

werk‘‘ gesteuert. Se 


Die durch den Vortragenden i in Deutschland seit 1936 betriebene Rechengeräteentwicklung | 


} hat folgende Merkmale: TEEN RP Se. 
; „ 1. Die Lösung einer Aufgabe wird in einzelne RSskenopetationen aufgeteilt, deren Gesamt- FE 
heit den Rechenplan bzw. das Programm bildet. S £ Lo 


2. Die Rechenmaschine arbeitet im Dualzahlensystem, wofür das erste Rechenwerk 198 
fertiggestellt wurde. Fa 


3. Um eine selbsttätige Bestimmung der Lage des Kommas durch die Maschine zu mh >= 
lichen, werden Zahlenwerte in der Maschine in halblogarithmischer Form därgestell ine 
Dezimalzahl y hat in der Maschine die Form y = 2%x b, wobei a BARZEEh NE: Ei 
zwischen 1 und 2 ist. 


4. Sämtliche Rechenschaltungen werden mit Relais ausgeführt. Drei 
bisher zum Aulbau der Rechengeräte verwendetä DI BERZ N 
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wurden, 


c) Röhrenrelais (nach Schreyer), welche als ele kt 
geschwindigkeiten zulassen. Bis 


ı) E. Bartholome, H.-J. Dres und En 
in et einer Wärmeflamme durch Lösen eines Kippe 
. Elektrochem. 54 (1950), 8. 246252. 
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Nach verschiedenen Modellkonstruktionen vor dem Kriege wurde 1941 die erste programm- 
gesteuerte Rechenmaschine in Betrieb genommen. Weitere Maschinen wurden gebaut, die jedoch 
infolge des Bombenkrieges nicht mehr zum Einsatz kamen. Gegenwärtig ist ein nach dem Kriege . 
behelfsmäßig fertiggestelltes Gerät im Betrieb, mit dem explizit gegebene Aufgaben gelöst wer- 
den können, welche auf die Operationen +, —, X, :, Y zurückgeführt werden können. Zu den 
zunächst versuchsweise gerechneten Aufgaben gehören: 

Determinanten, Systeme von linearen Gleichungen, verschiedene Differentialgleichungen 
(einfache und gekoppelte Differentialgleichungen erster Ordnung), Iterative Verfahren für Glei- 
chungen höheren Grades, Reihenansätze für transzendente und trigonometrische Funktionen, 
Näherungsansätze für 3. Wurzel u. dgl., Matrizenrechnungen (Matrizenmultiplikation, Rezi- 
proke Matrix), Eigenwertprobleme, komplexe Polynome höheren Grades. 


Die Rechengeschwindigkeit beträgt ca. 1 Multiplikation pro Sekunde. 


Entsprechend den Aufgabenbereichen haben sich im Laufe der bisherigen Entwicklung 

2 Gruppen von Rechengeräten herausgestellt, von denen die eine vornehmlich zur Durchführung 

wissenschaftlicher Rechnungen bestimmt ist, und die andere in erster Linie zur Lösung technischer 
Aufgaben eingesetzt wird. 


Die außerhalb Deutschlands erbauten programmgesteuerten Rechengeräte dürften nach 
dieser Einteilung größtenteils zu den ‚‚wissenschaftlichen‘‘ Rechenmaschinen gerechnet werden. 
Ihre wichtigsten Bauelemente sind Elektronenröhren. Die erzielten hohen Rechengeschwindig- 
keiten bedingen einen beträchtlichen materiellen Aufwand, sowie eine sehr komplizierte Hand- 
habung der Geräte, insbesondere ihrer Programmwerke. Man verzichtet meist auf eine automa- 
tische Bestimmung der Komma-Lage. 2 

Die deutschen Geräte sind hingegen zu den ‚‚technischen‘‘ Rechenmaschinen zu zählen. 
Das bedeutet, daß auch verhältnismäßig kurze Aufgaben und solche, die nicht oft zu iterieren 
sind, mit Vorteil maschinell gerechnet werden können. Aus dieser Aufgabenstellung heraus 
wurde besonderer Wert auf leichte Bedienbarkeit gelegt. Der Aufbau der Schaltungen wurde so 
getroffen, daß die Maschine schnell von einer Aufgabe auf die Rechnung einer beliebigen anderen 
umgestellt werden kann. Durch die automatische Bestimmung der Komma-Lage erübrigt sich 
eine vorherige Überschlagsrechnung zur Ermittlung der Größenordnung der Resultate. 

Die zukünftige Entwicklung wird voraussichtlich eine weitere Spezialisierung programm- 

' gesteuerter Rechenmaschinen in den bereits oben erwähnten Anwendungsgebieten bringen, zu 
denen in abgewandelter Form noch die programmgesteuerten Rechengeräte für das kaufmännische 
Rechnungswesen hinzutreten werden. 


Eine Bemerkung zum Thema „Schachmaschine“ 
3 Von Th. Fromme in Weil/Rh. 


1938 hat Shannon bereits die Äquivalenz zwischen logischen Operationen und den 
Schaltoperationen von Relaiskreisen gezeigt. Seither weiß man, daß jedes im Aussagenkalkül 
formulierbare Problem durch einen ‚elektrischen Schaltmechanismus darzustellen ist, also auch 
das Schachspiel. 

Verschiedene Autoren haben bereits zu der Programmgestaltung einer Schachmaschine 
beigetragen. Hier soll indessen nur auf ein grundsätzliches Problem hingewiesen werden, das 
bei der Schachmaschine besonders klar zutage tritt. 

Wir betrachten eine Maschine, die streng determiniert arbeitet und solgends Regeln erfüllt: 

1. Jede Stellung des Schachspiels sei darstellbar. 

2. Die nach den Spielregeln des Schachs in jeder Stellung möglichen Züge seien, beispiels- 
weise lexikographisch, abzufragen. 

3. Jeder Zug wird zurückgenommen, wenn es eine Antwort für den Gegner gibt, die einen 

 Gewinnzug bedeutet. Damit ist der zurückgenommene Zug ein Sperrzug geworden. 
NEN. 4. Ein Gewinnzug ist a) das Matt, b) ein Zug auf den nur noch Sperrzüge möglich sind, 


die übrigen Züge sind unbewertet. 
er. 00.5, Die Maschine spielt unter allen nicht gesperrien Zügen den lexikographisch kleinsten. 
RN 2 Der Kürze halber lassen wir die entsprechenden Regeln für das Remis fort. 


Gr: _ Damit ist in jedem Moment für die Maschine eine eindeutige Arbeitsweise festgelegt. 
Trotzdem kann die Maschine bereits folgendes, ‚‚kleine‘‘ Problem nicht lösen: 
Re Welches sind die Sperrzüge, wenn die Maschine von der Ausgangsstellung bis zum 15. Zuge 
a (9) en durchrechnet ? 
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Bis zum 15. Zuge von Schwarz, das sind 30 Züge insgesamt, müssen wir in jeder Stellung 
mindestens 10 Varianten rechnen, so daß die Maschine 10%° Varianten durchspielen müßte. 
Setzen wir für die Dauer eines Stellungswechsels die Zeit eines Elementarprozesses von 10713 sec, 
dann braucht die Maschine zum Durchrechnen mehr als: 1017 sec, d.h. etwa 1010 Jahre. 

Nun kennt aber jeder Großmeister Züge, die von der Eingangsstellung aus bis zum 15. Zuge 
nicht zum Verlust führen. Worin ist die Überlegenheit des Schachspielers über diese (sehr spe- 
zielle) Schachmaschine begründet ? 

1. Der Schachmeister schafft sich neue, teils unscharfe — gefühlsmäßige Begriffe — wie 
Position, Hemmung, Materialübergewicht, d.h. Stellungsbewertungen. 

2. Er operiert mit diesen Begriffen in einer durch die Erfahrung gesteuerten, aber durchaus 
intuitiven Weise. 

Wir formulieren den Sachverhalt: Selbst bei einem präzisen Problem ist es möglich, daß 
wegen Zeit- oder Raummangel präzise Lösungsmethoden versagen, daß dagegen dasselbe Problem 
unter Zuhilfenahme von unpräzisen Denk- oder Rechenmitteln für die Praxis genügend gut zu 
lösen ist. Dieser Sachverhalt soll darauf hinweisen, daß 1. das Studium der Denkmethoden der 
Schachmeister eine wissenschaftliche Bedeutung, zum Beispiel im Hinblick auf ähnlich umfang- 
reiche Probleme der Physik haben kann, 


daß 2.in den Rechenautomaten die Einführung von Mittelwertprozessen,Zufallsverteilungen 
und ähnlichem von Wert sein kann, ja, unter Umständen erst die Lösung gewisser Probleme 
ermöglicht }). 


Die Funktionsrechenmaschine’) 
Von K. Ramsayer in. Stuttgart # 


Mit den gebräuchlichen numerischen Rechenmaschinen (Bürorechenmaschinen) können 
wir nur Rechenoperationen ausführen, die sich aus den vier Grundrechnungsarten zusammen- 
setzen. Zu allen übrigen Rechenoperationen müssen wir Funktionstafeln zu Hilfe nehmen, in 
welche wir mit gegebenen oder mit der Maschine berechneten Argumenten eingehen und die 
zugehörigen Funktionswerte entnehmen, welche dann entweder von Hand in die Rechenmaschine 
übertragen oder als Ergebnis aufgeschrieben werden. Dies ist nicht nur eine unbequeme und 
zeitraubende Belastung für den Rechner, sondern vor allem eine Fehlerquelle, welche die an sich 
große Sicherheit des Maschinenrechnens ganz wesentlich vermindert. Es ist deshalb wünschens- 
wert, für häufig vorkommende Rechnungen, die bisher mit Rechenmaschine und Tafel ausgeführt 
werden, eine Rechenmaschine zu schaffen, welche bei gleichbleibender Genauigkeit die Tafel- 
rechnung entbehrlich macht und trotzdem noch handlich bleibt. 

Diese Forderungen lassen sich erfüllen, wenn die Rechenmaschine mit einem Speicherwerk 
versehen wird, in welchem einzelne Werte der darzustellenden Funktion sowie Faktoren für die 
Interpolation von Zwischenwerten gespeichert sind, und wenn Übertragungseinrichtungen vor- 
gesehen werden, welche die gespeicherten Werte in die Rechenmaschine zur Berechnung der 
gesuchten Funktionswerte übertragen. Eine derartige „mechanische Funktionstafel‘“‘ wird be- 
sonders einfach, wenn man sich auf lineare Interpolation beschränkt. Nach der Patentanmeldung 
von.1941 kann man dann beispielsweise eine ‚‚mechanische Sinustafel“ folgendermaßen äus- 
bilden: Von der Funktion sin x werden die Funktionswerte für = 0°, 1°,2°, „.., 89° (Grund- 
werte) und die zugehörigen Tafeldifferenzen 7°: 7T,, durch Blechschablonen mit verschieden 
tiefen Einschnitten gespeichert. Soll z.B. sin 26°54’ » sin 26° + 54 T,, berechnet werden, 
so wird in der Rechenmaschine 26°54’ eingestellt. Durch das Einstellen der Gradzahl wird die 
Tafel, z. B. durch Drehen, „aufgeschlagen“. Dann werden durch besondere Übertragungsglieder, 
z.B. durch Zahnstangen, die Werte sin 26° und 7’,, aus den Blechschablonen abgegriffen und in 
das Resultatwerk (sin 26°) bzw. das Einstellwerk (7,,) der Rechenmaschine übertragen. An- 
schließend wird 7,,, mit der Anzahl der eingestellten Minuten, also mit 54 multipliziert. Dann 
steht im Resultatwerk der gesuchte Funktionswert, der für eine anschließende Addition, Sub- 
traktion oder Division sofort verwendbar ist und für eine Multiplikation mit Hilfe einer Rück- 
übertragungseinrichtung in das Einstellwerk der Rechenmaschine gebracht wird. i Rn 

Die mechanische Funktionstafe] läßt sich mit jeder beliebigen Rechenmaschinentype ver- 
binden, wenn die Übertragungsglieder der Maschinentype angepaßt werden. Das Volumen des 


') Es ist dem Verfasser erst später der Auisatz von L. Ridenourin „J. appl. Physics 1950, April, 
zur Kenntnis gekommen, in der das von Lehmer und Ulam bearbeitete „Monte-Carlo-System“ zur 
Durchrechnung kernphysikalischer Probleme erwähnt wird, ein System, das im oben angeführten Sinn un- 
präzise Mittel benutzt. 4 FE a 

°) Ein ausführlicher Aufsatz über Fraktionsrechenmaschinen wird in einem späteren Hefte dr 
Zeitschrift erscheinen. Es n } 2 ze 


Z. angew. Math, Mech. i : i h i 
Bd.30 Nr. 8/9 Aug./Sept. 1950 @. Stochastik 235 


Speicherwerks ist sehr klein. Es beträgt z.B. für eine mechanische Sinustafel mit.90 Grund- 
werten und 90 Tafeldifferenzen ohne Übertragungsglieder und Einstellvorrichtung etwa 0,5 Liter. 
Es macht deshalb keine Schwierigkeiten mehrere Funktionen in das Speicherwerk einzubauen. 
Auch läßt sich das Speicherwerk leicht auswecehselbar machen. 


Die Genauigkeit einer mechanischen Funktionstafel mit linearer Interpolation hängt ab 
von der Intervallbreite 3 zwischen den unmittelbar gespeicherten Funktionswerten und der 
zweiten Ableitung der Funktion. Bei der üblichen linearen Interpolation würde der maximale 
Interpolationsfehler im Intervall , <z sı,-+h|1/8-h2-f’"(x,)| betragen. Ersetzen wir 
die Kurve in den einzelnen Intervallen nicht durch ihre Sehnen, sondern durch vermittelnde 
Geraden, so beträgt der maximale Interpolationsfehler noch |1/16 A? f’(z,)| und tritt am 
Anfang, etwa in der Mitte und am Ende eines jeden Intervalls auf. Eine Funktion eignet sich 
also um so besser je kleiner f’”’(x) ist. Der maximale Interpolationsfehler einer mechanischen 
Sinus- oder Cosinustafel neuer Teilung mit einem Argumentabstand der Grundwerte von 19 
beträgt demnach 1,5 10°. Noch günstiger liegen die Verhältnisse bei der Funktion arc tg «. 
Der maximale Fehler beträgt hier für = 0,01 und |tga| < 1 nur- 0,8” bzw. 2,6. Wir 
können also hier durch Speicherung von nur 100 Grundwerten und 100 Tafeldifferenzen eine 
Genauigkeit von mindestens’ 0,6 - 10% erreichen. Dabei läßt sich die arc tg-Tafel für jeden 
beliebigen tg-Wert benutzen, wenn für |tgx|>1 der reziproke Wert verwendet wird. 


Die Genauigkeit läßt sich durch Erhöhen der Zahl der unmittelbar gespeicherten Funk- 
tionswerte oder durch Anwendung einer mehrstufigen Interpolation steigern. Die erste Möglich- 
keit ist allerdings bald erschöpft, da ja das Volumen der Maschine nicht zu groß werden darf, 
wenn sie noch handlich bleiben soll. Die zweite Möglichkeit erfordert etwas mehr technischen Auf- 
wand, dafür ist sie aber wesentlich raumsparender. Hierauf kann jedoch im Rahmen dieses 
Vortragsreferates nicht näher eingegangen werden. ; 


Der Verfasser hat eine Funktionsrechenmaschine mit linearer Interpolation entworfen, 
. welche von der Rechenmaschinenabteilung der Deutschen Telefonwerke und Kabelindustrie A.-G. 
(Hamann-Rechenmaschinen)in den Jahren 1941—1945 konstruktiv entwickelt wurde, aber leider 
nicht mehr fertiggestellt werden konnte. Diese Maschine sollte die Höhe und das Azimut eines 
Gestirns aus dem Zeitwinkel, der Deklination und der geographischen Breite nach dem sphä- 
rischen Cosinussatz und dem Sinus-Cosinussatz bzw. Sinussatz vollautomatisch berechnen. Der 
Vorteil. dieser Konstruktion gegenüber anderen Lösungen für dieselbe Aufgabe war die über- 
ragende Rechensicherheit und die sehr große Genauigkeit bei nur geringen Anforderungen an 
die Feinmechanik. 


Die Funktionsrechenmaschine kommt selbstverständlich nicht nur für Spezialrechnungen 
in-Betracht. Die Verbindung zwischen mechanischer Funktionstafel und Rechenmaschine kann 
auch so gewählt werden, daß die Maschine auch für sich allein in der bisher üblichen Weise ver- 
wendet werden kann und darüber hinaus noch die Möglichkeit gegeben ist, besonders häufig vor- 
kommende Rechnungen, zu denen bisher gedruckte Funktionstafeln erforderlich waren, ohne 
dieselben auszuführen. So können z. B. mit einer Eunktionsrechenmaschine, in welcher mecha-' 
nische Tafeln für die Funktionen sin x, cos & und arctg x eingebaut sind, die meisten in der ° 
niederen Geodäsie vorkommenden trigonometrischen Berechnungen ausreichend genau aus- 
geführt werden. . 


Wie jetzt bekannt wird, werden auch bei den neuen amerikanischen Rechenautomaten 
Funktionswerte nach dem Interpolationsprinzip gebildet. Zwischen der amerikanischen und der 
deutschen Entwicklung, die völlig unabhängig voneinander vor sich gingen, besteht jedoch ein 
grundlegender Unterschied. Bei den amerikanischen Maschinen handelt es sich um riesige Rechen- 
anlagen mit voluminösen Speicherwerken und völlig neuen Rechensystemen, bei der deutschen 
Funktionsrechenmaschine wurde eine Erweiterung des Rechenbereichs der gebräuchlichen Klein- 
rechenmaschinen mit einfachsten Mitteln angestrebt, so daß die Anschaffung einer Funktions- 

 rechenmaschine für jedes technische Büro erschwinglich sein wird. 


\ Ä G. Stochastik 


Die eaawendung stochastischer Methoden in der Biologie 
Von Wilhelm Ludwig in Heidelberg 


x Eine ‚theoretische Biologie“, die sich zur beobachtenden Biologie ähnlich verhält wie die 
ä theoretische zur Experimentalphysik, ist schon seit 3 Jahrzehnten im Werden und hat bereits 
e} bedeutsame nee gezeitigt. Eine Übersicht über ihre Anfänge und die Vielfalt ihrer 
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Probleme habe ich kürzlich an anderer Stelle gegeben!). Fast alle mathematischen Disziplinen, 
einschließlich z. B. Matrizenrechnung, Zahlentheorie und Logistik, wurden der Biologie dienstbar 
gemacht, und nur darum liegen seit Jahren biomathematische Fragen brach, weil es namentlich 
in Europa an der Zusammenarbeit biologisch interessierter Mathematiker mit Biologen mangelt. 


Stochastische Verfahren finden in der Biologie verschiedene Anwendung: Entweder (Ia) 
kann auf eine biologische Frage nur eine Wahrscheinlichkeitsaussage gegeben werden bzw. (Ib) 
es liegt überhaupt eine reine Wahrscheinlichkeitsfrage vor, oder es handelt sich (II) um ‚‚Stati- 
stische Methodik“, d.h. um Auswertung von Massendaten mit der Grundfrage: Kann das Resultat 
noch Zufall sein? 


Zu Ia sei nur erinnert an die Kaskadenprozesse, auch Auslösungsfolgen (v. Schelling) u.ä._ 
genannt, die sich um das Grundproblem gruppieren, mit welcher Wahrscheinlichkeit ein Indivi- - 
duum önach » Generationen genau oder mindestens x Nachfahren besitzt, wenn die Wahrschein- 
lichkeiten, daß «im Mittel a; Kinder produziert, a; = 1 betragen; über die a; können beliebige 
Annahmen gemacht werden. R. A. Fisher hat 1930 bei einem Versuche der Neubegründung 
von Darwins Zuchtwahllehre dieses Problem zuerst aufgegriffen. Auch zu Ib möge ein Bei- 
spiel] genügen: Wächst auf einer Wiese, die man sich quadratisch gerastert denkt, die Pflanzen- 
art Xin Nestern ünd seien sowohl Nestzahl/Quadrat wie andererseits die Nestumfänge Poisson- 
verteilt, so läßt sich eine ‚‚doppelte Poisson-Verteilung‘‘ errechnen, die man in der Natur verifi- 
ziert vorfindet. — Die eigentliche statistische Methodik (II) folgt in der Biologie meist dem 
Schema, daß man zu einer Frage (1), die nur an Massenmaterial gelöst werden kann, siebungs- 
freie Daten sammelt (2), dann die Eingangsfrage in Ansehung dieses Materials vermittels einer 
‚„Nullhypothese“ HZ genau präzisiert (3), jetzt die notwendigen statistischen Größen errech- 
net (4) und schließlich (5) aus diesen Größen (%, 02, x?,t, F usw.) jene Wahrscheinlichkeit P be- 
rechnet, mit der die Unterschiede zwischen den beobachteten und den durch H geforderten 
Werten (oder noch größere Unterschiede) zufallsmäßig zu erwarten wären. Da nun in der Biologie, 
etwa im Gegensatz zur Industrie, das Material meist mühsam gewonnen werden muß, sind alle 
ihre statistischen Verfahren dahin ausgerichtet, auch bei kleinstem Probenumfang noch gültig 
zu sein und vor allem sparende Versuchsumfänge zu planen. Die ‚‚Großzahlforschung‘‘ scheidet 
daher praktisch aus; auf Abweichung von Normalverteilung wird durch g, und 9, geprüft. Der 
Mittelwertsvergleich wird immer mehr durch die Varianzanalyse abgelöst, die es gestattet, in - 
einem komplexen Versuch mehrere beeinflussende Faktoren und ihre Wechselwirkungen zu 
prüfen. In der Züchtungskunde ist hier fast eine Spezialwissenschaft entstanden. Der Varia- 
bilitäten-(o?)-Vergleich ist-heute auch für Proben wechselnden und kleinsten Umfangs (n >1) 
gelöst. Für das X2-Verfahren wurden immer neue Anwendungen geschaffen: seine Zerlegung, 
analog der vom o? bei der Varianzanalyse; zur Schätzung von Parametern; für den Vergleich 
von Häufigkeiten aus Proben auch kleinsten Umfangs (n < 5); für Order-Teste (z.B. obin 
einer Geburtenfolge $ und 2 zufallsmäßig wechseln); und schließlich (aber theoretisch noch 
keineswegs beiriedigend geBärz) zur Zusammenfassung mehrerer unabhängiger P, je durch Um- 
rechnung in ein 2-F@-x°, in ein Gesamt-P. Die Anwendung der Variationsbreite als Charak- 


‚ teristikum hat sich bisher nur für kleinste Proben als vorteilhaft erwiesen, wobei stets als Norm 


die von 1 Individuum gelieferte Information <= 1/no? dient. Die negative Binomialverteilung 
wird zur Prüfung auf Wahrscheinlichkeitsansteckung verwendet sowie neuerdings zur Planung 


-minimalster Versuchsumfänge. Zur Schätzmethode des ‚„‚währscheinlichsten Werts“ (Fishers 


m. of max. likelihood) wurden von diesem jüngst sehr leistungsfähige numerische Lösungsmetho- 
den entwickelt... Noch unbefriedigend gelöst ist das Problem der Gruppenunterschiede. Eine 
Einzelfrage Jautet etwa: Sind alle Individuen von y verschiedenen Gruppen auf je m (vermutlich 


untereinander korrelierte Maße) gemessen, welcher Gruppe ist das Neuindividuum X, dsen 
Maße bekannt sind, am ehesten einzuordnen ? Weder die mathematisch exakten amerikanischen 
Methoden (Grundfaktorenanalyse) noch die von Fisher stammenden Discriminantfunktionen” BT 
können praktisch ohne Besitz ganz moderner Rechenmaschinen angewendet werden. BB SEE 


P} - 


Aeeidungen der mathematischen Statistik auf medizinische Statistik 


Von Karl Freudenberg in Berlin —_ er BEE 
Die medizinische Statistik ist an sich schon ein ausgesprochenes Grenzgebiet. a no FAR 
die Anwendung der mathematischen Statistik hinzu, dann wird der Grenzgebietscharakter noch WITH, 


ausgeprägter. Man muß deshalb die Gefahr im Auge behalten, daß die mathematische Seite 
überspitzt wird, während sachlogische Forderungen vernachlässigt werden. Die Kritik des 
materials zeigt ohnehin meist, daß dieses zu ungenau ist, als ae a ee .. zu 
stochastischer Methoden gerechtiertigt wäre. | x 


!) Verh. D. Zool. Ges. (Kiel), 1948. 
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Soweit die mathematische Statistik zu Berechnungen über die Gesamtsterblichkeit ange- 
wendet wird, kann sie auch für die einzelnen Todesursachen herangezogen werden, und dann 
kann man für diese gegebenenfalls noch abhängige und unabhängige Sterbewahrscheinlichkeiten 
unterscheiden. 

Die Variationsstatistik, die Quetelet bei der Verteilung von Körpermaßen einführte, hat 
ihre Hauptdomäne noch immer auf diesem Gebiete. Die Körpermaße werden für alle einzelnen 
Altersklassen auf die Form der Verteilungskurven untersucht, die zusätzliche Benützung der 
3. und 4. Momente brachte weitere Aufschlüsse. Dieses Gebiet wurde auch auf körperliche und 
geistige Leistungen erweitert. Um die Beziehungen zwischen Körpergröße und -gewicht zu 
analysieren, wurde die Methode der Streuungszerlegung benützt. 

Die von Lexis gefundene Gesetzmäßigkeit, daß die natürliche Sterblichkeit durch die 
Verteilung der Sterbealter gemäß der Normalkurve charakterisiert ist, konnte unterbaut 
werden, indem die Todesursachenstatistik zur Abgrenzung der natürlichen Todesfälle heran- 
gezogen wurde. 

Ebenfalls von Lexis stammt die Anwendung seiner Dispersionstheorie auf das Geschlechts- 
verhältnis der Geborenen. Spätere Untersuchungen auf diesem Gebiete, prüften die Konstanz 
des Geschlechtsverhältnisses bei Aufspaltung der Zahlen nach dem Alter der Mutter und der 
Geburtenfolge. Weitere Untersuchungen befaßten sich mit dem Geschlechtsverhältnis bei Mehr- 
lingen, wobei Weinberg durch Anwendung der Wahrscheinlichkeitsrechnung der Rückschluß 
auf den genetischen Zusammenhang gelang. 

Die Korrelationsrechnung hat ihr Hauptanwendungsgebiet innerhalb der Medizin in der 
Erfolgsstatistik. Meistens handelt es sich um eine einfache alternative Fragestellung nach Erfolg 
oder Nichterfolg bei Anwendung oder Nichtanwendung eines bestimmten Heilverfahrens. Es 
können aber auch quantitative Zusammenhänge so erfaßt werden. Die Bedeutung mehrerer 
Faktoren für den gleichen Erfolg wird untersucht, indem man partielle Korrelationskoeffizienten 
berechnet. 

In manchen Fällen wird die Mortalität allein nicht als ausreichende Darstellung anerkannt, 
dann muß man Morbidität und Letalität getrennt berechnen. Hierbei ergeben sich weitere Korre- 


lationen. Doch muß man auch Kritik anwenden, z.B. an den Zahlen der Todesursachen. Diese 


wird: unterstützt, wenn die Korrelationsrechnung z.B. zeigt, daß nach Standardisierung der 
Zahlen in vielen Ländern die Sterbezilfern für Neubildungen und für Altersschwäche eine straffe 
negative Korrelation aufweisen, so daß offenbar ein ziemlich konstanter Teil der Todesfälle an 
angeblicher Altersschwäche in Wirklichkeit Neubildungen zuzuschreiben ist. 

Der Zusammenhang zwischen körperlicher und geistiger Entwicklung von Bndert wurde 
mit Korrelationskoeffizienten untersucht. Ob ein solches Merkmal von einem nicht-quantitativen, 
wie z. B. dem Beruf des Vaters abhängig ist, konnte mittels des Pearsonschen Korrelations- 
verhältnisses gepiüit werden. 


Ausgewählte Probleme der Wirtschafts- und soziologischen Statistik 
Von Hans Kellerer in München 


Eine Folge der bisher zu geringen Pilege der mathematischen Statistik in Deutschland ist, 
daß wir an den großen Fortschritten des Auslandes auf diesem Gebiete nur ungenügend teil- 
nehmen können. Nur wenn hier ein entschiedener Wandel durch entsprechende Ausgestaltung 
des statistischen Unterrichtswesens auf den Hochschulen, durch Sondeıkurse und durch ein 
eigenes Schrilttum zu diesem Gegenstand geschailen wird, können viele Pıobleme der amtlichen 


‚und nichtamtlichen Statistik in Deutschland einwandireier, schneller und billiger gelöst weıden. 


Nachstehend «einige Hinweise auf Beispiele, bei denen mathematische Veriahrensweisen ange- 
bracht sind: 

Wie etwa die Einkommens- und Vermögensverteilung zeigt, spielt im Wirtschafts- und 
Sozialleben die Normalverteilung nur eine geringe Rolle. Vor allem die italienischen Statistiker 
unter Führung von C. Gini haben daher besondere „Konzentrationsmaße‘ auf- 


- gestellt, um die hier vorliegende ausgesprochene Disparität und die damit zusammenhängende 


Konzentration z.B. des Gesamteinkommens zu messen. Einen recht anschaulichen Überblick 


über die Konzentration bietet auch die „Konzentrationskurve“: Ist z.B. y=/(x) die Ein- 


DR „ kommensverteilung aller n Einkommensteuerpflichtigen, so ist die Parameterdarstellung der 


zugehörigen er nonaien 


e =— „lu 1 x; = — = [#0 de; m — Gesamteinkommen). 
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Die Fläche, die von der Konzentrationskurve und der ‚‚Gleichverteilungsgeraden““. (Diagonale) 


umschlossen wird, gibt die Hälfte des Ginischen Konzentrationsmaßes!). Mit dem genannten 
Problemkreis hängt die Frage der Rückführung asymmetrischer Verteilungen auf eine 
Normalverteilung mit Hilfe einer Variablentransformation, etwa der Form 


z=aln(e —x,)+ b . 
zusammen. 
In der Indexrechnung überwiegt die sachlogische Seite. Wie nützlich aber auch 


auf diesem Gebiet die kritische Mitarbeit eines mathematischen Statistikers ist, zeigen drei 
neuere Arbeiten von O0. Anderson?). 


Die Korrelationsrechnung wird in der amtlichen Statistik noch zu wenig benutzt. Ein 


schönes Anwendungsbeispiel, das eine vertiefte Untersuchung verdient, sind die Wirtschafts- 
rechnungen. Hierbei führen eine Anzahl von Arbeiter-, Angestellten- und Beamtenhaus- 
haltungen Aufzeichnungen über die Einnahmen und Ausgaben. Die Ausgabengruppen 2, (z.B. 
für Nahrungsmittel, Miete, Heizung und Beleuchtung, Bekleidung usw.) hängen in ihrer Größe 
vor allem ab von dem Einkommen z und der Anzahl y der Haushaltmitglieder. Es gilt annähernd 


y»—=ytb'2+%:Y. 


Dabei ist es zweckmäßig, den Begriff der Haushaltsperson zu verfeinern, indem man mit sog. 
‚‚Konsumeinheiten“ (Vollpersonen) arbeitet ?). 


Bekannt ist das Verfahren, wirtschaftliche Wachstumserscheinungen durch eine logistische 
Kurve zu beschreiben. Aber auch Begriffe wie Sterbetafeln (Abgangsordnungen), mittlere 
Lebensdauer, weitere Lebenserwartung und die dazugehörigen Berechnungsmethoden können 
auf Vorgänge des Wirtschafts- und Soziallebens übertragen werden *). Als praktisch durch- 
geführte Beispiele seien genannt: Abgangsordnungen für Automobile, sonstige Fahrzeuge, Eisen- 
bahnschwellen, Telegraphenmasten, verschiedene technische Güter, Sterbetafeln für Einzel- 
handelsgeschäfte und Aktiengesellschaften, Abgangsordnungen für die Betriebszugehörigkeit von 


Arbeitern und Angestellten in großen Unternehmungen, die Aufgliederung der Mietparteien nach 


ihrer Mietdauer beim Auszug, die durchschnittliche Lagerdauer in Handelsbetrieben. 


Vielfach ist der Statistiker gezwungen, sich mit Teilerhebungen oder Teilaufbereitungen 
zu begnügen. In der deutschen Statistik überwiegt noch stark die bewußte Auswahl; hierbei 
stützt man sich auf vorhandene Kenntnisse, z. B. auf Zahlen der Volkszählungen. Gegenüber 
diesem -Vorgehen hat die Stichprobenmethode den Vorteil, daß die Wahrscheinlich- 
keitsrechnung angewandt und damit die Güte der Ergebnisse gemessen werden kann. Die Theorie 
und Technik des Stichprobenverfahrens ist in den letzten Jahrzehnten vor allem in den angel- 
sächsischen Ländern zu einer großen Vollkommenheit entwickelt worden. Davon geben z.B. die 
verschiedenen Auswahlverfahren: geschichtete, systematische, Flächen-, Klumpenauswahl, mehr- 
stufiges, mehrphasiges Auswahlverfahren, fortschreitendes Stichprobenverfahren (sequential 
analysis) Zeugnis ®). Eine der Hauptaufgaben besteht darin, die Auswahlmethode zu finden, die 
unter Berücksichtigung der. Kosten möglichst wirksam ist, d.h. den geringsten Fehler (mittlere 
quadratische Abweichung) aller möglichen Stichprobenschätzwerte liefert. Eine ebenso wichtige 
wie reizvolle Aufgabe des mathematischen Statistikers liegt darin, seinen Teil zu einer vermehrten 

„und richtigen Anwendung dieser Verfahren in der amtlichen und nichtamtlichen Statistik bei- 
zutragen. 


!) Siehe in der deutschsprachigen Literatur z.B. E.P. Billeter: Über die Messung der Einkommens- 
Konzentration. Bern 1949. 1348. H.Kellerer: Die Konzentrationskurven in der angewandten Statistik. 
Archiv für mathematische Wirtschafts- und Sozialforschung. Band 3, Heft 1, S. 57—67. Leipzig 1937. 

°) OÖ. Anderson: &) Über die Neuberechnung von Indexziffern der Lebenshaltungskosten in 
Deutschland. Weltwirtschaftliches Archiv, Band 62, Heft 2, S. 175--197. Hamburg 1949. * 

b) Mehr Vorsicht mit Indexzahlen. Allgemeines Statistisches Archiv, 33. Band, 4. Heft. $. 472—479. 
München 1949. > 


c) Und dennoch mehr Vorsicht mit Indexzahlen. Allgemeines Statistisches Archiv, 34. Band, 1. Heft, 
S. 37—47. München 1950. ’ 


*) Siehe z. B.: Haushaltrechnungen von Familien unselbständig Erwerbender 1936/37 und 1937/38. 


Sonderheft 1942 der, „Volkswirtschaft“, herausgegeben vom Eidgenössischen volkswirtschaftlichen Departe- 
ment Bern 1942. | 


*) Diesen Nachweis führt die noch nicht veröffentlichte Arbeit: H. Kellerer, Übertragung be. 


völkerungsstatistischer Begriffe und Methoden auf das Wirtschafts- und Sozialleben. 


°) Siehe z. B.: F. Yates, Sampling Methods for Censuses and Surveys. Griffin & Co., London 199, 
3188. — H. Kellerer, Elementare Ausführungen zur Theorie und Technik des Stichprobenverfahrens. 
Mitteilungsblatt für mathematische Statistik Jahrgang 1, S. 96—114; 203—218; Jahrgang 2,8.36-—49, 


München 1949 und 1950. 
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Mathematisch-statistische Fragen der Fabrikationskontrolle 
Von Maria-Pia Geppert in Bad Nauheim 


In Deutschland wurde die statistische Fabrikationskontrolle (st. F.K.) seit ca. 1922 von 
DaevesundBeckelin Form der sog. ‚‚Großzahlforschung‘‘ angebahnt. Die entsprechende 
Entwicklung im Ausland fand ihren ersten abschließenden Ausdruck inShewharts Standard- 
werk 1931, in welchem die systematische Anwendung der klassischen Zufallstestmethodik auf 
die F. K. ausgebaut wurde. Erheblichen Aufschwung und neue Ausrichtung erfuhr die st. F,-K. 
im zweiten Weltkrieg ebenfalls in England und USA. (Walds „Sequenz-Analysis‘ u.a.). 

Während die rein wirtschaftliche Betrachtung eines Betriebes als Einzelzelle größerer 
Wirtschaftseinheiten sich im wesentlichen mit den Methoden der beschreibenden Statistik be- 
gnügen kann, bedarf die eigentliche F.-K., d.h. die statistische Überwachung der Produktion, 
des Material-, Maschinen-, Energieverbrauchs, der Arbeitsleistung usw., sei es zu Forschungs- 
zwecken mit dem Ziele quantitativer oder qualitativer Verbesserung oder Verbilligung des Ar- 
tikels, sei es im Hinblick auf seine kommerzielle Verwertung, in stärkerem Maße stochastisch 
unterbauter statistischer Methoden. 

Im Gegensatz zur Verwendung statistischer Methoden für rein wissenschaftliche Forschungs- 
zwecke (etwain der Biologie, Psychologie) ist für die st. F.-K. der Gesichtspunkt der Rentabilität 
des Autwandes an Kosten, Afbeitszeit, fachlicher Vorbildung usw. charakteristisch. In der 
laufenden Produktionskontrolle interessiert zwar der Einfluß der Zeit, aber nicht in konsta- 
tierendem Sinn wie in der Konjunkturforschung, sondern in kritischerem Sinne: überzufällige 
Veränderung einer beobachteten Größe oder eines aus solchen abgeleiteten Verteilungspara- 
meters deutet auf neu hinzutretende Ursachen hin, die entweder als ungünstig ausgemerät, oder 
als günstig wirkend systematisch erforscht und» verwertet werden ; im Hinblick auf die Benutzung 
des Artikels im Zusammenhang mit anderen (z.B. Passungen) wird oft auf geringste Varia- 
bilität, maximale Homogenität, besonderer Wert gelegt. Die im Hinblick auf den Kunden 
erfolgenden Prüfungen von Warenposten anhand kleiner Stichproben sind wahrscheinlichkeits- 
theoretisch besonders reizvoll. 


Grundsätzlich finden alle heute üblichen klassischen Zufallskriterien und Schätzmethoden 
in der F.-K. Verwendung, z. B. Beurteilung und Vergleich von Häufigkeiten (Ausschußhäufigkeit), 
von Häufigkeitsverteilungen (x?-Test) und von Mittelwerten (i-Test, Varianzanalyse), Unter- 
suchung von Zusammenhängen (Regression, einfache und partielle Korrelation) usw. In seinem 
Standard-Werk nahm Shewhart 1931 Ausbau und Anpassung der klassischen Testmethodik 
— soweit sie damals existierte — für die laufende F.-K. vor in Form von handlichen Kontroll- 
Diagrammen (‚‚charts‘“). 


Daeves’,,Großzahlforschung‘“ ist, von der Annahme ausgehend, daß eine Beobachtungs- 
größe entweder selber, oder mindestens nach geeigneter, zumeist logarithmischer, Transformation 
annähernd normal verteilt sei, dem speziellen Problem der Prüfung auf Homogenität gewidmet. 
Erweist sich eine Stichprobe der (nötigenfalls transformierten) Variablen bei Eintragung in 
Wahrscheinlichkeitsnetz als nicht normal verteilt, so wird sie nach einem von Beckel ange- 
gebenen zeichnerischen Verfahren in einige annähernd normal verteilte Komponenten zerlegt 
und durch die entsprechenden Mittelwerte und Streuungen charakterisiert. Für Variablen, die 
der genannten Annahme wenigstens approximativ genügen, und bei hinreichend großem Stich- 
probenumfang, besitzt das Verfahren gegenüber den exakteren, zufallskritischen, mathematisch 
eindeutigen, aber rechnerisch anspruchsvolleren Methoden der klassischen Theorie heuristischen 
Wert. : 

Mehr noch als bei der laufenden F.-K. treten bei der F.-K. im Verkehr mit dem Kunden 
Theorie und Planung kleiner Stichproben in den Vordergrund. Preisvereinbarung und Abnahme- 
bedingung lauten etwa: ein Warenposten Py von N Stücken wird mit einem bestimmten Preis 
bezahlt, wenn in ihm die Defekthäufigkeit unterhalb einer bestimmten Grenze liegt; ob dies 
der Fall ist, wird aus der Defekthäufigkeit in einer P, zu entnehmenden n-gliedrigen Stichprobe 
© beurteilt. Die Analogie zur Prüfung einer Hypothese in der Neyman-Pearsonschen Theorie 
führt zur Identifizierung des Produzenten- bzw. Käuferrisikos mit den Fehlern erster bzw. 
zweiter Art. 

| Im zweiten Weltkrieg entwickelte W ald aus der klassischen Testtheorie die, letztere als 
Spezialfall umfassende, Sequenz-Analysis. Während die klassischen Kriterien mit festen Stich- 
probenumfängen operieren, wird bei den Sequenztesten die Stichprobe schrittweise vergrößert, 
bis eine Entscheidung erzwungen ist. Da die Sequenzteste effizienter sind als die klassischen 
Kriterien, bedürfen sie zur Erreichung gleich zuverlässiger Entscheidungen im Durchschnitt 
wesentlich kleinerer Stichproben. : 
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Zur Vereinfachung der laufenden F.-K. wurde von Desmond, Stevensu.a.de 
Eich-Methode (gauging) entwickelt, bei welcher die genaue Messung einer Variablen durch .Ver- Ba 
gleich mit einem oder zwei festen Maßstäben ersetzt wird und die Besetzung der beiden extremen $ 
Klassen als Indikator für Veränderungen des Mittelwertes oder der Streuung dient. Zur Prü- we 
fung von Zentralwerten in kleinen Stichproben wurden von W als h neue, höchst einfache Kr 
terien entworfen, die im Gegensatz zu den herkömmlichen Methoden keinerlei Voraussetzungen 
über Normalverteilung im Kollektiv erfordern. Wo Be 
Erwägenswert scheint eine Untersuchung, wie weit sich auf die F.K, auch spezielle Me- ag 
thoden der in besonderem Maße entwickelten Biostatistik mit Nutzen übertragen lassen, etwa 
die zur Untersuchung von Dosis und Wirkung wichtige Probit-Analysis, oder die Grundfaktor- # 


/ analysis (der Psychologie), oder die noch ausbaufähige Theorie der Diskriminantenfunktionen. 
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Bei der Schriftleitung sind folgende Bücher eingegangen (ausführiiche Besprechung bleibt vorbehalten) 
Dr, phil, V. Happach, Ausgleichsrechnung. Dr.-Ing. E. Schmidt (o. Prof. a. d. Techn. Hoch- 3 


Ein Lehrbuch der Fehlerausgleichung nach der Me- 
thode der kleinsten Quadrate in Wissenschaft und 
Technik. (Teubners Mathematische Leitfäden, Bd.18.) 
2. Aufl., 1048. mit 26 Abb. Leipzig 1950. B.G. 
Teubner Verlagsgesellschaft. Preis kart. 6,20 DM. 


Prof, Dr. Richard Gans, Vektoranalysis 
mitAnwendungaufPhysikundTech- 
nik. (Teubners Mathematische Leitfäden Bd. 16.) 
7. Auflg 119 S. mit 44 Abb. Leipzig 1950. B.G. 
Teubner Verlagsgesellschaft. Preis geb. 5,90 DM. 


Dr.-Ing. R. Zurmühl, Matrizen, eine Dar- 
stellung für Ingenieure. XV + 4278. mit 25 Abb. 
Berlin-Göttingen-Heidelberg 1950. Springer-Verlag. 
Preis geb. 25,50 DM. i 


schule Braunschweig), Einführung in die 
technische’Thermodynamik und in die 
Grundlagen der chemischen Thermodynamik. Vierte 
überarbeitete und erweiterte Auflage. XVI-+ 520 8. 
mit 244 Abb., 69 Tabellen und 3 Tafeln. Berlin- 
Göttingen-Heidelberg 1950. Springer-Verlag. Preis 
geb. 30,— DM. EL‘ 


H. Rödel (Studienrat a. d. Ing.-Schule Darmstadt), 
Mechanik für Ingenieure Bd. IV Dy- 
namik einschl. Schwingungslehre., 
(Westermanns Fachbücher der Ingenieurkunde) 198 
S. mit 242 Abb. Braunschweig-Berlun-Hamburg 1949. 
Georg Westermann-Verlag. Preis geb. 10,40 DM, 
kart. 9,20 DM. - j 


NACHRICHTEN 


25 Jahre Göttinger Institut für Strömungsforschung. 


Am 15. Juli 1950 erinnerte das Max-Planck-Institut 
für Strömungsforschung in einer schlichten Feier 
Schüler, Freunde und Institutsangehörige an seine 
vor 25 Jahren erfolgte Gründung als Kaiser-Wilhelm- 
Institut für Strömungsforschung. Das Göttinger In- 
‚stitut, das seine Tätıgkeit 1908 als Modellversuchs- 
anstalt der „Motorluftschiff-Studiengesellschaft‘‘ be- 
gann, verdankt seine Entstehung dem Wirken Felix 
Kleıns, des berühmten Organısatoıs der Göttinger 
angewandten Wissenschaften. Sein Weitblick war es, 
der die Wichtigkeit der Strömungsforschung in der 
damaligen Luftfahrttechnik, eıkannte und den Ideen 
Pref, Ludwig Prandtls zur praktischen Verwirk- 
lichung verhalf. Schon wenige Jahre später wurde der 


erste Plan eines größeren Instituts gefaßt, welches 


nicht nur Modellversuchen im Windkanalt), sondern 
in der Hauptsache der wissenschaftlichen Erforschung 
der hydroaynamischen Grundlagen im weitesten Sinne 
dienen sollte, Infolge‘ der Zeitverhältnisse gelang es 
jedoch erst 1923/24, nachdem Prof. Prandtl eınen 
Ruf an die Techn. Hochschule München erhalten hatte, 
diese Pläne wieder aufzunehmen und die Mittel zur 
Gründung eines ‚„Kaiser-Wilhelm-Instituts‘“ für Strö- 
mungsfoischung‘‘ aufzubringen. So konnte,das neue 
Haus schließlich am 16. Julı 1925 eingeweih® werden. 
In seiner Einweihungsrede sagte Prof. Prandtl nach 
einem Überblick über die neuen Versuchseinrichtungen 
am Schluß: ‚Ich erwarte von dem neuen Institut 
weniger aufsehenerregende Entdeckungen als vielmehr 


eine solide systematische Durchforschung des ganzen - 


:) Zu dem Windkanal von 1908 war während des Krieges 


/, an der Techn. Hochschule München ernannt. 


ein zweiter größerer hinzugekommen, 
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Gebietes“. Nun, die damit gewonnenen Erkenntnisse 
sind ein Teil des Lebenswerkes Ludwig Prandtls 
geworden. Und gerade die „solide systematische 
Durchführung‘, dıe die Göttinger Schule betrieb, 
erregte das Aufsehen der Fachgenossen in aller Welt. 
Hunderte von Veröffentlichungen geben Zeugnis 
davon, 

Die Feier wurde durch Begrüßungsworte des 
jetzigen Leiters Pıof. Albert Betz eıngeleitet. In 
seıncı Festiede gab Pıof. Prandtl eınen Bericht über 


die Entwicklung des Instituts und geaachte dabeı ım 


einzelnen der Leistungen seiner Schüler und Mıt- 
arbeiter. Der Präsident der Max-Planck-Gesellschaft, 
Prof. Otto Hahn, verkündigte dann, daß der wıssen- 


schaftliche Rat der Max-Planck-Gesellschaftt die 


Herren Prof. Ackeret, Zürich und Pıof. Busemann, 
USA., zu Auswärtigen WissenschaftlichenMitglieden, 
die Herıen Pıof. Tollmien und Prof. Vogelpohl, 


- beidein Göttingen, zu Wıssenschaftlichen Mitgliedern r 


des Max-Planck-Instituts für Strömungsforschung und 
damit auch der Max-Planck-Gesellschaft eınannt 
habe. Anschließend fanden eine Institutsbesichtigung 
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